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FOREWORD 

This  report  documents  the  effort  accomplished  at  Rocketdyne 
under  Contract  No.  F04611-67-C-0103  for  the  period  15  May  1967 
to  19  August  1968.  The  program  was  sponsored  by  the  Air  Force 
Rocket  Propulsion  Laboratory,  Research  and  Technology  Division, 
Edwards,  California,  Air  Force  Systems  Command.  United  States 
Air  Force.  Air  Force  Project  Engineers  were  Mr.  Donald  Saylak 
anc  Captain  Scott  Beckwith. 

The  program  to  study  the  applicability  of  viscoelastic 
fracture  equations  to  the  problem  of  cumulative  damage  in  solid 
rocket  propellants  was  a  portion  of  a  broader  study  of  cumulative 
damage  in  which  the  Air  Force  has  sponsored  several  parallel 
approaches  by  separate  contractors.  It  is  expected  that  the 
reports,  taken  together,  will  provide  new  insight  into  the  ap¬ 
plicability  and  appropriateness  of  the  several  techniques.  In 
addition  to  this  program,  Lockheed  Propulsion  Company  investi¬ 
gated  the  volumetric  response  mechanisms  of  the  propellant  and 
have  reported  their  work  (as  well  as  the  propellant  characteriza¬ 
tion  data)  in  AFRPL-TR-68-130.  Aerojet-General  Corporation 
(Sacramento)  investigated  various  damage  mechanisms  using  a 
linear  cumulative  damage  approach.  Their  work  is  reported  in 

AFRPL-TR-68-131. 

The  consulting  services  of  Dr.  M.  L.  Williams,  University 
of  Utah;  Dr.  P.  J.  Blatz,  North  American  Rockwell  Science  Center; 
and  Mr.  W.  B.  Jones,  also  of  the  University  of  Utah,  ore  grate¬ 
fully  achknowledged.  Their  intuitions  and  judgments,  constrained 
by  technical  understanding  gained  after  long  experience  with 
rocket  propellants,  have  provided  the  ideas  for  the  bulk  of  this 
program. 

This  technical  report  has  been  reviewed  and  is  approved. 


Donald  Saylak  (RPMCB) 
Project  Engineer 
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Scott  W.  Beckwith  (RPMCB) 
Captain,  USAF 
Project  Engineer 
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ABSTRACT 

This  report  documents  a  12-month  program  con¬ 
ducted  by  Rocketdyne  to  study  cumulative  structural 
damage  in  solid  rocket  propellant.  Tests  conducted 
and  expei'imental  data  obtained  therefrom  are  de¬ 
scribed,  and  conclusions  drawn  from  the  study  are 
presented.  The  theory  of  fracture  in  viscoelastic 
materials  provides  the  analytical  expressions  for 
evaluating  damage  to  propellant  during  loading  his¬ 
tories  comparable  to  those  which  might  be  experi¬ 
enced  by  rocket  motor  grains.  Poker  chip  and  biaxial 
8 trip  tensile  tests  were  used  to  simulate  the  states 
of  stress  typically  occurring  in  actual  motors. 

Simple  constant  rate  tests  were  used  to  evaluate  the 
constants  necessary  to  apply  the  fracture  equations 
to  histories  involving  sequences  of  loadings. 
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SECTION  I 
INTRODUCTION 

This  report  describes  the  experimental  work,  the  analytical  calcu¬ 
lations,  and  the  resulting  conclusions  of  the  portion  of  the  Cumulative 
Damage  Study  Program  that  was  performed  at  Rocketdyne.  The  program  was 
funded  by  the  Air  Force  Rocket  Propulsion  Laboratory  on  Contract  No. 
F04611-67-C-0103. 

The  program  was  designed  specifically  to  test  the  applicability  of 
the  recently  developed  viscoelastic  fracture  techniques  to  the  predic¬ 
tion  of  crack-type  propellant  failures  following  loading  histories  char¬ 
acteristically  experienced  by  rocket  motor  grains.  Such  histories  usu¬ 
ally  involve  sequences  of  loadings,  each  differing  from  the  other,  until 
the  motor  either  satisfactorily  performs  its  mission  or  prematurely  fails. 
Current  failure  analyses,  however,  are  only  valid  following  elementary 
loadings  such  as  constant  loads  or  constant  loading  rates  or  constant 
amplitude  and  frequency  vibrations.  The  customary  procedure  is  to  select 
the  dominant  loading  and  analyze  the  structural  capability  of  the  grain 
following  its  application.  Then,  judgment  is  used  to  modify  the  solution 
to  account  for  damages  incurred  by  earlier  loadings.  It  was  felt  that 
the  viscoelastic  fracture  techniques  might  well  provide  an  analytical 
vehicle  for  making  engineering  evaluations  of  the  effect  of  the  earlier 
loadings . 

Accordingly,  the  experimental  test  phase  of  this  program  was  oriented 
to  evaluate  this  technique.  It  should  be  noted,  however,  that  all  test 
results  showing  sensitivity  to  damage  accumulation  are  invaluable  for  the 
study  of  the  applicability  of  any  analysis  procedure  and  consequently,  are 
reported  in  descriptive  detail. 

SUMMARY  OF  THE  THEORY 

A  theory  of  fracture  applicable  to  linear  viscoelastic  materials  has 
been  formulated  by  Williams  and  his  coworkers  (l,  2,  3)»*  The  approach 
utilizes  an  energy 'balance  as  originally  proposed  and  demonstrated  by 
Griffith  (4)  for  brittle  fracture  in  elastic  materials. 

Griffith  postulated  that  fracture  always  originates  at  defects  or 
flaws  which  exist  throughout  real  materials.  This  explanation  readily 
accounted  for  the  fact  that  glass  and  similar  brittle  crystalline  mate¬ 
rials  always  fail  at  stresses  much  smaller  than  those  known  to  be  necessary 

*Numbers  in  parentheses,  i.e.,  (l,  2,  3)>  refer  to  the  tabulated  refer¬ 
ences  on  page  235- 
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to  rupture  the  interatomic  or  intermoleeular  bonds.  He  was  successful 
in  showing  that  the  rupture  stress  in  glass  increased  in  proportion  to 
the  inverse  of  the  existing  flaw  size. 

Griffith's  approach  utilized  a  mathematical  solution  derived  by 
Inglis  (5)  for  the  stresses  and  strain  energy  in  an  initially  cracked 
thin  plate.  The  initial  flaw  located  at  the  center  of  an  infinite  plate 
was  assumed  to  be  ellipiically  shaped  but  with  the  major  axis  taken  to 
be  very  long  in  comparison  with  the  minor  axis,  giving  the  effect  of  a 
line  crack.  This  specimen  is  illustrated  in  Fig.  1. 


o  a 


a 


Figure  1.  An  Infinite  Tensile  Specimen  With  An 
Initial  Elliptically  Shaped  Crack 

4ccording  to  Griffith's  plane  stress  analysis,  fracture  will  occur  if 
and  when 

(V  +  Ay)  =  0  (1) 

where  V  is  the  potential  energy  of  the  system  (see  (6,  p  382)  for  a  def¬ 
inition  of  potential  energy),  A  is  the  free  surface  area  of  the  specimen 
which  increases  as  the  crack  grows,  and  Y  is  the  free  surface  energy  per 
unit  of  new  surface  area.  The  form  taken  by  the  terms  of  Eq.  1  depend 
upon  the  boundary  conditions.  For  the  important  case  vhere  the  applied 
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stress,  ( t ,  (some  distance  away  from  the  crack)  is  held  constant  and  the 
crack  length  is  2c 


d_ 

dc 


(-  -  4Yo) 


<  0 


(2) 


Here  the  first  term  reflects  the  lowering  of  the  strain  energy  due  to 
the  increase  in  the  crack  length  while  the  second  reflects  the  increase 
in  surface  energy  due  to  the  increase  in  the  free  surface  area.  A  plot 
which  schematically  shows  these  two  terms  is  given  in  Fig.  2.  When  the 
derivative  is  evaluated,  one  finds  the  most  common  and  well  known  form 
of  the  Griffith  criterion 


a 

cr 


(3) 


Figure  2.  Energy  Criterion  For  a  Plate  With  a  Crack 
of  Length  2c 

The  suggestion  that  the  energy  balance  equations  might  prove  valuable 
for  predicting  the  onset  of  fracture  in  elastomeric  materials  has  been 
studied  and  tested  extensively.  A  brief  survey  of  the  monograph  edited 
by  Rosen  (7)  on  the  subject  reveals  the  extent  of  the  variety  and  success 
of  these  explorations . 
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During  the  1950's,  exploratory  efforts  to  apply  the  energy  balance 
approach  to  amorphous  polymers  were  jr ' necree  by  RivVn  and  Thomas  (8). 
Their  studies  clearly  indicated  the  exigence  ox  «  characteristic  energy 
for  the  initiation  of  crack  growth  defined  {.}«£•».  i:.  Eq.  l)  as: 


where  F  is  the  total  stored  strain  energy  t  •  -o‘"icr:'.pt  u  indicates 
that  the  derivative  is  evaluated  with  the  '>r,ru  317  «Ut  pin  some?  ;a  held 
constant.  But,  taking  the  lead  established  by  Ijvi,:  Orowan  (9,  10), 

Rivlin  and  Thomas  emphasized  the  point  that  y  miwt  cub,  for  these  mate¬ 
rials,  necessarily  be  interpreted  as  a  free  surface  o.  ergy. 

Later  Greensmith  and  Thomas  (ll)  hypothesize-..  as  ■  result  of  exten¬ 
sive  experimentation,  that  the  characteristic  energy  depended 

st'  *».  gly  upon  both  the  temperature  and  the  tiirn.  s.  ,e  associated  with  the 
st'fa.ing  history.  However,  constant  moduli  (in  time)  were  used  in  the 
calculations  supporting  this  conclusion. 

The  generalization  of  the  energy  balance  concept  to  the  prediction 
of  crack  type  fai?wes  in  viscoelastic  materials  is  possible  primarily 
tocotv o  sim*?.'*ied  flaw  geometries  are  hypothesized.  The  assumption  that 
the  flaws  are  either  spherical  or  cylindrical  in  shape,  centered  in  a 
spherical  or  cylindrical  body  of  incompressible  material  has  been  made 
to  find  a  tractable  expression  for  incipient,  flaw  growth.  Because  of 
the  time-dependent  nature  of  the  material  properties,  a  balance  of  the 
time  rate  of  change  of  the  energy  (power)  rather  than  a  simple  energy 
balance  is  necessary.  This  balance  leads  to  criticality  expressions  quite 
similar  in  nature  to  Eq.  3  except  that  viscoelastic  material  properties 
replace  the  constant  Young's  modulus.  The  form  with  which  these  time- 
varying  properties  enter  the  equation  is  dependent  upon  the  loading  his¬ 
tory,  which  allows  the  balance  theory  to  be  used  in  evaluating  cumulative 
damage.  The  time  to  failure  following  essentially  any  loading  history, 
including  sequences  of  completely  different  types  of  loadings  (such  as  a 
step  strain  followed  by  a  superposed  vibratory  load)  can,  therefore,  be 
conveniently  predicted. 

The  formulation  of  an  energy  balance  criterion  with  expressions 
reflecting  viscoelastic  material  properties  is  especially  attractive  be¬ 
cause  it  provides  the  expressions  necessary  to  evaluate  the  influence  of 
the  time  or  rate  of  loading,  the  order  of  loading,  the  material  proper¬ 
ties,  the  state  of  stress,  and  perhaps  even  temperature  on  crack  initia¬ 
tion  in  a  viscoelastic  material.  This  unifying  idea  motivated  the  Rocket- 
dyne  approach  to  cumulative  damage. 
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Derivations  of  the  analytical  expressions  for  criticality  following 
practical  loading  histories  are  shown  in  Section  II.  Those  used  in  this 
report  are  summarized.  These  criticality  equations,  pp  30  and  33  present 
the  unified  theory  of  failure  for  viscoelastic  materials. 

Dr.  Paul  Blatz  of  the  North  American  Rockwell  Science  Center,  has 
conducted  supplementary  research  activity  while  serving  as  a  consultant 
to  this  program.  A  portion  of  his  efforts  is  reported  in  Appendix  I. 

These  studies  have  been  directed  toward  determining  the  importance  of 
finite  strains  relative  to  the  interpretation  of  the  fracture  results, 
interpretation  of  the  triaxial  poker  chip  test  data  after  the  many  growing 
internal  flaws  have  effectively  rendered  the  specimen  into  a  sponge,  and 
whether  or  not  the  fundamental  assumptions  on  which  the  flaw  growth  ex¬ 
pressions  are  based  can  be  modified  to  improve  or  refine  the  resulting 
predictions. 

As  the  applicability  of  a  fracture  theory  is  widened,  either  by 
providing  a  more  fundamental  basis  or  by  experimentally  extending  the 
range,  its  value  to  the  practicing  stress  analyst  is  likewise  enhanced. 

So  a  fracture  criterion  capable  of  evaluating  many  if  not  all  of  the  in¬ 
fluences  on  the  initiation  of  a  propellant  grain  crack  is  a  very  impor¬ 
tant  contribution  to  the  rocket  motor  industry. 

SUMMARY  OP  THE  TESTING 

The  laboratory  experimental  investigations  were  designed  specifi¬ 
cally  to  test  the  applicability  of  the  spherical  and  cylindrical  flaw 
theories  to  the  prediction  of  cracking  in  propellants.  Those  parameters 
which  were  believed  to  be  the  most  influential  were  selected  to  vary  as 
extensively  as  possible.  The  tests  provided  a  wide  variation  ic  loading 
rates  and  times,  the  state  of  stress,  and  the  various  kinds  of  loadings, 
including  temperature.  The  propellant  utilized  for  the  tests  was  the 
STV  propellant  furnished  by  Lockheed  Propulsion  Company.  It  was  an  86.5/6 
solids-loaded  CTPB  formulation  typical  of  the  state  of  the  art.  A  sum¬ 
mary  of  the  testing  program  presented  in  Table  I. 

The  tests  shown  in  the  table  fall  into  two  categories:  first,  the 
single  environment  or  elemental  tests  and,  second,  multiple  environment 
tests.  The  elemental  tests  (conducted  early  in  the  program)  were  designed 
specifically  to  provide  the  'ata  necessary  for  a  numerical  evaluation  of 
the  arbitrary  constants  required  by  the  analytics.  Following  these  ex¬ 
periments,  multiple  loading  tests  were  conducted  to  check  the  feasibility 
of  using  the  theory  to  predict  cracking  following  practical  loading  his¬ 
tories.  The  consistencies  between  the  various  elemental  tests  also  pro¬ 
vided  a  check  on  the  reliability  and  applicability  of  the  energy  balance 
expressions  to  propellant  cracking. 


5 


^c<ryri 


TABLE  I.  Continued 
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TABLE  I.  Concluded 
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The  triaxial  tests  of  the  program  which  are  most  applicable  to  the 
evaluation  of  the  spherical  flaw  portion  of  the  theory  were  conducted 
using  the  poker  chip  specimens.  These  circular  thin  flat  specimens, 
bonded  between  two  platens,  were  pulled  at  constant  crosshead  rates  in 
Instron  machines.  Load  cells  were  installed  at  the  center  of  one  platen 
in  each  test  and  were  used  to  monitor  the  center  bond  stresseg. 

As  shown  in  Table  I,  16  of  these  tests  were  conducted,  some  at  an 
intermediate  crosshead  rate  and  room  temperature  and  others  across  a 
spectrum  of  rates  and  temperatures.  Forces,  platen  displacements,  and 
center  load  cell  outputs  were  recorded  during  all  of  the  tests.  The 
center  load  cell  output  provided  a  sensitive  measure  of  the  time  of 
failure  during  the  tests.  In  addition  to  the  constant  3tress  rate  tests 
of  the  poker  chip,  .three  constant  load  tests  of  the  poker  chip  specimens 
were  also  conducted.  These  three  tests  represent  multiple  environmental 
conditions  in  that  they  essentially  had  the  load  applied  instantaneously 
and  then  held  constant  for  the  duration  of  the  test.  This  provided  a 
limiting  case  of  a  bilinear  stress  rate  type  loading  history.  A  signi¬ 
ficant  amount  of  time  elapsed  from  the  time  the  load  was  applied  until 
failure  occurred.  It  was  possible  to  X-ray  the  cross  section  (X-ray  ex¬ 
posures  take  up  to  5  minutes  on  the  poker  chip)  to  observe  the  internal 
mechanisms.  These  X-rays  were  most  informative;  they  enabled  a  sound 
technical  explanation  of  the  constant  load  test  results.  In  addition, 
they  provided  insight  into  the  meaning  of  the  center  load  cell  readings 
during  constant  stress  rate  tests. 

Because  cracks  usually  originate  in  the  center  ports  of  rocket  grains 
and  because  these  center  ports  are  frequently  in  essentially  the  same 
state  of  stress  as  the  biaxial  strip  tensile  specimen,  major  attention 
was  given  to  tests  of  this  nature.  Constant  rate  tests  of  biaxial  strip 
specimens  were  run  during  the  propellant  characterization  by  Lockheed. 

The  data  from  these  tests,  both  failure  stress  and  failure  strains,  were 
used  as  elemental  tests  in  the  Rocketdyne  program. 

In  addition,  vibration  fatigue  tests  witu  both  displacement-  and 
force-specified  boundaries  were  conducted  using  the  new  Cumulative  Damage 
tester  at  Rocketdyne.  This  testing  machine  provided  the  capability  of 
imposing  sinusoidal  platen  displacements  on  the  biaxial  strip  specimen. 

The  capability  of  pre-programming  any  required  thermal  history  (-75  to 
250  F)  was  provided  with  an  environmental  box  surrounding  the  specimen. 
Other  testing  capabilities  of  the  machine  included  bi-linear  rates 
(strain  or  stress),  combinations  of  fatigue  tests  involving  different 
frequencies  or  different  amplitudes,  and  the  very  practical  test  involv¬ 
ing  simultaneous  changes  in  strains  and  temperatures.  As  shown  in  Table  I, 
most  of  these  capabilities  were  exploited  on  the  program. 
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A  sequence  of  four  tests  involving  small,  hollow  propellant  cylin¬ 
ders  were  performed  by  Hilzinger  at  Rocketdyne  Research  (Canoga  Park). 
Photomicrographs  showing  the  behavior  of  binder  and  failure  during  these 
constant  pressure  rate  tests  were  taken.  Close-up  movies  were  also  made. 
These  pictures  provided  a  basic  understanding  of  the  propellant  behavior 
before  and  during  the  actual  fracturing  processes. 

Appendix  II  summarizes  work  performed  by  J.  E.  Hilzinger.  He  was 
able  to  capture  on  film  the  initiation  and  propagation  of  cracks  in  the 
STV  propellant  following  a  constant  stress  rate  loading  history.  The 
specimens  used  for  these  studies  were  small,  hollow  cylinders  subjected 
to  internal  pressurization.  Perhaps  the  most  significant  finding  was  the 
basic  difference  in  the  behavior  of  the  specimens  with  cast  surfaces 
from  that  of  specimens  with  machined  surfaces. 
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SECTION  II 


THEORY  OF  VISCOELASTIC  FRACTURE 

Until  about  I960,  the  primary  contribution  to  the  study  of  failure 
in  amorphous  viscoelastic  media  (e.g.,  polymers)  was  the  work  of  Rivlin 
and  Thomas  (8)  which  extended  the  Griffith  fracture  hypothesis  (see 
Section  i)  to  the  rupture  of  rubber.  The  rupture  (or  tearing)  was  shown 
to  result  from  strain  energy  being  released  at  a  greater  rate  than  it  can 
be  transferred  into  surface  free  energy  or  other  dissipative  mechanisms. 
Their  results  showed  the  change  in  strain  energy,  F,  with  respect  to  the 
area  of  a  cut  in  a  thin  sheet,  A,  computed  with  fixed  displacement  over 
the  boundary  where  forces  are  prescribed  is  given  by 

(&F/&A)  |  u  =  y  (5) 

where  y  is  the  characteristic  energy  of  tearing.  Although  the  form  of 
the  result  is  similar  to  Griffith's  criterion,  y  includes  other  dissipa¬ 
tive  mechanisms  than  surface  free  energy.  In  particular,  at  long-time 
rubbery  and  short-term  glassy  conditions,  surface  energy  effects  will  be 
predominant.  In  the  transition  region,  a  significant  amount  of  dissipa¬ 
tion  and  resultant  heat  flow  will  absorb  the  energy.  These  phenomena 
qualitatively  agree  with  the  observations  of  Greensmith  and  Thomas  (ll) 
that  the  characteristic  tearing  energy  is  strongly  dependent  on  the 
tearing  rate  and  temperature. 

In  a  review  article  (l)  in  1963,  Williams  related  the  findings  to 
date  to  include  evaluation  of  the  effects  of  load  and  crack  tip  geometry 
at  fracture  initiation,  the  prediction  of  crack  acceleration  and  stable 
propagation,  and  the  results  observed  in  triaxial  stress  fields.  Further, 
he  discussed  an  approach  to  cumulative  damage  based  on  superposition 
analogous  to  Miner's  hypothesis  for  fatigue  or  vibration  in  metals. 

Williams  later  proposed  (2)  a  qualitative  expression  of  the  Griffith 
critical  stress  formulation  as 

°crit  =  k  (\  +  Yd  <■  \  +  ...  (6) 

where  E  is  the  material  modulus;  l,  the  crack  length;  and  the  subscripts 
b.  d,  v,  .  .  .  stand  for  the  brittle,  ductile,  and  viscoelastic  dissipa¬ 
tive  processes,  respectively.  Pursuing  the  application  of  the  thermo¬ 
dynamic  theory  of  fracture  presented  earlier  (l),  equations  of  criticality 
were  established  for  general  conditions  of  loading  and  specimen  geometry. 
The  following  development  of  criticality  formulas  is  based  on  this  ap¬ 
proach. 
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THE  ENERGY  APPROACH  TO  FAILURE 

The  elemental  nature  of  the  Griffith  approach  to  failure  is  intrin¬ 
sically  attractive  since  it  overcomes  the  indeterminancy  of  the  infinite 
strain-stress  state  at  the  crack  tip  (i.e.,  an  infinitesimal  area).  Of 
further  benefit  is  the  fact  that  energy  is  a  scalar  quantity  and  may  be 
summed  directly  as  in  Eq.  6.  Considering  the  case  of  the  viscoelastic 
materiel  and  neglecting  kinetic  energy,  which  will  be  small  for  slow 
velocities  of  crack  propagation,  one  has 

I  =  F  +  2D  +  SE  (7) 

• 

where  I  is  the  rate  at  which  work  is  being  done  on  the  body  (power  input), 
F  is  the  rate  of  increase  of  strain  energy,  2D  is  the  viscous  dissipation 
(work  converted  to  heat);  and  SE  is  the  rate  of  increase  of  the  surface 
energy. 


Neglecting  work  done  by  the  body  forces,  and  assuming  a  volume  con¬ 
stant  body  of  volume  t  and  surface  S,  (6) 


(I) 

dt 


i.  ♦ «.  **V 


where  T.  are  the  surface  forces  prescribed  over  the  surface  and  the 
Uj  are  the  displacements  prescribed  over  the  remaining  surface  Sg. 


(a)  ■  Sr  f 


where  Y  is  the  characteristic  tearing  energy.  Finally,  due  to  the  viscous- 
elastic  nature  of  the  material,  the  remainder  of  the  power  input,  typ¬ 
ically  denoted  as  free  energy  must  be  accounted  for  in  terms  of  elastic 
strain  energy  and  viscous  dissipation.  Thus 


(F)  ■  St  f  °ieidtdT  - 


vhere  P  is  the  stored  strain  energy  and 


I  f  w 


't  0 


is  the  volume  integral  of  the  strain  energy  density.  Using  this  formu¬ 
lation  of  xhe  energy  approach,  it  is  now  pertinent  to  inquire  into  the 
conditions  for  which  I  exceeds  F  +  2D  and  resulting  flaw  growth  occurs. 
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ELASTIC  ANALYSIS  FOR  VISCOELASTIC  MEDIA 


The  problem  of  analytical  determination  of  the  stress-strain  state 
in  even  a  simple  geometry  is  complicated  by  the  time-temperature  depen¬ 
dence  of  the  material  properties  when  the  material  is  of  a  viscoelastic 
nature.  Historically,  two  approaches  to  the  solution  of  viscoelastic 
problems  are  available.  A  thorough  discussion  of  the  representation  of 
the  material  properties  is  presented  in  (12).  Of  interest  here  for  sub¬ 
sequent  application  is  the  analogy  between  viscoelastic  and  elastic 
response . 

In  principle,  with  zero  initial  conditions,  the  Laplace  or  Fourier 
time-transformed  linear  viscoelastic  field  equations  and  boundary  con¬ 
ditions  are  equivalent  with  the  equations  of  an  elastic  body  with  the 
same  geometry.  This  analogy,  essential  to  the  pursuit  of  viscoelastic 
stress  analysis,  was  expressed  by  Lee  (13)  for  isotropic  materials  and 
Biot  (14)  for  anisotropic  materials.  Of  special  importance  to  the  appli¬ 
cation  of  this  concept  are  the  requirements  for  linear  viscoelasticity, 
special  stress  distribution  independent  of  material  properties,  and 
separation  of  space  and  time  dependence  of  prescribed  loads  and  displace¬ 
ments.  Byway  of  illustration,  take  the  case  of  an  elastic  uniaxial 
specimen  loaded  at  constant  rate  of  displacement 

1  -  1  (11) 

Now,  consider  a  constant  strain  test  of  a  viscoelastic  uniaxial  specimen 

O  , (t)  =  E  ,(t)  s  (12) 

rel  rel'  7  0  ' 

Applying  the  Laplace  transform, 

a  "  (p)  =  E  .  (p)  ~  (13) 

rel'r/  rel'r/  0 

or 

^  <u) 

Next,  consider  the  constant  strain  rate  test  (Rate  =  R)  on  the  visco¬ 

elastic  specimen.  The  ratio  of  the  transform  of  the  strain  input  to  that 
for  the  constant  strain  test  is 
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Prom  the  thermodjuamic  considerations  of  linear  viscoelastic  materials 
(15)  one  obtains 


tens 


6tens  R 


or 


"  Er.l<p> 


Patens 

R 


(16) 


which  may  be  written  as 


°teD>>  *  Erel<p>  f  =  Erel(p)  p  e(p) 


(  7) 


which  is  the  transform  of  the  viscoelastic  stress  response  to  the  con¬ 
stant  strain  rate  input.  To  complete  the  solution, 


16 


SV/ifc 
'  «V*A'W-4i 

'■’  <-M 


J 5  K Ufvrtl -fi  $,  ‘vV>K'»<ir  * 


It  is  worthwhile  to  note  that  the  elastic  analysis  is  valid  to  the  point 
of  obtaining  a  transformed  solution.  The  final  step,  inversion  of  the 
transformed  solution  to  the  real-time  plane,  is  generally  not  directly 
available  and  must  be  accomplished  through  approximate  methods  as  de¬ 
scribed  in  (l6).  In  the  present  analysis,  direct  inverse  transformation 
is  adequate . 


A  final  area  of  note  in  the  treatment  of  viscoelastic  materials  is 
that  of  time-temperature  interdependence.  From  qualitative  consideration 
of  the  molecular  behavior  of  amorphous  polymers  (17),  it  is  possible  to 
infer  an  equivalence  relation  on  the  time-temperature  mode  of  material 
response.  From  empirical  information  obtained  from  the  reduction  of 
relaxation  data,  Tobolsky,  (l8)  suggested  that  those  data  could  be  ade¬ 
quately  represented  by 


l0g  Erel^  "  l0g  Ee 


log  Eg  -  log  Ee 


1  -  erf  (h  log  ) 


where  K  (t)  is  the  value  of  f(T)  at  which  the  reduced  time  parameter, 
t/f  (t)  equals  unity  and  h  is  standard  deviation.  Plotting  K  (t) 

7* 

against  T  revealed  an  inflection  point  at  the  temperature  T^,  usually 

about  10  F  above  the  glass  transition  temperature.  A  useful  approximation 
of  Tobolsky’ s  result  is  given  by 


-16  (T  -  Tq) 


100  +  T  -  Tt 


for  T  in  degrees  Farenheit.  A  different  representation  was  given  by 
Williams,  Landel  and  Ferry  (19)  which  resulted  in  placing  the  glassy 
relaxation  behavior  near  the  unit  reduced  time  t/a_,  =  1.  Subsequent 
trends  in  application  have  shown  that  it  is  desirable  to  determine  the 
shift  factor  a^,  such  that  a,p  =  1  occurs  at  77  F,  the  point  at  which  most 

laboratory  data  are  obtained.  The  form  of  the  equivalence  remains 


-Cj  (T  -  Ts) 
!°g  aT  -  +  (T  _ 


where  Tg  is  the  temperature  at  which  a^,  =  1.  Using  Eq.  21  it  is  possible 

to  handle  temperature  dependency  as  a  time  function  in  viscoelastic 
analysis . 


VISCOELASTIC  FRACTURE  ANALYSIS 


Equipped  with  the  energy  concept  and  the  elastic-viscoelastic  nry 

(and  keeping  in  mind  the  limitations  and  assumptions)  the  analysii 
fracture  initiation  in  viscoelastic  bodies  will  now  be  considered, 
this  point,  the  analysis  will  be  further  limited  to  infinitesimal  .  n 
theory. 

Three  principal  geometries  of  interest  are  illustrated  in  Fig.  3- 
The  spherical  flaw  subjected  to  uniform  loading  at  infinity  is  at  once 
a  simple  geometry  for  analysis  as  well  as  one  of  practical  interest  in 
the  analysis  of  void  growth  in  triaxial  stress  fields.  The  cylindrical 
flaw  geometries  in  both  plane  strain  and  plane  stress  fields  are  analogous 
to  rocket  motor  geometries  as  well  as  laboratory  test  specimens.  The 
subsequent  derivations  follow  the  methods  reported  in  (2  and  20). 

STRESS  BOUNDARIES 

Consider  a  long,  circular  cylinder  with  a  concentric  circular  flaw. 

The  material  is  incompressible,  linearly  viscoelastic,  and  will  respond 

with  infinitesimal  strain  to  the  prescribed  loadings  and,  e  =0. 

z 


For  a  boundary  condition  specified  by  o(b,t) 

u(b't)  "a  „TjL  i 


bE 


2  2 
a *  b“ 


Oof(t)  one  has 


(22) 


Since 


r  u(r,t)  =  b  u(b,t) 


u(r,t)  =  J 


v(t) 


rE 

1_  Jj 
2  ~  ,2 

a  b 

(23) 


From  elastic  theory, 
or(r.t).a 


1  - 


r2  , 


(24) 
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a.  Spherical  Flaw  Geometry 


Uniform  Boundary  Stress,  a f(t) 


Uniform  Boundary  Displacement, 

*  «(t) 


b.  Cylindrical  Flaw  Geometry;  Plane  Strain 


Uniform  Boundary  Stress, 
*  *<*> 


Uniform  Boundary  Displacement, 
«.  i(*) 


Cylindrical  Flaw  Geometry;  Plane  Stress 


Uniform  Boundary  Stress, 
°  *(*) 


Uniform  Boundary  Displacement, 
uo  *(*) 


Figure  3.  Initial  Flaw  Geometries 
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Farther  definiug 


,.\  12  a  -  CLa  2f(t) 

*H  -ni  k  -r  ;  R(t)  •  -  2  “~  ' 

*•  o  J  l  -  ot 


one  has 


"r  -  °e  ■  ^2  R<*>  & 

Yom  Eq.  V  and  considering  geometry  and  boundary  conditions  we 
note  v/iat 

e 

1  =  2"b  0.(b,t)  u  (b,t)  (32 


F  +  29 


rr 


2Qrrr  dtdr 


where  Q  =  stored  strain  energy  per  unit  volume  or 

Q  =  a  e  +  o  e  +  o  i 
r  r  0  0  z  z 


SE  =  2naY  (35) 

Tius,  for  an  elastic  body,  Eq.  7  becomes 

°  "  +  dt J^(t)  J  t  Qr  dtdr  "  b0r^b,t^  ”  (*»*)  (36) 

To  facilitate  the  solution  of  Eq.  36  for  the  viscoelastic  body,  we 
employ  the  Laplace  transform  analogy  so  that 

°(p)  ~  tffl(p)  =  — %  R  (p)  (37) 
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^ v— -j.  .  «*-,,  «♦*«*-  *  -**  -•<'*•  x  ■».’.'■?» wrsfr<tv 


and 


likeviae 


so 


er(p)  =  -  e  (p)  »  “§  p  Dcrp  (p)  R  (p) 
r 


U  (b,p)*--”  P  Dcrp  (p)  ft  (p) 


p  Q(p)  =  (y(p)  -  a0(p))  er(p) 


4a 


=  — 7-  P  D  (p)  r(p)‘ 

^4  r  crpvr/ 

For  obtaining  the  inverse  transform 
a 


■1  01 

p 


Q(p)  -  <  R  (p)  Jp  R  (p)  R  (p)J 


and 


■C  q(o-^*  r  Mof?  [v(5)  +c  d5 


(38) 


(39) 


(40) 

(41) 


(42) 


Then,  the  total  stored  energy  per  unit  length  is 

r  p 

-/.{')  Jo 


Qdt  2nrdr 


r .  .  n 

r\ 

l  iD  „(§  -  t) 

K>| 

“»h 
< _ 

|X  xs)tg 

y(?)  +Jo 

■iff^rr  *<T>  dTJ 

d§ 


(43) 
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And  differentiating  with  respect  to  time, 


y  +  2D  =  ■ 


-JjP  »(5>|s  v(5)  +  P  n 


(5-t) 


alT^r) 


R(t)  dr 


i  1  ^  f  ft  3D  (t  -  t) 

-  h  J  “<*>  t7[v(t)  d  -ift-'-T)"  r(t)  dTJ' 


Similarly,  for  the  power  input,  Eq.  32  gives 


(p)  »  -2rr  oo2  [”p  7  (p)  D_rp  (p)  R(p) 
L 


and  inverting 


I  >  2W  ( 
0 


»  f  ft  aD  ft  -  t) 

(*)  It  [Y(t)  +  Jo  iftr-77-  dT 


which  exactly  equals  the  second,  term  in  Eq.  44. 

Combining  Eq.  44,  46,  and  35, 

0  -  a  jv  -  ^§-  ft  R(?)  jf  Dg  »(s)  +  05  -3-  R(T)dTld? 


Up  to  initiation  of  flaw  growth,  a  =  0.  Thus,  the  condition  of 
criticality  is 


4o  2  rt  f  f 

rf rj,  xofe  »,««)*]. 


?  6D  (?  -  t) 

B(T)  dT  d? 


to  jm**^*^*.tays* ti*  w 


Further,  since  a  =  0,  it  is  possible  to  simplify  Eq.  48  to 

Y  *  }  ^7f_  20 '{Jo °  ,(?) «  > f(?)  +r  %^rr  f(T)  dT  d?) (49) 


3  0  -  k2)2  ^  "  2oo2  0  f(?)  j|  L_1  ^  Dcrp  (p)  7  (p) 


Equation  50  is  especially  useful  in  that  it  may  be  shown  that  the  right 
side  of  the  equation  is  independent  of  geometry  (within  the  three  cases 
being  considered)  and  the  term  multiplying  y/a0  is  a  function  only  of 
the  geometry.  The  use  of  the  inverse  transform  notation  will  be  shown 
to  be  practical  when  considering  various  prescribed  loading  functions. 

For  the  plane  stress  cylindrical  flaw  problem,  a  process  similar 
to  that  above  yields 

a  (1  -  K2)2  ^  -  2002  J0to  f<?)  L-1  |p  Dcrp  (p)  7  (P  j  d?  (51 ) 

while  for  the  spherical  flaw  problem  one  has 


If  (1  -  k 3)2  ^  -  2002 


DISPLACEMENT  BOUNDARIES 


j/  '(5>jfi-l|pDcrp<p)7(p)jdS 


We  will  now  follow  a  similar  process  to  that  described  for  the  stress 
boundary  to  a  plane  strain  cylinder  with  applied  displacements  at  the 
boundaries. 

As  before, 

u(r,t)  -  ^  u(b,t)  -  7  u  g(t)  (53) 


with  the  strains 


er(r,t)  -  -uog{t)  ~j 


en(r>t)  •  »0«(0  \  ■*€r(r>t) 
«  r 


The  principal  stress  difference  is  given  by 


4u  2 

Cr(r,t)  -  Og(r,t)  -  E  g(t)  b  -j 


”  "(f)  u0  E  «(*)  ~2 


Utilization  of  the  Laplace  transform  and  the  inversion  to  obtain  Q  results 


,  u  2b2  ft  aE  ,  (t  -  t) 

- f-v  ^  [*,  «<*> + Jo  rit'-'v)  •  i(t)  dT_ 


fb  A 
j«(t)  Jo 


Qdt  2nrdr 


4tu  2b2 
0 


-1 — 1-1  P  AtiSl  Te  *(?)  +  P 
z(t)  ba  j  ^  ®5  L  •  Jo 


3(S  -  t) 


g  (t)  dx  d? 


Differentiating  with  respect  to  time, 


2  2 
4nu  b 


F  t  2D  »  1 


^  r  ip  [v<!>  *  r  «<’>  *] « *[>  -  [v<*> 


*  I  li:  -■?}  * w 


• .  -vaixJtiBmiaiBaaSb 


The  rate  of  energy  input, 


I  »  2nbc 


2  2  r  C 

2T7bor(b,t)  a  (b,t)  -  —A.  i_ .  i_  iflii  /  E^g(t) 

+ L  «  (t)  *4 


which  again  is  equal  to  the  second  term  of  Eq.  59. 
And  finally, 


se  -  2ti *y 


Substituting  into  the  power  balance  equation, 


O  -  i  (  Y 


As  before,  up  to  flaw  growth  a  =  0  and  the  condition  for  criticality 
may  be  expressed  as 

2  6a  ‘  2^>  I  ^  V(s)  +  ^  t)  }«(t)  dT  d5  (63) 

The  right  side  of  Eq.  63  is  independent  of  the  geometry  whereas  the 
multiplier  of  y/a0  is  entirely  a  fraction  of  geometry.  Following  the 
methods  above  one  may  deduce  for  the  plane  stress  cylindrical  flaw 


,  <■  - kS,!  ;  ■ 2  £>2  j?  ^  V,!)  *  X 


!  «„.(S  -  t) 
3  15  -  t) 


and  for  the  spherical  flaw 


|  (1  -  k3)  *  -  2  (■ 


!  W„,(S  -  t)  ] 

Hs  -t)  dT 


.av.'>sn5S<B3B©*2 


Utilizing  the  expression  for  stress  and  displacement  boundaries  we  will 
now  consider  the  correlation  of  criticality  for  a  variety  of  loading 
conditions . 

CONDITIONS  OF  CRITICALITY  FOR  FLAW  GROWTH 
Stress  Boundaries 

Utilizing  the  expressions  of  Eq.  50  and  52,  one  may  deduce  the  "energy" 
available  for  creation  of  surface  for  stress  boundary  loadings.  It  is 
necessary  only  to  evaluate  the  right  hand  side  of  each  expression  to 
determine  the  "criticality"  for  all  three  geometries.  In  practice,  it  is 
necessary  to  determine  the  value  of  the  "energy"  as  a  function  of  time  and 
infer  the  time  to  failure  by  intersecting  the  resultant  curve  with  the 
le*t  side  of  the  expression  for  each  geometry. 

One  case  of  special  interest  is  that  of  the  step  input 

o(b,t)  -  ag  i(t)  (66) 

From  above,  we  wish  to  evaluate 


2ao2  Io°  f(?>35  L'^P  D  (p)  *  (p)| 


which  is 


P° '  ><«>  fi  [»,  k?)  -  P  ~f fer  1(t)  dT] d5 
-  2oo8  Jo°  S*  fe  d?  +P  1(5  -  T)dT]  d? 

%^d?]  *  <  -  »J 

=  0  2  L)  (t  )  -  -il 

0  L  crP  0  2  J 


Let  us  next  consider  a  bilinear  stress  function 


rVi1 

b,t)  -\oo[c,t  ♦ 


(C2  -  C,)  (t  -tj)J 


,  t  <  tj 
,  t  >  tj 


***  ■  -t  *. 
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From  Eq.  67  we  have 


2on2  jXs)  l'1  {^P  Dcrp(p)  7(p)J  -l  °  L-*^p  Dcrp(p)  f(P)j  d?  ^ 


and  for  the  loading,  Eq.  69 


’  cl  j[  D(r)dr  +  (c2  -  Cj)  Jt  D(t  -  tj)dT 


and  we  obtain 


a°oa  (cia « 11  »<*>*  L'1  -  cia  lh  r n(T)dT 
+  [el5  +  (c2  -  cl}^  -  *l>j  [C1  I  D<t>  +  <c2  -  ci)  It *  D{t  -  tl)dT]  |ti  • 

*  c2  Xj  °  C1  II  D(T)dT  ‘  c2(c2  ‘  el>  I j  °  X,?  D(t  '  (72- 

■  [Vo  +  <c2  *  Cl>(t.  *  V]  [C1  D<1)  (to>  +  <c2  *  cl>  d(1>  ^0  *  41>] 

-  clc2  d(2)  <4o>  *  CI  -  c2>  d(2)  <V  ’  ‘2  <c2  '  cl>  D(2)  <4o  ’  4I>) 

(n) 

where  the  D'  1  (tQ  -  tj)  notation  follows  that  demonstrated  by  Jones 
and  others  (2l).  For  the  simpler  case  of  a  ramp  function, 


o(b,t)  «  ooct 


merely  set  =  Cg  =  c  in  Eq.  72  and  obtain 


02  (c2  t/1'  (t0)  -  cV2)  (to)J 
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yA'  c*sr?  ,fft ,  iMl  -»  t.  »$<-<*; 


For  a  rczp-to-constant  loading,  typical  of  an  approximation  of  a  step 
loading  in  a  real  test,  set  =  0  in  Eq.  72  with  the  result 

2oo2  [ci\  [d(1)  (*.)  -  d(1)  (to  -  *1>]  -  ci2l)(2)  <v}  (75) 


where 


o(b,t) 


\ViY 


t  <  t 
t  >  t 


1 

1 


These  results  and  those  for  other  loadings  of  interest  are  summarized 
in  Table  II. 

Displacement  Boundaries 

The  evaluation  methods  used  above  are  now  to  be  applied  to  Eq.  63 — 
65  for  prescribed  displacement  boundaries.  The  equation  of  interest  is 


°Erel(T> 


or  in  alternate  notation 


4>2{f  ^  L_1  (P  Erel<P>  7(P)}  d? 

Consider  now  the  constant  displacement  boundary, 

u(b,t)  »  uQ  l(t),  i.e,  g(t)  =  1 


(76) 


(77) 


(78) 


which  gives  in  Eq.  76 

2(ir)2  (f  if1  V  {z)  d?  +fo  ^  f  -  T> dTd?} 

-  2<f >2 (a  Eg +  [E<*>  -  Eg] -  [E<*>  -  Eg] ■  (r>2 Eg 


(79) 
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TABLE  II.  CRITICALITY  CRITERIA  FOR  STRESS  LOADINGS 


Utilizing,  Eq.  77, 


(84) 

the  final  expression  being  the  one  of  greatest  applicability  in  reducing 
actual  tost  data  where  the  failure  displacement  and  ratio  of  constants, 
c2  C1 »  are  known.  Again,  these  results  are  summarized  along  with 
others  of  interest  in  Table  III. 

The  applications  of  the  conditions  of  criticality  derived  here  are 
reported  in  subsequent  chapters  in  the  interpretation  of  test  results. 

It  is  pertinent  at  this  point  however  to  internret  and  evaluate  the 
significance  of  the  theoretical  results. 
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CRITIQUE  OF  VISCOELASTIC  FRACTURE  THEORY 

Due  to  the  limitations  and  assumptions  required  to  develop  the  solu¬ 
tions  presented,  it  is  worthwhile  to  investigate  the  implications  of  some 
of  these  limitations  to  the  application  of  the  results  to  real  materials 
in  typical  test  situations. 

CRACK  GEOMETRY 

The  spherical  and  cylindrical  flaws  analyzed  here  for  the  purpose 
of  obtaining  a  simple  stress  distribution  are  somewhat  limited  in  analogy 
with  typical  crack  type  failures.  However,  the  ability  to  separate  en¬ 
ergy  and  geometry  for  these  limited  cases  suggests  the  applicability  of 
the  results  to  crack-type  flaw  geometries.  Also,  for  the  elastic  case, 
it  has  been  shown  (2)  that  up  to  a  geometrical  constant,  the  critical 
stress  for  the  crack  geometries  with  crack  length  a0  is  of  the  same  form 
as  the  spherical  or  cylindrical  flaws  of  radius  a. 

Using  quasi-elastic  analytical  techniques,  Rocketdyne  (22),  has 
demonstrated  the  applicability  of  the  fracture  mechanics  technique  to 
predicting  thermal  bore  cracking  in  propellant  grains.  Similar  approaches 
have  been  pursued  by  Swanson  (23)  and  others. 

EFFECTS  OF  FINITE  STRAINS 

Although  the  effects  of  finite  strains  on  the  viscoelastic  solution 
are  presently  beyond  the  state  of  analytical  capability,  some  insight 
into  the  effects  of  large  strains  on  the  criticality  for  a  spherical 
flaw  have  been  obtained  (24)  for  a  body  of  neo-Hookean  material.  The 
essence  of  the  result  of  that  analysis  is  that  for  small  cavities  one 
needs  to  account  for  the  finite  strains.  Williams  has  suggested  an 
extension  of  the  finite  elastic  strain  result  through 

o(\,t)  =  EreJ  (t,)  r  (x)  (85) 

where  X  =  1  +  e  for  small  strains.  These  effects  are  currently  under 
investigation  (24).  , 

Flaw  Growth 

Since  body  forces  were  neglected,  a  degree  of  limitation  was  placed 
on  the  analysis  for  flaw  growth.  Blatz  has  investigated  the  result  for 
both  Hookean  and  a  neo-Hookean  body  with  a  spherical  flaw  including 
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inertial  effects  (25)  and  found  (for  the  Hookean  material)  that  the 
initial  crack  velocity  is  linear  in  time  and  proportional  to  the  square 
of  the  shear  wave  speed  in  the  material. 

An  additional  limitation  of  the  present  analysis  arises  from  the 
flaw  geometry.  The  propagation  of  the  flaw  by  uniform  increase  of  the 
radius  could  be  of  significant  quantitative  difference  from  that  by 
cracking  although  qualitative  similarities  may  exist  (l). 

SURFACE  ENERGY  AND  CUMULATIVE  DAMAGE 

In  the  foregoing  results,  we  were  able  to  express  the  criticality 
criteria  for  a  given  geometry  and  boundary  condition  in  terms  of  an 
equality  between  "constant"  terms  and  time-temperature-loading  history 
dependent  terms.  From  a  purely  theoretical  standpoint,  it  is  not  pos¬ 
sible  to  deduce  a  dependency  on  time  or  temperature  on  the  characteristic 
surface  energy,  In  any  event,  it  would  appear  sufficient  to  establish 
Y  for  the  failure  condition  only  since  its  role  in  the  power  balance  up 
to  failure  is  nil.  Thus,  the  "energy"  term  invoking  the  time  dependent 
material  properties  and  the  loading  history  will  account  for  the  cumula¬ 
tive  effects. 
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section  in 

POKER  CHIP  TESTS ‘and  ANALYSES 

The  first  tests  run  on  the  Cumulative  Damage  program  were  tensile 
tests  in  which  poker  chip  specimens  were  used.  Specimens,  4-inch  diameter 
by  l/3-inch  thick,  are  bonded  between  two  rigid  steel  platens.  When  the 
platens  are  separated,  a  tensile  strain  is  generated  in  the  direction  of 
the  thin  specimen  dimension,  and  an  essentially  hydrostatic  stress  field 
is  induced  in  the  center  portion  of  the  incompressible  material  (26). 

This  is  a  particularly  appealing  test  because  the  hypothesized  stress 
field  of  the  spherical  flaw  derivations,  a  uniform  triaxial  tension, 
matches  that  of  the  specimen. 

The  experim.  jtal  technique,  including  data  collection,  is  described 
in  this  section;  and  data  interpretations  and  quantitative  comparisons 
with  the  analytical  expressions  are  made. 

EXPLitlMENTAL  r»J)GRAM 

The  poker  chip  specimens,  each  with  an  oversized  diameter,  were 
milled  flat  on  a  specially  modified  Var.  Norman  12  Universal  Milling  Ma¬ 
chine  from  bulk  propellant  furnished  by  Lockheed.  Loose  oxidizer  par¬ 
ticles  were  removed  from  the  bond  surfaces  by  impregnating  a  layer  of 
surgical  gauze  with  an  epoxy  resin  in  contact  with  the  propellant  surface. 
After  a  30-mnute  cure,  the  epoxy  resin  was  easily  removed  from  the  speci¬ 
men,  and  the  loosely  imbedded  oxidizer  crystals  adhered  to  the  epoxy  resin. 

It  is  normally  possible  to  prime  the  bonded  surfaces  at  this  point 
and  subsequently  to  assemble  the  specimen  to  the  steel  platens  using  the 
same  epoxy  adhesive.  However,  it  was  necessary  to  expose  this  particular 
propellant  to  pure  oxygen  for  48  hours  before  an  adequate  bond  could  oe 
obtained.  Specimens  for  all  tests  reported  herein  were  exposed  to  oxygen 
in  this  manner  before  they  were  r-’med  with  Eastman's  adhesive.  The 
Eastman  910  primer  layer  was  air  ar^ed  for  12  hours  followed  by  1  hour  at 
170  F.  The  specimen  was  then  bonded  to  the  steel  platen  while  being 
aligned  in  a  special  bonding  fixture  such’ as  that  shown  in  Fig.  4.  A 
light  pressure  was  maintained  on  the  assembly  while  the  epoxy  adhesive 
was  curing.  The  assembly  was  placed  in  instrumentation  rings  (Fig.  5)  after 
the  edges  of  the  propellant  had  been  trimmed.  The  specimen  was  then  sub¬ 
jected  to  tensile  loads  on  a  20,000-pound  Model  TTD  Instron  tester  at  the 
desired  rate  and  temperature. 
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A  typical  data  acquisition  trace  is  shown  in  Fig.  6.  These  data 
are  recorded  as  a  function  of  time.  The  three  LVDT  traces  indicate  the 
amount  of  displacement  in  the  direction  of  the  pull  at  three  points 
equally  spaced  around  the  specimen.  These  traces  must  be  parallel  tc 
ensure  achieving  a  uniform  separation  during  the  test.  The  output  from 
the  Instron  load  cell  is  indicated  by  the  over-all  load  trace  which  repre¬ 
sents  the  total  load  exerted  over  the  entire  cross  section  of  the  speci¬ 
men.  The  fifth  trace  shown  is  the  output  from  a  center  load  cell  located 
at  the  center  of  one  of  the  steel  end  platens.  This  cell  has  an  area  of 
0.1  sq  in.  and  is  bonded  to  the  propellant  specimen.  The  element  is 
flush  with  the  surface  of  the  platen  and  is  very  stiff  so  there  are  essen¬ 
tially  no  modifications  in  the  stress  distribution  due  to  the  presence 
of  the  cell.  The  output  from  the  center  load  cell  reached  a  peak  well 
ahead  of  the  peak  in  the  "over-all  load"  as  indicated  by  the  Instron  load 
cell . 


This  phenomenon  has  been  noted  in  many  constant  stress  rate  tests 
with  a  variety  of  composite  propellants  over  the  past  few  years.  In 
general,  during  constant  stress  rate  tests  it  has  been  noted  that  the 
center  load  cell  stress  is  greater  than  the  average  stress  P/A,  usually 
by  a  factor  just  less  than  2:1.  An  elastic  analysis  relating  this  ratio 
to  the  material  properties  and  specimen  geometry  has  been  given  by  Lindsey, 
et  al  (26),  and  Messner  (27)*  The  ratio  appears  to  hold  essentially  con¬ 
stant  to  the  point  where  microvoids  begin  to  grow,  generally  at  internal 
points  distant  from  the  edges,  resulting  in  changes  in  the  load  distri¬ 
bution.  The  point  of  maximum  output  on  the  center  load  cell  trace  has 
previously  been  interpreted  as  the  poirt  of  failure  (see  Fig.  6,  Point  a). 
However,  during  the  work  with  constant  loads  (Page  53) >  X-rays  of  speci¬ 
mens  revealed  that  the  growth  of  numerous  small  voids  on  the  order  of 
0.040- inch  diameter  led  to  the  decrease  in  the  load  cell  reading.  Based 
on  these  findings,  the  triaxial  failure  point  has  been  defined  at  a  point 
earlier  on  the  curve  where  it  just  begins  to  break  from  a  straight  line, 
Point  b,  Fig.  6. 

There  was  no  visible  indication  of  failure  on  the  outside  of  the 
specimen  at  the  point  previously  discussed.  However,  the  test  was 
continued  until  the  specimen  was  pulled  completely  apart.  A  typical 
surface  after  complete  separation  is  shown  in  Fig.  7;  the  propellant 
layer  is  thin  over  most  of  the  top  platen.  Such  separations  were  en¬ 
countered  quite  frequently;  however,  flaw  growth  failures  near  the 
center  of  the  specimen  in  other  tests  produced  essentially  the  same 
terminal  results.  Therefore,  the  actual  plane  of  separation  is  not 
deemed  pertinent  to  the  test  results. 
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Axial  Deformation,  in. 
a.  Typical  X-Y  Plot  (Moseley  Autograph 


Recorder).  These  Data  Obtained 
Simultaneously  with  That  Collected 
with  the  CEC  Recorder 


Figure  6.  Date  Traces  Collects  During  Poker  Chip  Tests 
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Figure  7.  Typical  Surface  Following  Complete 
Separation  of  Platens 
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DATA  INTERPRETATION 

Experimental  test  results  from  the  elemental  poker  chip  tests  are 
summarized  in  Table  IV.  These  experimental  data  can  be  analyzed  in 
several  ways  insofar  as  the  use  of  the  spherical  flaw  theory  is  con¬ 
cerned.  In  Tables  II  and  III  (a  summary  of  the  criticality  equations) 
it  is  seen  that  either  stress  or  displacement  strain  histories  can  be 
used.  These  equations  have  been  derived  from  both  linear  and  quadratic 
(in  time)  histories.  When  an  attempt  is  being  made  to  evaluate  a  new 
theory,  it  would  be  most  desirable  to  use  the  simplest  case  possible  for 
the  experimental  input,  namely  the  expression  for  linear  displacement. 
Two  factors,  however,  influenced  the  selection  of  the  expression  for  the 
quadratic  stress  history  for  determining  the  failure  parameter,  y/aQ. 


TABLE  IV.  SUMMARY  OF  EXPERIMENTAL  TEST  RESULTS 
FOR  THE  ELEMENTAL  POKER  CHIP  TESTS 


Test 

No. 

Test 
Temp, 
jdeg  F 

Nominal 
Crosshead 
Rate, 
in. /min 

Linear 
Stress 
Rate, 
psi/ sec 

Time  to 
Failure, 
sec 

Center 
Stress  at 
Failure 

Average 
Strain 
at  Failure, 
in. /in. 

35 

77 

0.02 

1.9 

92.0 

175 

0.0080 

17 

77 

0.2 

17.5 

12.5 

218 

0.0084 

18 

77 

0.2 

20.1 

11.5 

223 

0.0082 

29 

77 

0.2 

16.5 

14.5 

200 

0.0069 

31 

77 

0.2 

17.6 

12.5 

220 

'll 

32 

77 

0.2 

19.2 

11.9 

228 

„  j084 

33 

77 

0.2 

19.2 

10.8 

208 

0.0070 

34 

77 

20.0 

2950.0 

0.120 

354 

0.0041 

42 

0 

0.02 

2.3 

127.0 

286 

0.0050 

41 

0 

0.2 

24.5 

15.3 

374 

0.0044 

43 

0 

20.0 

4690.0 

0.145 

650 

0.0023 

38 

140 

0.02 

1.7 

75.0 

124 

0.0090 

36 

140 

0.2 

18.5 

8.2 

152 

0.0081 

40 

140 

0.2 

19.2 

8.1 

156 

0.0088 

37 

140 

20.0 

2235.0 

0.135 

302 

0.0085 

39 

140 

20.0 

2435.0 

0.103 

250 

0.0069 
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First,  it  was  noted  that  the  strain  history  imposed  on  the  poker 
chip  specimens  in  the  Instron  tester  were  not  linear  as  one  would  expect 
to  result  from  a  constant  rate  of  crosshead  displacement.  It  was  deter¬ 
mined  that  this  was  caused  by  the  relatively  low  compliance  of  the  test 
machine  and  linkage  as  compared  to  the  stiffness  of  the  propellant 
specimen.  (The  detailed  investigation  into  this  problem  is  reported  in 
Appendix  IV.)  The  extent  of  this  influence  is  shown  in  Fig.  8  and  9» 
typical  stress  and  strain  time  relationships.  Dotted  lines  represent 
the  linear  rate  corresponding  to  the  failure  points  that  would  have  been 
obtained  if  such  a  history  had  been  possible. 

The  stress  history  was  then  selected  because  the  geometric  parameter, 
k,  can  conveniently  be  omitted  from  the  calculation  without  causing  ex¬ 
cessive  error.  From  Tables  II  and  III,  the  exprejsion  for  criticality 
of  the  spherical  flaw  following  a  stress-prescribed  history  can  be 
written  in  the  form  Critical  Stress  =  (l  -  k3)  f (mechanical  properties 
and  geometry);  but  following  a  displacement-prescribed  history,  critical 
Strain  =  (k*)  g(mechanical  properties  and  geometry). 

Because  k,  the  ratio  of  the  initial  flaw  size,  a0,  to  the  outer 
radius  of  the  sphere,  b,  used  for  the  mathematical  derivations  is  ex¬ 
pected  to  be  rather  small  (<  l),  it  is  readily  apparent  that  it  may  be 
neglected  in  the  equation  for  stress  loading  but  not  for  displacement 
loading. 

Therefore,  the  first  analysis  of  poker  chip  test  results  for  y/80  * 
utilizes  the  equation 


(1  -  2  V) 
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the  expression  for  the  time,  tQ,  to  incipient  flaw  growth  following  a 
stress  history  of  the  form 

CT(t)  =  Oo  (Cjt  +  C2t2)  (87) 

The  meaning  of  the  various  terms  is  as  described  in  Section  II. 

Table  V  shows  corresponding  calculated  values  for  y/ao  for  the  16 
tests.  These  data  reflect  temperature  and  rate  changes  in  the  experimental 
condition  from  0  to  140  F  and  from  0.02  to  20.0  in ./min  nominal  :rosshead 
speeds.  Values  assigned  to  the  creep  compliance  integrals,  (t0), 

were  calculated  using  uniaxial  relaxation  data  supplied  by  Lockheed. 

(These  calculations  are  described  and  results  are  shown  in  Appendix  V.) 
Times  to  failure  were  shifted  using  the  ap  values  also  taken  from  the 
uniaxial  data  mentioned.  From  the  values  calculated  for  y/a0,  it  appears 
possible  that  temperature  and  rate  influence  the  parameter.  For  example, 
at  0  F  the  average  y/a0  was  64.8  in.-lb/m.  ,  while  at  77  F  the  average 
value  had  dropped  to  48.3  in.-lb/cu  in.,  and  at  140  F  the  average  was  down 
to  40.4  in.-lb/cu  in.  The  influence  of  rate  effects  is  not  quite  so  evi¬ 
dent,  but  there  seems  to  be  a  tendency  for  the  higher  rates  to  exhibit  a 
higher  y/a0  a+'  eac6  temperature.  It  should  be  noted  that  utilization  of 
the  uniaxial  stress  relaxation  modulus  data  in  the  triaxial  (or  biaxial 
near  the  flaws)  stress  fields  in  arriving  at  the  values  for  y/a0  may 
actually  be  responsible  for  these  apparent  effects 

LINEAR  VS  NONLINEAR  AN \ LYSIS 

It  was  observed  that  the  loading  rate  for  these  tests  varied  as 
shown  in  Fig.  8  and  9.  Because  these  histories  are  probably  typical 
for  poker  chip  tests,  it  was  desirable  to  determine  the  influence  of  the 
variation  if  a  linear  history  to  failure  had  been  assumed  end  the  data 
interpreted  accordingly.  Therefore,  the  failure  data  were  substituted 
into  the  equation  (from  Table  II ): 


o  ct 
o  o 

1-k3 


(88) 


which  was  then  solved  for  Y/a0.  (Again,  the  various  terms  are  described 
in  Section  II.) 
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Table  VI  shows  the  results  which  vary  only  slightly  from  those 
computed  in  Table  V.  The  conclusion  must  be  that  the  slightly  different 
loading  histories  involved  in  these  tests  do  not  essentially  influence 
failure  stresses  and  strains  or  the  parameter  Y/a0- 

The  critical  valuer,  for  Y/a0  represent  the  failure  criteria  for  the 
corresponding  rates  and  temperatures.  Therefore,  a  composite  of  these 
data  would  represent  the  failure  envelope  over  the  range  of  the  tests. 
Figures  10  and  11  show  these  failure  envelopes  for  different  Y/a0  values. 
Figure  10  shows  failure  predictions  using  an  average  Y/a0  of  51.1  (the 
average  value  for  the  77  F  data-linear  analysis).  The  use  of  this  value 
assumes  that  there  is  no  temperature  effect  on  Y/a0.  Therefore,  the 
experimental  points  for  0  and  140  F  do  not  fall  on  the  failure  envelope. 
The  curve  indicates  the  prediction  error  expected  as  a  result  of  tempera¬ 
ture  effects  on  Y/a0.  Figure  11  presents  the  failure  envelopes  at  dif¬ 
ferent  temperatures  using  the  corresponding  Y/a0  for  the  temperatures 
shown.  As  expected,  the  experimental  data  points  fit  the  envelope 
better  since  the  values  for  the  corresponding  Y/a0  were  determined  using 
averages  from  the  data  at  each  temperature. 

Having  obtained  Y/aQ,  it  is  now  possible  to  calculate  the  value 
for  k.  The  equation  for  criticality  following  a  second  order  displace¬ 
ment  history  provides  the  necessary  tool.  Using  the  average  Y/a0  for  the 
77  F  data  (51. 1  in. -lb/in. 5)  and  the  second  order  strain  history  equa¬ 
tion,  values  for  k  were  calculated.  These  values  are  shown  in  Table  VII. 
If  the  effective  outer  sphere  radius,  b,  is  assumed  to  be  one-half  the 
poker  chip  thickness,  0.167  inch,  an  estimate  for  the  characteristic 
flaw  size  can  be  deduced  from  the  calculations  of  Table  VII. 


a 

Y  =  0.3  b  =  0.167 


a  =  0.055  inch 
o 


It  seems  logical  that  the  assumed  value  for  b  is  an  upper  limit  (for 
the  poker  chip  specimen),  implying  that  the  calculated  value  for  the 
flaw  size  is  also  an  upper  limit.  This  flaw  size  seems  to  be  reasonable 
because  the  larger  solid  particles  within  the  propellant  are  of  the 
order  of  600  microns  (0.024  inch)  in  diameter.  The  small  variation  be¬ 
tween  the  several  calculated  k  values  also  adds  to  the  credibility  of 
the  calculations. 
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TABLE  VII 


On->  further  interesting  calculation  can  be  made  using  the  value  for 
a0  found  above.  Again  using  Y/a0  -  51*1  in. -lb/in. 3  and  a0  -  0.055  inch, 
one  finds  the  characteristic  strain  energy  release  rate  to  be 

Y  =  (51.1  in. -lb/in.3)  (0.055  inch) 

=  2.81  in. -lb/in. 3 

which  one  should  also  interpret  as  an  upper  limiting  value  based  on  the 
rather  gross  assumptions  that  have  been  described. 

CONSTANT  LOAD  CONFIRMATORY  TESTS 

Three  constant  load  poker  chip  tests  were  conducted  to  confirm 
failure  prediction  based  on  critical  values  for  Y/a0.  These  tests  repre¬ 
sent  a  more  complex  loading  history  than  the  simple  pulls  used  in  deter¬ 
mining  the  failure  parameter.  Loads  of  900,  750,  and  600  pounds  were 
oonlied,  relatively  instantaneously  (loading  time  <  15  seconds),  across 
the  specimens.  The  loads  exerted  an  average  stress  of  P/A  as  indicated 
in  Fig.  12  for  each  test  Data  collected  during  these  tests  included 
the  output  of  the  center  load  cell  and  displacement  across  the  specimen 
in  the  direction  of  the  pull.  Center  load  cell  output  was  used  to  cal¬ 
culate  the  stress  at  the  center  of  each  specimen. 

The  rapid  drop  early  in  the  900-  and  750-pound  tests  indicated  that 
failure  had  occurred.  The  specimen  was  X-rayed  perpendicular  to  the  axis 
of  symmetry  during  the  750-pound  test  to  verify  this  interpretation. 

Data  revealed  by  X-ray  are  presented  schematically  as  a  function  of  time 
in  Fig.  13.  It  is  generally  considered  that  flaws  0.01  of  the  propellant 
thickness  (in  this  case  0.04  inch)  are  the  smallest  that  can  be  detected 
by  X-ray.  By  extrapolating  back  to  the  smaller  (undetectable)  flaws, 
one  can  surmise  that  the  initial  flaw  size  must  have  been  on  the  order 
of  0.025  inch  which  is  consistent  with  previous  findings.  The  flaws 
continued  to  grow  in  both  size  and  number  for  the  duration  of  the  test. 

Figure  14  shows  test  results  compared  to  those  predicted.  The  pre¬ 
dicted  values  of  0.25  and  2.4  minutes  agree  quite  well  for  the  900-  and 
750-pound  tests,  respectively.  The  failure  point  was  taken  to  be  the 
time  when  the  slope  of  the  center  load  <vll  output  first  became  negative. 
The  prediction  for  the  time  to  failure  in  the  600-pound  test  did  not 
materialize — the  specimen  had  not  failed  after  3000  minutes,  and  the 
test  was  discontinued. 
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One  other  significant  aspect  of  the  curves  in  Fig.  12  involves  the 
plots  resulting  from  the  600-pound  test.  The  output  of  the  center  load 
cell  continued  to  rise  throughout  this  test.  The  more  incompressible 
the  specimen  material  is  (i.e.,  K»  G),  the  larger  the  center  load  cell 
output  will  be  (28).  The  limiting  value  will  be  2  x  P/A  when  K  *-*  co  or 
G  -*  0.  The  continuing  increase  in  the  center  load  cell  reading  in  the 
600-pound  test  was  interpreted  to  indicate  that  no  flaw-type  failures 
were  present  and  that  the  shear  modulus,  G,  was  relaxing  at  a  decidedly 
faster  rate  than  the  bulk  modulus,  K. 

Calculations  of  predicted  times  to  failure  during  these  tests  were 
based  on  the  criticality  formula  for  the  time  to  failure,  t0,  following 
a  constant  stress  test 


(t)  -  D 


2  D 


crp 


g 


(90) 


although  the  stress  history  as  recorded  by  the  center  load  cell  was 
certainly  not  constant.  However,  this  calculation  should  provide  a 
reasonable  estimate  of  the  times  to  failure. 
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SECTION  IV 

STATIC  LOADINGS  ON  BIAXIAL  STRIP  SPECIMENS 

In  the  previous  discussion  and  interpretation  of  the  poker  chip 
results  it  was  shown  that  small  deviations  from  a  linear  loading  history, 
of  either  the  stress  or  strain,  e<ert  little  influence  on  the  failure 
stresses  and  strains  predicted  on  the  basis  of  the  straightline  history. 
This  provided  the  license  necessary  to  draw  conclusions  about  the  values 
of  theoretical  parameters  with  only  the  point  values  of  failure  stresses 
and  strains  following  constant  strain  rate  tests  of  biaxial  strips. 

Several  of  these  tests  were  conducted  and  reported  by  Lockheed.  Inter¬ 
pretation  of  the  resulting  data  to  obtain  the  fracture  parameters  appro¬ 
priate  for  the  biaxial  specimen  will  be  discussed  first.  Then  by  using 
the  results  of  the  elemental  tests  the  analytical  expressions  of  Tables 
II  and  III  make  it  possible  to  predict  the  time  to  failure  following  the 
bilinear  strain  histories.  Two  tests  of  this  nature  were  performed,  and 
the  measured  and  predicted  failures  are  indicated. 

CONSTANT  RATE  TESTS 

The  failure  data  on  which  the  calculations  of  these  sections  are 
based  were  collected  by  Lockheed  during  the  characterization  phase  of 
the  Cumulative  Damage  program.  The  small  biaxial  strip  specimens,  which 
measured  7-5  x  0.25  x  1.5  inches  long,  were  tested  in  tension  in  an 
Instron  machine.  Strain  rates  were  varied  over  a  4-decade  range  by  set¬ 
ting  the  crosshead  speeds  to  0.01,  0.1,  1.0  and  10.0  in ./min.  Their 
magnitudes  were  calculated  by  dividing  these  speeds  by  the  specimen  length, 
1.5  inches.  The  tests  were  run  at  nine  temperatures  varying  from  -100 
to  200  P.  Some  of  those  tests  run  at  1  in ./min  (three  replicate  speci¬ 
mens  at  each  temperature)  were  stopped  when  the  strain  level  reached 
2. 5}£, and  the  resulting  force  relaxation  was  monitored  with  time.  After 
a  pause  of  about  10  minutes,  the  constant  rate  loading  was  resumed  and 
continued  until  failure.  No  variation  could  be  detected  between  the 
results  cf  the  continuous  and  the  interrupted  tests. 

In  general,  the  maximum  stress  was  used  as  the  criterion  for  failure. 
Observers  at  Lockheed  reported  only  slight  differences  from  the  time  of 
the  observation  of  the  first  crack  and  the  maximum  stress  level.  The 
reported  failure  stresses  and  failure  strains  were  plotted  against  tem¬ 
perature  in  Fig.  15  which  was  taken  directly  from  (29). 
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DATA  INTERPRETATION 

Knowing  both  the  failure  stream ?s  and  strains,  it  is  possible,  pro¬ 
viding  the  loading  history  is  known,  *o  use  the  fracture  equations  of 
Tables  II  and  III  to  solve  for  both  y/a0  f.ad  k.  The  exact  time  history 
of  neither  the  stress  nor  the  strain  was  available.  However,  because  the 
strains  were  imposed  by  a  constant  crosshead  movement  rate,  because  the 
strip  specimen  is  small,  and  because  the  propellant  is  relatively  compli¬ 
ant  it  seems  reasonable  to  assume  a  constant  strain  rate  history.  The 
force  response  to  a  linear  strain  history  on  a  biaxial  strip  at  room 
temperature,  shown  in  Fig.  6  ,  Appendix  IX,  reveals  that  under  this  con¬ 
dition  the  viscoelastic  effects  are  not  too  great  and  the  force  is  rela¬ 
tively  linear  in  time.  Linear  rates  were  assumed  for  both  stress  and 
strain  because  slight  variations  from  a  linear  history  are  believed  to 
make  only  slight  differences  in  the  results  (see  Section  III).  Hence, 
the  formulas  for  criticality  of  a  flaw  in  a  plane  stress  body  (Tables  II 
and  III)-  following  linear  histories  were  applied. 

Since  both  failure  stresses  and  strains  are  reported,  two  equations 
can  be  written  for  each  test,  i.e.,  equations  having  the  two  unknowns, 
k  and  \/aQ.  By  assuming  that  the  time -temperature  shift  was  appropriate 
for  the  material  properties,  data  at  all  temperatures  were  used.  The 
values  used  in  the  calculations  are  summarized  in  Table  VIII  for  the  strain 
boundary  and  Table  IX  for  the  stress  boundary.  The  two  simultaneous 
equations  originate  with  the  criticality  expressions  of  Tables  II  and  III. 
The  expression  for  the  incipient  growth  of  a  cylindrical  flaw  in  a  plane 
strain  body  following  a  linear  strain  history  is  given  by 
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BLE  VIII.  CALCULATIONS  USING  BIAXIAL  STRIP 
FAILURE  STRAINS 
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Solving  Eq.  91  and  92  for  the  unknown  k,  one  finds 


which  permits  solving  either  Eq.  91  or  92  for  y/a0. 

These  numerical  calculations  are  outlined  in  Table  X  and  the 
properties  of  relaxation  modulus  and  creep  compliance  used  herein  were 
determined  by  biaxial  tests  and  reported  in  Appendix  V. 

The  results  of  the  calculations  tre  summarized  in  Fig.  16  and  17 
with  k  and  y/a0  plotted,  respectively,  showing  a  definite  influence  of 
both  temperature  and  strain  rale.  The  geometrical  parameter,  k  (=  a0/b), 
is  seen  to  be  fairly  constant  m  all  tests.  However,  y/a0  fluctuates 
widely  indicating  a  definite  sensitivity  to  both  rates  and  temperatures. 
Because  the  geometrical  parameter  does  not  vary  widely  it  was  suspected 
that  the  initial  flaw  size  a0  was  also  relatively  constant;  so  most  of 
the  variation  was  credited  to  the  characteristic  strain  energy  release 
rate,  y.  Because  of  the  prevalence  of  opinion  that  this  term  should  be 
a  constant,  a  brief  examination  of  the  influences  that  might  make  it 
appear  to  be  time  and  temperature  dependent  was  made.  These  considera¬ 
tions  included  the  difference  in  the  strens  state  in  a  biaxial  stri  » 
from  the  uniform  stress  assumed  for  the  model,  the  effects  o'  iir.ite 
strains,  and  the  error-producing  effects  that  the  assumption  of  incom¬ 
pressibility  would  give  to  the  apparent  y.  Each  of  thes ;  is  discussed. 
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Figure  16.  The  Values  Calculated  for  k  Using  Prcpellair 
Failure  Data  Furnished  by  Lockheed.  (The 
biaxial  strip  tests  were  run  with  the  cross- 
head  set  to  move  at  a  constant  rate.) 
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Figure  17-  The  Values  Calculated  for  >/a©  l*  leg 
Propellant  Failure  Data  Famissed  Sr 
Lockheed.  (The  biaxial  strip  tests 
vere  run  vith  the  crosshead  set  zc 
sore  at  a  cods  tact  r-ite.) 
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UNIFORM  TENSION  ONLY 

The  basic  model  from  which  the  mathematics  were  derived  was  con¬ 
sidered  in  an  effort  to  determine  a  plausible  explanation  for  the  varia¬ 
bility.  The  axially  symmetric  geometry  (see  Fig.  18a)  is  subjected  to 
a  uniform  tensile  stress.  The  fact  that  the  resulting  stress  field  is 
axially  symmetric  provides  the  basis  for  a  mathematically  tractable  solu¬ 
tion.  It  should  be  noted,  however,  that  the  actual  state  of  stress  in  a 
biaxial  strip  specimen  is  quite  different  from  the  hypothesized  model. 

If  the  strip  is  considered  to  be  very  wide  so  the  strain  in  the  x-direc- 
tion  (i.e.,  e.x  distant  from  the  edges)  is  taken  as  zero,  the  state  of 
stress  resulting  from  a  tensile  displacement,  v,  in  the  y-direction  can 
be  deduced.  To  preserve  constancy,  Poisson's  ratio  (v)  is  set  equal  to 
0.5*  Then  using  Hooke's  law  for  the  plane  stress  (thin)  specimen  (0  =  0, 

where  x  is  the  coordinate  taken  through  the  thickness  of  the  specimen) 
one  can  deduce  using 


that 


CT  a  V  0 

*  y 

or  with  v  =0.5 


(94) 


(95) 

(96) 


This  is  the  3tate  of  stress  illustrated  in  Fig.  18b  and  is  quite  different 
from  the  uniform  stress  of  the  hypothesized  model.  Just  how  different  can 
best  be  visualized  by  separating  the  stresses  of  Fig.  18b  into  the  uniform 
portion  and  the  shear  portion  (see  Fig.  19).  This  suggests  that  the 
results  of  the  criticality  calculations  might  be  improved  if  only  the  uni¬ 
form  tension  component  of  the  stress  state  in  the  biaxial  strip  were  used. 

Using  this  approach,  the  calculations  were  simplified  by  assuming 
k  =  0.6.  This  is  a  reasonable  assumption  in  view  of  the  values  plotted 
in  Fig.  16.  The  resulting  calculations  are  summarized  in  Table  XI.  While 
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TABLE  XI.  CALCULATION  OF  y/&0  USING  ONLY 
THE  UNIFORM  STRESS  FROM  THE 
BIAXIAL  STRIPS 


Temperature, 
deg  F 

Strain 

Hate, 

in ./ in ./min 

Failure 
Stress , 

ViV 

psi 

Uniform 
Stress , 

■?  a  c.t  , 

4  0  1  o’ 

psi 

2  (7-)  (t  )  -  (t  ), 

xtQ'  crp  '  o'  ^  2  crp  v  o' 

0 

psi 

y/«v 

in.-lb/cu  in. 

-100 

6.67 

1,455 

1,090 

0  10 

63.8 

0.667 

932 

700 

o.o„„oi66 

39.7 

0.0667 

578 

434 

0.0000332 

30.5 

0.00667 

850 

637 

0.0000500 

99.1 

-65 

6.67 

836 

627 

0.0000355 

68.1 

0.0667 

679 

509 

0.000103 

130 

0.00667 

514 

385 

0.000172 

124 

-40 

0.667 

516 

388 

0.000186 

136 

O.O667 

418 

314 

0.000312 

150 

0 

6.67 

505 

378 

0.000312 

218 

O.667 

368 

276 

0.000455 

169 

0.0667 

278 

208 

0.000810 

171 

0.00667 

195 

146 

0.00105 

109 

40 

0.667 

201 

151 

0.00C710 

79 

0.0667 

158 

118 

0.00117 

79.5 

0.00667 

121 

90.8 

0.00166 

66.8 

77 

6.67 

183 

137 

0.000750 

68.7 

O.667 

123 

92 

0.00102 

42.2 

0.0667 

109 

81.7 

0.00152 

49.5 

0.00667 

89.3 

67 

0.00218 

47.8 

100 

6.67 

163 

122 

0.000820 

59.6 

0.667 

120 

90 

0.00132 

52.2 

O.O667 

88.6 

66.5 

0.00167 

36.1 

O.OO667 

77 

57.7 

0.00235 

38.2 

150 

6.67 

123 

92 

0.00142 

58.7 

0.667 

115 

86.2 

0.00204 

74 

0.0667 

78.3 

58.6 

0,00255 

42.7 

0.00667 

63.5 

47.6 

0.0033 

36.5 

200 

6.67 

121 

90.5 

0.00253 

101 

0.667 

102 

76.5 

0.0033 

94.3 

0.0667 

72.7 

54.6 

0.00386 

56.2 

0.00667 

62.7 

47 

0.0045 

48.5 
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it  xs  evident  that  this  approach  yields  v/a0  values  more  constant  than 
those  indicated  in  Table  X,  it  is  readily  obvious  that  the  apparent 
characteristic  strain  energy  release  rates  are  not  nearly  constant  enough 
for  practical  use. 


Idealized  Model  on  which  the 
Mathematics  are  Based  (Uniform 
Tensile  Stress) 


^ — ►  x,u 


Specimen 


Platens 


b.  Model  Appropriate  for 
Biaxial  StriD  (2:1 
Stress  Ratio) 


Figure  18.  States  of  Stress  for  Idealized  and  Realistic 
Models 
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LARGE  STRAINS 

Another  factor  certain  to  influence  the  results  is  the  inapplicability 
of  the  linear  viscoelasticity  at  the  higher  strain  levels  inherent  in  the 
biaxial  strip  tests.  The  importance  of  accounting  for  the  finite  strains, 
even  in  the  poker  chip  specimen,  is  discussed  in  detail  by  Blatz  in  Appen¬ 
dix  I.  The  great  uignificance  of  the  finite  strains  in  the  biaxial  strip 
where  the  gross  strain  at  failure  is  on  the  order  of  JOfi  is  readily  ap¬ 
parent,  especially  when  it  is  contrasted  with  the  poker  chip  where  the 
gross  strain  at  failure  is  only  about  2$. 

One  approach  to  determining  the  significance  of  this  effect  is  to 
contrast  the  example  y'a  calculated  using  Hookean  and  neo-Hookean  formu¬ 
lations  in  Fig.  5  and  4  in  Appendix  I.  When  the  modulus  is  low  (higher 
strains)  the  solution  which  takes  finite  strains  into  account  gives  a 
significantly  different,  and  higher,  value  for  the  characteristic  strain 
energy  release  rate.  Consequently,  it  is  believed  that  large  deformations 
are  important  to  fracture  predictions  in  propellants,  especially  at  free 
grain  boundaries  such  as  the  inner  port  surface.  This  influence  might 
be  a  primary  one  in  creating  the  apparent  rate  and  temperature  sensitivity 
of  y. 

COMPRESSIBILITY 

Blatz  also  discussed  one  other  factor  which  may  cause  y  to  appear 
time  and  temperature  dependent  in  these  tests.  That  is  the  assumption 
of  incompressibility  that  was  made  to  formulate  the  criticality  equations 
of  Tables  II  and  III. 

It  is  universally  agreed  that  extensive  dewetting,  hence  dilatation, 
is  present  in  the  biaxial  strip  under  large  tensile  strains.  Blatz 
suggests  that  to  account  for  the  volume  change  in  an  elastic  material 
one  can  simply  replace  the  modulus,  E,  by  E/[2(l  -  v)]  .  If  the  dewetting 
process  is  a  function  of  time  and  temperature  it  is  clear  that  a  calcu¬ 
lation  of  y  which  disregards  this  effect  will  yield  values  which  appear 
to  be  functions  of  these  variables  even  though  it  night  be  completely 
independent  of  them. 


SUMMARY 

The  implications  of  each  of  the  three  influences  are  obviously  impor¬ 
tant.  But  it  is  clear  that  definitive  statements  concerning  what  appears 
to  be  a  time  and  temperature  dependence  by  the  characteristic  strain 
energy  release  rate  as  deduced  from  biaxial  strip  tensile  tests  cannot 
be  made  until  a  careful  study  of  these  effects  has  been  completed. 


72 


•♦**-*•* •owai 


-  >s5j»jw*4' tt»  rji 


BILINEAR  TESTS 

Two  bilinear  strain  rate  tests  were  conducted  using  a  biaxial  strip 
specimen.  This  specimen,  13*25  x  3/8  x  2.8  inches  long,  is  somewhat 
larger  than  that  employed  at  Lockheed;  but  the  proportions  are  not  greatly 
different.  The  cumulative  damage  testing  machine,  which  can  be  prepro¬ 
grammed  to  generate  two  ramp  rates  in  sequence  either  for  the  force  or 
for  the  platen  displacement,  was  used  for  the  tests. 

The  two  rates  were  chosen  about  4  decades  apart  iu  both  tests.  One 
was  run  slow-to-fast,  and  the  other  was  run  fast-to-slow.  Even  greater 
rate  changes  could  have  been  studied,  but  these  variations  appeared 
sufficiently  large  for  study  purposes,  yet  provided  a  high  rate  still  slow 
enough  to  permit  the  visual  observation  of  failure. 

Strain  and  maximum  stress  histories  of  the  two  tests  are  given  in 
Fig.  20.  Both  tests  were  performed  with  the  specimen  at  140  F. 

FAILURE  PREDICTION  (BILINEAR  HATES) 

Time-to-failure  predictions  were  made  using  the  appropriate  expression 
for  criticality  and  the  y/a0  and  k  values  found  in  the  constant  rate  tests. 
The  expression  for  a  plane  stress  body  subjected  to  a  bilinear  strain 
rate  is 


gX-  (X  +  3  k2)  -  <-§)  [.,  (c,  -  C2)  e£>  (t,)  (9?) 

+  C1  c2  Erel  bt)  +  °2  (c2  ~  cl^  ELl  "  V 

where  t^  represents  the  time  at  which  the  strain  buildup  changes  rate. 

The  one  problem  is  which  value  should  be  appropriately  used  for  y/a0 
in  this  calculation.  By  examining  the  values  deduced  from  the  constant 
rate  tests  and  plotted  in  Fig.  17  one  can  sketch  a  curve  of  y/a0  vs 
loading  rate  as  shown  in  Fig.  21.  This  curve  can  then  be  used  to  deduce 
values  corresponding  to  the  terminal  loading  rates  in  the  two  bilinear 
tests .  They  were  found  to  be 
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Test  No.  114: 

C2Uo 

■  2,*76  in. /in. /mm,  =  160 

C 

X. 

in. -lb 

b 

in.3 

Test  No.  116: 

C2Uo 

b 

=  O.OOt'S^S  ir.  /m./im?.  = 

'  a 

in. -lb 
44  3 

When  substituted  into  Eq.  97  the  value*  i  '  iV'ilure  times  of 

Test  No.  114:  tf  -  tj  «  0.1121  *._M 
Test  No.  116:  t^  -  t^  =  460  miu 

these  predictions  are  compared  with  the  *hee~*'«d  failure  strains 

and  times  in  Fig.  22. 

Correlation  of  the  limited  observations  t*.  the  predicted  failures 
;e  net  clearly  indicated.  However,  the  experiences  of  the  associated 
contractors  seemed  to  indicate  a  definite  relationship  between  the 
failure  stresses  (or  characteristic  strain  energy  release  rate)  and  the 
time  a  specimen  surface  has  been  exposed,  particularly  to  humidity. 

^cosure  variations  of  a  similar  order  of  magnitude  as  that  of 
toe  prediction  errors  are  expected  merely  from  the  data  scatter  and  the 
aging  effect. 
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Figure  20.  Stress  and  Strain  Histories  of  the 
Two  Bilinear  Strain  Rate  Tests 
(biaxial  strip  specimen) 
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Test  No.  116 


Figure  20.  Concluded 
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Figure  21.  Variation  of  y/a  vs  Strain  Kate 
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Figure  22.  Concluded 
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SECTION  V 
VIBRATION  FATIGUE 

The  principal  portion  of  Rocketdyne's  effort  on  the  Cumulative  Damage 
program  was  given  to  vibration  fatigue  testing  of  biaxial  strip  specimens 
of  propellant.  These  tests,  which  were  conducted  on  a  Cumulative  Damage 
Tester  fabricated  by  the  MTS  Corporation  especially  for  this  program, 
are  described  below.  Data  from  tests  having  a  displacement-controlled 
input  (rather  than  force)  are  also  reported.  (The  Cumulative  Damage 
Tester  and  its  capabilities  are  discussed  in  Appendix  VI,  and  Data  from 
force-controlled  tests  are  summarized  in  Appendix  VII.) 

BIAXIAL  FATIGUE  TESTING  PROGRAM 


i 


i 


I 


Three  replicate  biaxial  fatigue  tests  were  conducted  at  each  of 
three  vibration  frequencies  and  three  amplitudes  to  find  the  parameters 
y/a0  and  k  necessary  to  use  the  viscoelastic  fracture  expressions  for 
crack  initiation.  The  plan,  Table  XII,  required  replicate  testing  to 
give  a  measure  of  the  reproducibility  of  such  tests. 

TABLE  XII .  PLANNED  VIBRATION  FATIGUE 
TESTS  OF  BIAXIAL  PROPELLANT 
STRIPS  PLATEN  DISPLACEMENTS 
OF  THE  FORM 


Platen  Displacements  of  the  Form: 
u(t)  =  uQ  +  uQ  sinut;  t  >  0 
=  0  ;  t  <  0 

Frequencies , 

cps 

Displacement,  in. 

5 

15 

50 

0.084 

3 

3 

3 

0.098 

3 

3 

3 

0.112 

3- 

_J _ 

_ 

SPECIMEN  PREPARATION 

The  biaxial  strips,  Fig.  23,  were  14  x  2.8  x  3/8  inch,  providing  a 
length- to-width  ratio  of  5»  However,  because  of  the  flared  bond  surface 
(the  cross  section  was  shaped  exactly  like  a  tab-ended  JANAF  specimen), 
it  was  necessary  to  machine  both  sides  of  the  specimen. 
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Figure  23.  Original  Biaxial  Strip  Specimen  Design 
and  the  Revised  Design  Actually  Used 


Preliminary  tests  revaled  that  cracks  would  more  often  than  not 
originate  from  a  site  within  the  body  of  the  propellant  when  the  end  edges 
of  the  specimen  were  reinforced  with  a  soft  rubber,  3M'a  EC1949,  and  con¬ 
toured  as  shown  in  Fig.  23b.  Preliminary  tests  also  showed  that  when  the 
rubber  was  used  to  form  fillets  at  the  bond  surfaces  with  the  metal  platens 
the  expensive  and  tedious  machining  operations  could  be  eliminated.  Con¬ 
sequently,  specimens  used  on  this  program  were  sawed  directly  from  the 
bulk  propellant  blocks  furnished  by  Lockheed.  Only  a  slight  and  inexpen¬ 
sive  milling  operation  was  then  required  to  form  flat,  parallel  sides. 

Each  specimen  was  bonded  to  the  platens  in  a  specially  constructed 
fixture  which  provided  for  proper  alignment  and  dimensions.  While  the 
specimen  was  in  the  bonding  fixture,  a  small  thermocouple  was  placed  on 
the  flat  propellant  surface  and  covered  with  a  thin  layer  of  the  EC1949 
which  was  feathered  as  thin  as  possible  near  its  edges.  The  thermocouple 
was  always  placed  a  slight  distance  from  the  specimen  centerlines,  both 
horizontally  and  vertically,  in  an  attempt  to  deter  crack  initiation  near 
the  thermocouple . 


Figure  24  shows  a  specimen  mounted  between  the  two  platens  ready  for 
testing  in  the  cumulative  damage  tester. 
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Figure  24.  Biaxial  Strip  Specimen  Ready  for  Vibra¬ 
tion  Fatigue  Test  (Edges  are  reinforced 
by  thin  layers  of  soft  EC1949.  The  small 
spot  of  rubber  covers  a  surface-mounted 
thermocouple  from  which  the  temperature 
was  controlled.) 


TESTING 

The  specimen,  with  the  platens,  the  thermocouple,  and  the  edge¬ 
reinforcing  EC1949  in  place,  was  mounted  in  the  loading  frame  of  the 
tester.  Ambient  temperature  within  the  cabinet  was  adjusted  from  the 
output  of  the  surface-mounted  thermocouple  to  bring  the  specimen  tempera¬ 
ture  to  that  required  for  the  test.  During  all  vibration  fatigue  tests, 
this  specimen  surface  temperature  was  held  constant  at  140  F.  Coarse 
observations  of  thermocouple  output  indicated  that  the  temperature  varied 
as  much  as  +  3  F  between  the  cycling  of  the  environmental  controls. 

Figure  25  shows  the  Cumulative  Damage  Tester  including  the  loading 
frame  and  environment  cabinet.  Note  the  observation  window  and  the  two 
ports  for  hand  access. 


Figure  25.  Cumulative  Damage  Tester 

LOADINGS 

Since  these  biaxial  strip  specimen  tests  were  conducted  to  evaluate 
the  usefulness  of  the  criticality  expressions  based  on  energy  balance  for 
predicting  the  service  life  of  propellant  subjected  to  fatigue-type  load¬ 
ings,  it  was  desirable  to  apply  a  strain  history  that  reproduced  the 
strain  history  assumed  for  the  theoretical  derivations  as  closely  as  pos¬ 
sible.  One  of  the  general  displacement  histories  for  which  a  tractable 
solution  for  the  growth  of  a  spherical  or  cylindrical  flaw  is  available  is 

u(b,t)  =  Up  l(t)  +  uq  sin  urt  (98) 

where 

l(t)  =0  t  <  0 
=  1  t  >  0 

which  is  as  shown  in  Fig.  26. 
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Figure  26.  Vibratory  Tensile  Displacement 
Input  for  which  a  Viscoelastic 
Fracture  Formula  is  Available 

Preprogramming  the  exact  displacement  magnitudes,  u  and  u0,  was  more 
difficult  than  had  been  expected;  but  precise  values  could  be  monitored. 
Therefore,  approximate  valutj  were  preprogrammed,  and  calculations  and 
data  interpretations  were  based  on  the  displacements  that  were  actually 
measured. 

Acceptable  data  were  recorded  from  tests  conducted  at  5  and  15  cps. 
But  at  higher  frequencies,  the  vibration  of  the  loading  frame  and  its 
foundation  caused  erroneous  data  to  be  recorded.  For  example,  when  the 
lower  platen  was  vibrated  at  50  cps  at  a  displacement  corresponding  to 
8^  strain  with  no  specimen  in  place,  the  force  cell  above  the  upper  platen 
gave  an  oscillating  output  corresponding  to  forces  of  up  to  40  pounds. 
Therefore,  testing  was  suspended  until  this  problem  was  solved. 

FOUNDATION  PROBLEM 

To  provide  a  seismic  foundation  for  the  loading  frame,  a  11,000- 
pound  mass  of  concrete  was  poured  into  a  cork-lined  cubicle  pit.  This 
structure,  which  was  believed  to  be  highly  damped,  has  a  calculated 
natural  frequency  approximately  midway  between  15  and  50  cps,  both  fre¬ 
quencies  at  which  tests  were  to  be  run.  (These  calculations  were  sub¬ 
sequently  confirmed  by  tests.)  As  discussed  above,  the  frame  and  founda¬ 
tion  performed  well  for  tests  conducted  at  5  and  15  cps;  but  at  higher 
frequencies  they  vibrated  at  unacceptable  levels. 
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Three  approaches  for  alleviating  this  problem  were  considered.  The 
first  was  to  increase  the  mass  of  the  foundation,  thus  lowering  the  natural 
frequency  of  the  foundation-frame  system.  The  second  possible  approach 
was  to  design  a  tuned  energy  dissipater  that  would  absorb  the  input  vibra¬ 
tion  energy  during  testr  run  at  the  precise  frequency  to  which  it  was 
tuned.  Theoretically,  this  would  leave  the  foundation  and  attached  load 
frame  relatively  vibration  free.  Both  of  these  approaches  were  discarded — 
the  first  because  of  the  expense  and  time  involved,  and  the  second  be¬ 
cause  of  the  uncertainty  associated  with  its  predicted  performance. 

The  solution  adopted  was  to  replace  the  cork  lining  with  pneumatic 
mountings  to  provide  the  vertical  support  for  the  concrete  foundation. 
Calculations  indicated  this  should  reduce  the  effective  spring  constant 
as  required  to  lower  the  natural  frequency  of  the  spring-mass  system 
(foundation  and  frame)  from  about  35  to  about  2  cps.  These  calculations 
proved  to  be  realistic;  the  amplitude  of  the  sympathetic  responses  to  a 
50  cps  input  of  the  moving  platen  was  reduced  by  a  factor  of  4  or  5* 

Figure  27  shows  a  photograph  of  one  of  the  pneumatic  springs. 


8-2-68-992 

Figure  27*  One  of  the  Four  Pheumatic 
Springs  Used  to  Soften  the 
Support  Structure  for  the 
Loading  Frame  of  the  Cumu¬ 
lative  Damage  Tester 


86 


^  i.Miif^'v.t^r-^-i.Wir^  .aca&tttt  iA^wwtiw'r 


t  ^ii*^  VS^*^^****  ^  ,*s*r«sv*itza,Vj  >  .,-w  *f  J5r^s?ffi*jva*i(5a.v  > 


TEMPERATURE 


| 

» 

t 


There  was  concern  about  the  cyclic  energy  input  and  its  effect  on 
thermal  gradients  across  the  thickness  of  the  specimen  at  the  higher  fre¬ 
quencies  and  strains.  Since  the  specimen  was  pre-conditioned  to  a  stable 
temperature  of  140  F  before  the  test  was  started,  cyclic  inputs  would 
obviously  increase  the  specimen  temperature.  To  partially  compensate  for 
this  effect,  a  thermocouple  was  attached  to  the  surface  of  the  specimen; 
and  the  temperature  of  the  box  was  controlled  automatically  so  the  speci¬ 
men  surface  temperature  remained  at  140  F.  This  did  not  ensure  that  the 
entire  specimen  temperature  would  be  held  to  this  temperature,  however. 
Since  heat  was  flowing  from  the  specimen  to  the  air,  the  temperature  of 
the  internal  parts  of  the  specimen  could  become  significantly  higher  than 
the  140  F  surface  condition.  Therefore,  the  following  calculation  was 
performed  to  determine  the  maximum  temperature  gradient  across  the  speci¬ 
men  during  one  of  the  higher  frequency  and  strain  tests. 


Figure  28  shows  a  typical  Lissajous  pattern  for  the  force-displace¬ 
ment  relationship  during  a  47  cps  and  6 %  strain  test.  By  measuring  the 
area  of  such  patterns  and  converting  to  energy  density,  in.-lb/in.^,  it 
is  possible  to  establish  a  point  of  equilibrium  for  the  energy  input  as 
shown  in  Fig.  29.  The  equilibrium  energy  input,  q,  can  then  be  used  to 
determine  the  maximum  temperature  gradient  in  the  expression 


T  =  T  + 


2k  v  0 


(Xn2  -  X,2) 


Temperature .  T 


(99) 


where 

T  =  Temperature 
T  =  Surface  temperature,  140  F 

q  =  Equilibrium  energy  input, 
Btu/in.3  hr 

k  =  Conductivity  of  propellant, 
~  0.25  Btu/hr(ft2)(deg  F/ft) 
X  =  Reference  plane  located  at 
center  of  specimen,  t/2  = 
0.1875  inch 


Specimen  Thickness 

2X1 


The  calculation  revealed  that  the  maximum  temperature  rise  for  the  47  cps, 
6^  strain  test  was  13-8  F,  as  the  maximum  temperature  at  the  center  was 


153-8  F. 


The  importance  of  such  thermal  gradients  to  the  time  or  number  of 
cycles  to  failure  during  such  tests  remains  unknown.  The  numbers  calcu¬ 
lated  above  seem  to  indicate  that  gradients  arising  during  5  and  15  cps 
(5  F)  would  be  negligibly  small.  It  seems  equally  likely,  however,  that 
the  temperature  difference  might  cause  measurable  differences  in  results 
obtained  from  47  cps  tests. 
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Figure  29*  Deterioration  in  Energy  Lose  Measured 
During  a  Biaxial  Strip  Vibration  Test 

DATA  COLLECTION 


Several  items  of  data  including  anbient  temperature  :n  the  condi¬ 
tioning  box,  platen  forces,  and  platen  displacements  were  collected,  all 
as  a  function  of 'tine,  during  each  test.  A  counter  on  c  console  indi¬ 
cated  the  number  of  cycles  which  was  also  recorded  as  a  'unit  ion  of  time. 
The  time  the  crack  first  appeared  and  the  observed  rate  at  which  it  grew 
were  also  recorded.  Most  of  the  data  were  collected  by  frequent  observe-., 
tions. 


The  environmental  cabinet  control  operates  from  a  feedback  system 
that  provides  the  ambient  temperature  necessary  to  maintain  the  surface- 
mounted  thermocouple  at  a  constant  temperature  (140  F).  A  thermometer  was 
used  to  monitor  the  ambient  temperature  of  the  box.  The  difference  between 
ambient  and  surface  temperatures  is  a  measure  of  energy  or  heat  loss  from 
the  specimen.  Ambient  temperature  was,  therefore,  recorded  at  each  ob¬ 
servation. 

Platen  forces  and  displacements  were  measured  by  load  cells  in  the 
stationary  platen  support  and  an  LVDT  on  the  moving  platen  ram,  respec¬ 
tively.  Output  voltages  were  fed  to  an  oscilloscope  where  their  magni¬ 
tudes  could  be  displayed  separately  as  a  function  of  time  or  compared  to 
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provide  the  so-called  Lissajous  ellipse  on  a  force-displacement  plot. 

These  traces  vere  photographed  periodically  through  every  test,  indicating 
both  storage  and  loss  properties. 

Observations  for  cracking  vere  made  periodically  throughout  each  test. 
After  the  beginning  of  crack  growth,  the  point  termed  failure  in  this 
report,  more  frequent  observations  vere  made  to  correlate  crack  length 
with  the  nuab^r  of  cycles.  These  data  vere  noted  until  the  sequence 
caaera  was  installed  to  provide  a  photographic  record  of  crack  growth. 
Results  of  observations  of  crack  initiation  are  plotted  in  the  three 
graphs  of  Fig.  30*  Plotted  points  are  keyed  to  the  box  number  from  which 
speciasns  vere  cut  to  ascertain  whether  box-to-box  variations  vere  evident 
froa  these  tests.  They  vere  not.  Arrows  pointing  away  from  some  data 
points  indicate  that  the  exact  number  of  cycles  at  which  the  crack  began 
is  unknown  but  ie  bounded  .m  the  direction  indicated. 

Actually,  33  tests  rather  than  the  27  shown  in  T'lble  XII  were  suc¬ 
cessfully  completed  and  are  plotted.  The  extra  tests  were  performed, 
primarily  at  the  5  cps  frequency,  in  an  effort  to  narrow  the  spread  of 
the  scattered  data. 

The  ordinates  of  Fig.  30  are  actually  the  sum  of  the  vibration  dis¬ 
placement  amplitude,  uQ,  and  the  superposed  constant  displacement,  up 
(see  Eq.  9k).  It  was  intended  that  they  be  equal  in  all  tests.  However, 
as  shown  in  Table  XIII,  they  varied  slightly  from  this  intended  loading. 

After  actual  separation  of  the  specimen,  the  fracture  surface  was 
carefully  examined  to  determine  whether  the  mechanism  of  failure  could 
be  detected  visually.  Following  only  one  or  two  tests,  it  was  possible 
to  link  the  failure  to  an  actual  flaw  or  void.  Photographs  of  an  example 
flaw  from  which  a  failure  is  believed  to  have  originated  are  shown  in 
Fig.  31. 

CAMERA  COVERAGE 

A  70  am  Hulcher  sequence  camera  was  installed  to  monitor  crack 
initiation  and  growth,  (it  was  necessary  to  replace  the  original  door 
of  the  conditioning  box  with  a  fully  transparent  plexiglass  door  to 
facilitate  the  use  of  the  camera.)  The  shutter  of  the  camera  opened 
after  a  pre-set  time  interval,  and  the  associated  strobe  light  fired  when 
the  displacement  reached  the  next  maximum  value.  Consequently,  any  cracks 
present  in  the  specimen  were  in  the  strained  (open)  position  as  the 
shutter  opened. 
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TABLE  XIII.  SUMMARY  OF  THE  DISPLACEMENT 
BOUNDARY  FATIGUE  TESTS 


Temperature  -  140 

F 

u(t)  -  up  -  uo 

sm  <j>t 

Test 

No. 

Frequency, 

cps 

Nominal 

Strain, 

* 

Constant 
Displacement, 
Up*  iD- 

Amplitude, 

V 

in. 

u  ♦  u 

D  O 

Initial 
Crack  inf, 
Cycles 

SpeciBKn 

Separation, 

Cycles 

Box 

K  umber 

2 

in. 

57 

15 

7 

0.09375 

0.1062 

0.0999 

117,500 

133,700 

15 

59 

5 

7 

0.09575 

0.09575 

0.09575 

315,600 

420,300 

14 

60* 

5 

6 

0.0760 

0.0813 

0.0786 

243,100 

257,100 

14 

61 

5 

8 

0.1050 

0.114 

0.1090 

14,000 

22,700 

14 

62 

15 

6 

0.0810 

0.0863 

0.08360 

590,000 

636,500 

14 

63 

15 

7 

0.0920 

0.0991 

0.0955 

87,000 

105,300 

14 

64 

15 

8 

0.1080 

0.1680 

0.138 

21,800 

24,700 

14 

65 

5 

6 

0.077 

0.077 

0.077 

- 

184,700 

16 

66 

5 

7 

0.0980 

0.1015 

0.099 

- 

82,100 

16 

67  ' 

5 

8 

0.1080 

0.125 

0.116 

17,800 

31,900 

16 

68* 

15 

6 

0.0815 

0.0850 

0.0832 

287,400 

321,200 

16 

69 

15 

7 

0.0930 

0.107 

0.100 

70,700 

80,300 

16 

70 

15 

8 

0.1120 

0.1176 

0.1145 

23,200 

26,500 

16 

71 

5 

6 

0.0845 

0.088 

0.086 

262,000 

315.700 

13 

72 

5 

7 

0.0955 

0.106 

0.100 

586,200 

767.700 

13 

73 

5 

8 

0.1110 

0.1217 

0.116 

56,000 

131,300 

13 

74b 

15 

6 

0.0825 

0.1465 

0.1145 

331,900 

392,900 

13 

75 

15 

7 

0.0915 

0.0969 

0.0922 

52,500 

68,700 

13 

76 

15 

8 

0.106 

0.106 

0.106 

42,800 

49,800 

13 

77 

5 

6 

0.078 

0.078 

0.078 

- 

166,800 

18 

84 

5 

6 

0.081 

0.081 

0.081 

250,000 

267,200 

3 

85 

5 

7 

0.090 

0.090 

0.090 

49,800 

56.400 

3 

86 

5 

6 

0.082 

0.082 

0.082 

172,100 

229,300 

3 

87 

5 

7 

0.090 

0.090 

0.090 

53,100 

67,200 

3 

88 

47 

8 

0.107 

0,100 

0.104 

- 

4,300 

14 

91b 

47 

6 

0.095 

0.095 

0.095 

16,500 

26,800 

14 

92 

47 

6 

0.082 

0.092 

0.087 

27,000 

31,900 

14 

93 

47 

6 

0.079 

0.079 

0.079 

46,000 

53,900 

13 

94 

47 

7 

0.095 

0.095 

0.095 

22,100 

27,900 

13 

95b 

47 

8 

0.113 

0.099 

0.106 

20,000 

22,000 

13 

96 

47 

6 

0.078 

0.085 

0.082 

26,200 

32,800 

16 

97c 

47 

7 

0.090 

0.108 

0.099 

8,200 

10,500 

16 

98d 

47 

8 

0.107 

0.107 

0.107 

4,800 

6,400 

16 

a.  Electrical  Shut  Dovn 

b.  Data  Corrected  for  Electronic  Load  of  Camera  on  Bead  Out 

c.  Flaw  in  Specimen  (Not  risible  at  outset) 

d.  Theraal  Control  Problems 
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Figure  30.  Cycles  to  Failure  in  Biaxial  Strips 
Following  Tensile  Displacement 
Vibration  Tests 
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Figure  31.  Photomicrographs  of  a  Small  Void 

From  Which  One  Biaxial  Test  Speci¬ 
men  Failed 

Photographic  records  were  made  only  during  the  last  third  of  the 
fatigue  testing  program.  Data  collected  on  film  are  summarized  in  the 
length  of  crack  vs  number  of  cycle  curves  in  Appendix  VIII. 

The  Hulcher  camera  is  not  necessarily  the  best  equipment  for  monitor' 
ing  the  initiation  of  a  crack.  However,  the  availability/schedule  im¬ 
pact  necessitated  its  use.  A  limiting  condition  of  the  70  mm  camera  is 
that  if  it  is  set  far  enough  away  for  the  entire  specimen  to  be  within 
the  field  of  view,  the  photograph  scale  is  so  large  that  small  initial 
cracks  (0.25  inch)  can  neither  be  seen  nor  measured.  This  problem  was 
solved  by  plotting  the  longer  era ck-vs -cycle  curve  then  extrapolating 
the  curve  backward  until  it  intersected  the  zero  axis.  The  number  of 
cycles  at  this  point  was  considered  the  number  of  cycles  to  failure. 
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ANALYSIS  OF  RESULTS 


Early  in  the  analysis  phase  of  the  program  it  was  established  that 
one  can  hope  to  predict  the  initiation  of  fracture  reliably,  only  if 
stresses  and  strains  can  first  be  reliably  calculated.  A  failure  predic¬ 
tion  cannot  bc>  better  than  the  stress  analysis  on  which  it  is  based. 

Therefore,  calculations  were  made  to  ascertain  just  how  close  analyt¬ 
ical  stress/strain  calculations  (based  on  linear  viscoelasticity)  came  to 
actual  measured  forces  in  the  biaxial  strip  under  oscillatory  displace¬ 
ment  loading.  Material  properties  determined  with  uniaxial  tests  were 
first  used  in  these  calculations  and  were  found  to  predict  a  quite  dif¬ 
ferent  force  response  from  that  measured  on  the  biaxial  strip.  This 
finding  led  to  the  exploration  reported  in  Appendix  IX  where  the  effective 
differences  between  properties  determined  using  uniaxial  specimens  are 
contrasted  with  those  using  biaxial  specimens.  It  was  concluded  that  it 
would  be  best  to  use  biaxially  determined  properties  for  all  failure  pre¬ 
dictions  concerning  biaxial  strip  specimens. 

Figure  32a  and  b  compare  predicted  and  measured  vibration  responses 
shortly  after  the  start  of  a  test  and  after  50,000  cycles,  respectively. 

Interpreting  these  propellant  fatigue  data  and  the  extension  to 
other  loadings  and  sequences  of  loadings  requires  that  they  be  fitted 
to  the  equations  for  the  criticality  of  a  flaw.  The  energy  balance  ex¬ 
pression  for  the  initiation  of  the  fracture  of  a  spherical  flaw  following 
a  loading  of  this  type  has  been  derived  by  Jones,  Wagner,  and  Williams 
(ll).  The  derivation  incorporates  a  viscoelastic  relaxation  modulus  ex¬ 
pressed  in  the  form  of  a  Prony  series.  The  analogous  expression  for  a 
cylindrical  flaw  in  a  plane  stress  (very  thin)  body  has  been  formulated 
as  shown  in  Table  II.  That  expression,  probably  most  applicable  to  the 
biaxial  specimen, is  written 
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The  right-hand  side  of  this  expression  reflects  an  energy  buildup,  with 
time,  under  the  influence  of  the  two  elemental  loadings  of  Eq.  97?  the 
instantaneous  step  displacement,  Up,  and  a  sinusoidally  vibrating  dis¬ 
placement  of  amplitude,  u0.  The  material  properties,  namely  the  tensile 
modulus,  enters  through  the  use  of  coefficients,  E^  and  Ti>  chosen  to 
give  the  best  fit  to  the  measured  relaxation  modulus.  The  surface  temper¬ 
ature  of  the  specimen  was  held  constant  at  140  F  during  all  the  tests.  It 
was  assumed  that  it  was  permissible  simply  to  shift  the  property  data  in 
accordance  with  the  biaxial  time-temperature  equivalence  factors  reported 
in  Appendix  V. 

A  plot  of  the  calculated  value  of  the  right-hand  side  of  Eq.  99  vs 
time  for  the  three  loading  frequencies  of  the  Rocketdyne  tests  are  shown 
in  Fig.  32.  The  frequencies  1.288,  0.412,  and  0.1374  at  77  F  are  equiva¬ 
lent  to  47,  15,  and  5  cps  at  140  F.  The  time  scale  (the  value  oi  tf  in 
Eq.  99)  is  unchanged  by  temperature  and  represents  laboratory  time. 

The  hypothesis  is  that  when  the  value  of  the  right-hand  side  becomes 
just  equal  to  the  constant  of  the  left-hand,  side  the  flaw  begins  to  grow. 
The  constant  value  assigned  depends  upon  the  geometry  and  displacement 
magnitudes  and  is  directly  dependent  upon  the  characteristic  strain  energy 
release  rate,  y,  divided  by  the  initial  flaw  radius,  a0.  It  was  hoped 
that  the  values  of  these  parameters,  determined  in  the  static  tests  de¬ 
scribed  in  the  preceding  section,  could  be  used  to  calculate  the  constant. 
For  example,  if  one  uses  the  typical  values 

f-  -  50  ;  k  =  0.5  ;  ~  =  14.3 

Cl  ,  U 

0  o 

then  the  value  of  the  constant  is 

8^~  (1  +  3k2)2  ~  2  =  3.9  x  103  (101) 

o  o  in. 

But  this  value  is  much  too  small  to  compare  to  the  curves  of  Fig.  33* 

At  least  one  reason  for  this  problem  was  the  assumption  of  an  instan¬ 
taneously  applied  displacement,  up  (=  the  amplitude,  u  ,  in  these  tests), 
which  leads  to  the  large  term  on  the  right  hand  side  of  Eq.  100.  This 
single  term  is  precisely  the  term  for  the  critical  energy  level  during  a 
step  displacement  loading  without  the  sinusoidal  vibration.  This  led  to 
the  idea  of  replacing  the  term  with  the  equivalent  expression  for  criti¬ 
cality  following  a  ramp  loading  history.  That  is,  replace  the  term 
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E  with 

E 

of  the  reap 
of  a  cycle, 
times,  with 
in  Appendix 


(2) 

2  Er#j  (tj)  where  t^  ia  the  reap  duration  time.  An  estimate 


tiaea  can  be  made  by  assuming  their  duration  was  one-quarter 
This  leads  to  the  ramp  times  shown  in  Table  XIV.  From  these 
the  plot  of  the  biaxial  relaxation  modulus  and/i^s  integrals 
IX,  one  can  evaluate  the  critical  energy,-^ 


id  .its  int< 
E(2)(tl). 


TABUS  XIV.  APPROXIMATE  RAMP  LOADING  TIMES 
AND  THE  CORRESPONDING  CRITICAL 
ENERGY  LEVELS  AT  77  F 


Reduced 

Frequency, 

u/aT, 

cps 

Ramp  Loading 
Time,  t*, 
min 

.  2  rel 

psi 

(t.), 

1.288 

3.24  x  10"3 

2800 

0.412 

10.1  x  10'3 

2270 

0.1322 

31.5  x  10~3 

1810 

The  curves  of  Fig.  33  are  lowered  significantly  when  the  critical 
energy  corresponding  to  a  ramp- type  loading  is  used  in  place  of  that  for 
a  step  loading.  These  calculations  resulted  in  the  curves  of  Fig.  Jk, 
However,  only  a  quick  calculation  indicates  that  the  parameters  determined 
during  biaxial  static  tests  do  not  fit  the  vibration  fatigue  data. 

As  a  consequence,  earlier  parameters  were  completely  disregarded; 
and  new  values  were  chosen  to  give  a  reasonable  fit  to  the  test  data. 

The  predicted  failure  lines,  for  47,  15,  and  5  cps,  at  140  F  are  compared 
to  the  measured  failure  points  in  Fig.  35-  Several  interesting  observa¬ 
tions  can  be  made  from  Fig.  35.  First,  although  the  data  show  the  wide 
scatter  characteristic  of  fatigue  failure  tests,  there  is  still  a  rather 
definite  indication  of  a  slope  of  about  -l/4  on  the  log-log  plot.  (Broken 
lines  indicating  this  trend  were  not  calculated  numerically.  This  is  in 
direct  contrast  to  theoretical  predictions  which  have  a  slope  of  about 
-l/2. 
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T  =  140  F  /  ' 

W  =  94.5  rad/sec  f 

Specimen  Dimensions:  / 

13— l/4  x  3/8  x  2.8-in  long  /  I 

Loading:  '  ' 

u(t)  =  0.09374  (l  +  sin  94. 5t)  /  I 

Measured  ^  y/^  /l 

/  /: 

^  /  j' — Predicted 

//  ■  /  1  /  1  1 

//  0.05  /  0.1  '  0.15  0.2 

‘  f  /  / 


Displacement,  in. 


a.  After  a  few  cycles  of  tbe  test  (<30) 


Figure  32.  A  Comparison  of  the  Biaxial  Force-Displacement 
Curves  as  Measured  and  Predicted  (using  the 
biaxial  characterization). 


T  =  140  F 

W  =  94.5  rad/sec  Predicted 


b.  Alter  50,000  cycL  of  the  test 
Figure  32.  Concluded 


100 


Figure  35.  Fatigue  Bata  Collected  by  Subjecting  Biaxial 
Strips  to  Vibrating  Displacements  of  the  Form 
u(t)  =  u  +  u  sin  cot 

'  A  A 


5  cps 


47  cps 


47  cps 


*  Tests  were  interrupted 


Test  #99* 


Test  #100 


Test  #101 


Test  #103 


Test  #104* 


Test  #105 


Test  #106 


Figure  36.  Summary  of  the  Dual  Amplitude,  Constant  Frequency 
Vibration  Fatigue  Tests  (Biaxial  Strip  Specimens) 
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10,000  Cycles 
47  cps 


51,700  Cycles 
5  cps 


wiry  A. 


Test  No.  112 
Temperature  =  140  F 


Figure  37-  Summary  of  the  Dual  Freque 
Constant  Amplitude  Vibrati 
Fatigue  Tests  (Biaxial  Str 
Specimens) 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  1-year  program,  conducted  in  association  with  three  other  con¬ 
tractors,  and  coordinated  by  the  AF  Rocket  Propulsion  Laboratory,  opened 
for  public  view  the  limitations  and  capabilities  of  current  grain  stress/ 
strain  and  failure  analyses.  The  concensus  must  surely  be  that  there  are 
more  limitations  than  capabilities. 

The  early  decision  by  the  Air  Force  to  use  a  combine  of  companies 
exploiting  the  particular  specialties  of  each  proved  to  be  a  wise  one. 

Due  to  frequent  coordination  conferences,  close  and  interested  Govern¬ 
ment  surveillance  and  a  spirit  of  free  exchange  of  data  (both  favorable 
and  unfavorable)  among  the  several  researchers,  maximum  use  and  benefit 
were  derived  from  the  volumunious  quantity  of  data  and  information  gen¬ 
erated.  Surprisingly  few  contradictions  were  found  in  test  results 
originating  from  the  different  organizations,  all  of  which  adds  to  the 
significance  of  the  conclusions  which  can  be  made. 

Two  salient  conclusions  are  most  obvious  and  important.  The  first 
is  that  if  the  stresses  or  even  the  strains  cannot  be  calculated  to 
within  acceptable  engineering  accuracy,  then  any  attempt  to  meaningfully 
predict  failure,  regardless  of  the  loading  history  represents  merely  an 
analytical  exercise.  Accordingly  the  linear  viscoelastic  analyses  based 
on  the  elemental  properties  provided  and/or  derived  proved  inadequate 
even  for  the  most  elementary  geometries  and  loadings.  The  variation  in 
the  classical  properties  due  to  state  of  stress  and  the  type  of  loading 
were  clearly  shown  by  the  data  on  several  of  the  tests.  Neither  these 
variations  nor  their  causes  are  currently  understood.  However,  their 
relative  importance  to  the  working  stress  analyst  must  surely  provide 
motivation  for  thoir  further  study.  The  second  is  that  environmental 
effects  on  propellant  surfaces  are  important,  perhaps  to  the  point  of 
overshadowing  or  masking  the  effects  of  accumulated  damage  due  to  prior 
structural  loadings.  The  susceptibility  of  the  CTPB  propellant  failure 
properties  to  change  as  a  result  of  the  environment,  especially  high 
humidity,  is  real  and  must  be  taken  into  account  in  motor  grain  analyses. 

A  multitude  of  specific  conclusions  are  evident  from  the  Rocketdyne 
portion  of  the  program.  Some  concern  the  testing  procedures  and  are 
listed  only  to  assist  other  experimenters  working  with  similar  specimens 
or  tests.  Others  are  more  general  and  will  serve  as  a  guide  for  those 
attempting  to  predict  failures  using  the  energy  balance  techniques. 
Ordered  somewhat  by  importance  they  are: 
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1.  The  criticality  equations  for  the  growth  of  a  spherical 
or  cylindrical  flaw  in  a  viscoelastic  media  can  be  de¬ 
rived  for  a  wide  variety  of  practical  loading  histories. 
These  equations  at  heat  can  be  used  to  predict  the  load¬ 
ing  conditions  which  will  lead  to  propellant  cracking, 
at  worst,  can  be  used  to  evaluate  the  significance  of 
the  influential  parameters  on  propellant  cracking. 

2.  The  triaxial  poker  chip  test  results  fit  the  predictions 
of  the  viscoelastic  fracture  equations  better  than  those 
from  any  of  the  various  experimental  tests.  This  is 
reasonable  because  the  state  of  stress  closely  simulates 
that  hypothesized  lor  the  derivations  of  the  spherical 
flaw  equations  and  because  the  problem  with  surface 
effects  is  essentially  eliminated. 

3.  The  time  to  failure  in  the  vibration  fatigue  tests,  like 
the  creep  tests,  is  an  extremely  sensitive  measure  and 
thus  exaggerates  the  statistical  spread  of  failures. 
However,  it  is  important  to  the  analyst  that  this  varia¬ 
tion  be  quantitatively  understood. 

4.  A  reasonable  amount  of  evidence  became  available  to 
indicate  that  the  characteristic  strain  energy  release 
rate  for  propellants  is  a  function  of  both  loading 
rate  and  temperature. 

5.  During  the  vibration  fatigue  tests,  the  propellant 
response  properties  change.  This  change  must  be  taken 
into  account  if  the  crack  initiation  is  to  be  predicted. 

6.  The  drop  off  in  the  center  load  cell  output  which  has 
long  been  observed  to  precede  the  ultimate  poker  chip 
loading,  during  constant  rate  tests,  is  indicative  of 
the  formation  and  growth  of  small  scattered  internal 
flaws. 

7.  The  cumulative  damage  testing  machine  can  be  used  to 
subject  biaxial  propellant  strips  to: 

a.  extremely  slow  stress  or  strain  rates 

b.  extremely  fast  stress  or  strain  rates 

c.  simultaneous  preprogrammed  straining  and  cooling 
histories 
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d.  bilinear  strain  rates 

e.  vibration  fatigue  tests  (either  force  or  displace¬ 
ment  controlled) 

8.  Premature  failures  near  the  edges  and  bond  lines  of  the 
biaxial  strip  specimens  can  be  precluded  by  careful 
reinforcement  with  a  soft  layer  of  EC  1949- 

These  conclusions,  while  generally  encouraging,  remain  nebulous  and 
unsatisfying  to  the  practicing  engineer  calculating  the  margin  against 
failure  of  his  particular  motor  design.  For  this  reason  it  is  important 
that  careful  consideration  be  given  to  the  direction  that  similar  or 
related  future  activities  might  take.  The  great  expense  of  the  carefully 
controlled  laboratory  tests  of  the  completed  program  dictates  that  their 
full  meaning  be  gleaned  and  presented  in  a  usable  and  concise  format. 

Such  a  format  would  certainly  include  formulas,  graphs,  tables,  and  ref¬ 
erences  to  the  more  detailed  origins  of  the  information.  The  importance 
of  the  statistical  distributions  to  the  reliability  of  the  resulting 
calculations  should  be  quantitatively  evaluated. 

Further  testing  activities  must  be  promoted  to  determine  more  precise 
information  concerning  surface  effects  and  environment  upon  the  propel¬ 
lants  susceptibility  to  cracking  failure.  It  is  desirable  that  further 
tests  designed  to  correlate  more  directly  with  actual  motors  and  STV's  be 
designed  and  perfected.  This  is  especially  true  for  loadings  involving 
thermal  changes. 

But  of  first  order  importance  is  the  necessity  to  devise  a  technique 
adequate  for  calculating,  from  laboratory  characterizations,  meaningful 
stress  and  strain  magnitudes  in  any  stress  field  following  practical 
loading  histories. 
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APPENDIX  I 


FINITE  ELASTICITY  AND  FOAM  MECHANICS 
by  P.  J.  Blatz 


FOAM  MECHANICS 


When  a  thin  disc  (poker  chip)  of  propellant  is  subjected  to  par¬ 
allel  plate  loading  and  held  at  some  fixed  strain,  the  bonded  specimen 
gradually  relaxes  the  applied  load  and  simultaneously  dewets  and  becomes 
a  foamed  or  bodied  structure.  This  latter  phenomenon  was  actually 
observed  by  X-ray  photography  (see  Poker  Chip  Tests  and  Analysis, 

Section  III), 


After  dewetting,  the  foamed  material  should  evince  a  stress-strain 
behavior  characteristic  of  a  foamed  material.  This  should  imply,  for 
example,  a  value  of  Poisson's  ratio  equal  to  l/4.  This  latter  statement 
is  consistent  with  the  observation  that  after  relaxation  ensues  and 
while  the  material  is  dewetting  the  peak  stress  2  P/A  gradually  relaxes 
to  ?/A.  Accordingly,  a  constitutive  equation,  first  used  by  Ko  and 
Blatz  (30)  was  postulated  to  represent  the  behavior  of  this  material, 
name  ly: 


where 


X  =  Extension  ratio 
J  =  Volume  ratio,  V/Vo 
G  =  Shear  Modulus 

=  True  stress,  ith  component 

The  problem  of  a  spherical  shell  with  a  small  spherical  flaw  at 
its  center  was  solved  for  this  particular  constitutive  equation.  In 
the  meanwhile,  a  triaxial  stress-stretch  curve  was  procured  on  the 
foamed  material.  When  the  data  (Fig.  38)  became  available,  it  was 
noted  that  the  range  of  strain  over  which  the  material  was  tested  was 
only  a  few  per  cent,  leading  to  the  conclusion  that  it  was  no  longer 
necessary  to  approach  the  characterization  in  terras  of  finite  strain 
theory,  but  rather  that  Eq.  104  could  take  the  form: 
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(105) 


where 


=  Dilatation  ~  +  e2  +  e3 


=  ith  component  of  strain 

which  corresponds  to  Hooke's  law,  with  v  =  l/4 

Furthermore,  once  the  material  is  foamed,  deformation  between  the 
parallel  plates  is  purely  uniaxial  rather  than  hydrostatic.  Therefore: 


ex  =t?  -  c 


e2  -  e3  .  0 


(106) 

(107) 


(108) 


tm* 


(109) 


From  the  slope  of  the  data  in  Fig.  38,  it  is  concluded  that  the 
value  of  the  sheer  modulus  of  the  foamed  material  is  ~3500  psi.  Since 
this  value  is  significantly  higher  than  that  of  the  shear  modulus  of 
the  undewetted  material,  the  assumption  that  the  material  has  completely 
and  uniformly  dewetted  is  obviously  erroneous. 

ka  alternative  procedure  is  to  return  to  the  spherical  shell  model 
(i.e.,  assume  that  the  dewetting  is  highly  localized).  Under  this  as¬ 
sumption,  the  effective  stress-strain  behavior  of  the  spherical  shell 
is  given  by: 


(no) 


where 

a  -  Flew  Radius 

b  -  Outer  Radius  of  Spherical  Model 

If(l-2v)  ia  now  neglected  by  the  assumption  that  v"*0.5  in  the  numerator, 
Eq.  110  becomes: 
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(m) 

or 

(112) 

Taking  G  to  be  1000  psi,  which  is  typical  of  Lockheed's  STV  propellant 
0064-61E,  leads  to  the  assumption  that  about  3/1  of  the  material  has 
voided,  i.e.: 


b3 


(113) 


X-ray  photographs  indicate  that  this  figure  is  varying  with  time  as 
would  be  predicted  using  Eq.  110  if  values  for  G  and  v  were  varying  with 
tine . 

GRIFFITH  CRACK  THEORY 

After  Griffith  proposed  that  the  critical  condition  for  the  growth 
of  a  crack  can  be  determined  by  a  power  balance  between  the  rate  at  which 
mechanical  work  is  done  on  a  body  and  the  rate  at  which  this  work  is 
taken  up  by  an  increase  in  strain  energy  together  with  an  increase  in 
surface  energy,  cracks  of  various  geometries  in  bodies  of  various  shapes 
were  analyzed  according  to  this  criterion.  In  such  calculations,  two 
precautions  must  be  observed.  First,  the  integral  of  the  strain  energy 
density  over  the  body  must  exist.  This  can  be  ensured  by  using  a  super¬ 
position  principle  in  which  one  adds  together  two  problems,  one  of  which 
involves  zero  loading  at  the  outer  boundary  but  locally  variable  through¬ 
out  and  the  other  involves  finite  loading  at  the  outer  boundary  and 
locally  constant  stress  throughout.  The  latter  problem  is  trivial  and 
cau  be  added  to  the  solution  of  the  first  problem  to  give  a  desired  re¬ 
sult.  The  first  problem  is  usually  one  involving  pressure  in  the  crack. 
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The  latter  problem  is  trivial  and  can  be  added  to  the  solution  of  the 
first  problem  to  give  a  desired  result.  The  first  problem  is  usually 
one  involving  pressure  in  the  crack. 

The  second  precaution  deals  with  the  question  of  whether  the  power 
balance  applied  to  the  infinite  body  with  an  infinitesimal  crock  gives  the 
same  result  as  the  power  balance  applied  to  the  finite  body  with  a  finite 
crack,  after  which  the  double  limit  is  taken  as  the  outer  boundary  moves 
to  infinity  and  the  finite  crack  becomes  a  flat  or  infinitesimal  crack. 

The  answer  to  this  question,  shown  in  the  following  section,  is  clearly 
no.  The  Griffith  stress  is  derived  for  a  finite  body  between  confocal 
ellipses.  First,  when  correctly  done,  the  result  is  independent  of  the 
position  of  the  outer  boundary  which  can  then  be  taken  to  be  infinite. 
Second,  in  the  limit  as  the  elliptical  crack  becomes  a  flat  crack,  the 
Griffith  stress  becomes  infinite. 

IMPORTANCE  OF  FINITE  ELASTIC  EFFECTS  IN  THE  INSTABILITY 
OF  A  SPHERICAL  FLAW  IN  A  THICK  SPHERICAL  SHELL 

In  stretching  a  thin  disc  bonded  to  two  rigid  parallel  plates  in  the 

limit  of  infinite  diameter-to-thickness  ratio  (so  that  the  strains 

c  =  e  =  0),  infinitesimal  axial  stress  is  given  by: 
a  y 

o  «  (K  +  ~  G)e  (114) 

Z  V  3  Z 


where 


K  =  Bulk  Modulus 
Since,  for  rubberlike  materials 

K  »  G  (115) 


the  stress  field  is  very  nearly  hydrostatic.  On  this  basis,  it  is  con¬ 
venient  to  imagine  a  spherical  shell-like  region  between  the  two  plates, 
at  the  center  of  which  is  initiated  a  spherical  flaw  that  eventually 
grows  to  produce  the  fracture  (see  sketch  below). 


a  =  radius  of  hole 
D  =  diameter  of  plate 
b  =  outer  radius  of  sphere 


t\ 
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CO 
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(K  +  4/3G)  e  ~  Ke,  .  .  hydrostatic 
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Elaatost^  cic  Analva  is  of  Parallel -Plate  Tenaion  Data 

For  thia  geometry  it  is  possible  io  do  an  exact  general  analysis  of 
the  instability  criterion,  for  \.hich: 

•7* **  >  W  +  K  (ll6) 

A  —  a  a 


where 


A  =  Unde formed  Radius  of  Spherical  Flaw* 

V  =  Strain  Energy  +  Dissipation  Rate  Density 

Y  =  Characteristic  Strain  Energy  Release  Rate 

Since  Eq.  116  shows  that  surface  energy  depends  directly  on  the 
strain  energy  density  at  the  flaw  interface,  it  might  be  expected  that 
the  surface  energy  calculated  from  finite  elastic  theory  would  differ 
significantly  from  that  calculated  from  infinitesimal  elastic  theory. 

This  statement  is  poignant  because  it  is  precisely  at  the  flaw  boundary 
that  large  displacement  effects  are  encountered.  For  example,  it  will 
not  be  uncommon  for  the  deformed  radius  to  be  several  times  larger  than 
the  unde  formed  radius,  which  means  a  strain  of  several  hundred  per  cent. 
Therefore,  y  vas  calculated  from  both  theories.  Kinetic  energy  effects 
were  not  included  in  the  calculation.  For  moderately  slow  strain  rates, 
the  inertial  term  is  5 — 7  orders  of  magnitude  less  than  the  strain  energy 
term . 


Assuming  that  a  strain  and  stress  at  the  point  of  instability  (growth 
•  initiation)  of  the  spherical  flaw  are  provided  and  taking  the  outer 

radius  of  the  spherical  shell  to  be  0.5  the  disc  thickness,  the  two  impor¬ 
tant  material  properties,  A/B  and  y/B,  can  be  calculated?*  Pertinent 
equations  are: 


*The  capital  A  in  this  equation  corresponds  to  a0  as  used  in  the  body 
of  this  report. 

**B  is  the  undeformed  outer  radius  of  the  spherical  model. 
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In  evaluating  these  quantities,  various  values  of  the  Young's 
modulus  were  chosen  to  show  how  sensitive  the  computed  values  of  y  and  A 
are  to  that  parameter.  Tables  XV  and  XVI  show  lie  results  for  a  Neo- 
Hookean  and  a  Hookean  material,  respectively.  When  E  is  500  psi,  which 
corresponds  to  the  actual  test  material,  the  y  calculated  from  Hookean 
behavior  is  only  30^  lower  than  that  calculated  from  Neo-Hookean  behavior; 
whereas  A  is  overestimated  by  a  factor  of  2.  This  says  that  in  a  tri- 
axial,  nearly  hydrostatic,  state  of  stress,  geometrical  nonlinearities 
are  not  important  in  an  incompressible  material.  On  the  other  hand, 
further  attention  will  be  paid  to  the  effect  of  geometrical  nonlinearities 
in  uni-  or  biaxial  states  of  stress.  As  an  expediency,  the  linear  theory 
will  presently  be  applied  to  these  problems  of  instability. 


TABLE  XV.  SOLUTIONS  USING  NEO-HOOKEAN  FORMULAS 


ab = 

225  psi 

B  =  0.3  inch 

Sj,  =  0.00861 

E,psi 

270 

300 

313 

337.5 

373 

450 

500 

“rib 

2.50 

2.25 

2.14 

2.00 

1.80 

1.50 

1.35 

*a 

CO 

8.00 

5.60 

4.02 

2.88 

2.06 

1.81 

V 

0 

0.0000505 

0.000148 

0.000404 

0.00108 

0.00297 

0.00434 

A,  in. 

0 

0.00437 

0.00626 

0.00913 

0.0121 

0.017 

0.0192 

Y,  in. -lb/ in. 2 

03 

14.1 

10.2 

7.81 

5.33 

3.65 

3.01 

ffija 

G 

00 

64.0 

31.4 

16.1 

8.28 

4.18 

3. 18 

119 


*’4B.  f!  VVI 


SOLUTIONS  Os  INC,  -:»f:  T>'r  FORMULAS 


4 


E,  psi 

270 

300 

315 

337.5 

375 

450 

500 

0rllA 

2.50 

2.250 

2.140 

2.00 

1.80 

1.50 

1.350 

ea 

0.625 

0.563 

0.535 

0.500 

0.450 

0.375 

0.338 

70^  low 

V 

0.0*38 

0.015 

0.0161 

0.0172 

0.0191 

0.0230 

0.0255 

6  fold 

A,  in. 

Y,  inrib 

0.0283 

0.0292 

0.0293 

0.0304 

0.0315 

0.0335 

0.0347 

in. -lb 

in.2 

3.05 

3.35 

2.77 

2.66 

2.48 

2.2 

2.02 

{’  .  2 
in. 

!ak 

G 

3.75 

3.178 

.3.21 

3.00 

2.70 

_ 

2.25 

2.03 

50^  low 

FOAM  MECHANICS 


When  the  parallel-plate  tensile  test  under  constant  load  was  con¬ 
ducted,  it  vas  observed  that  the  disc-shaped  specimen  gradually  became 
full  of  voids  from  25— 50  mils  in  size.  For  a  foamed  material,  Blatz  (30) 
has  shown  that  a  useful  stress-strain  law  is  given  by: 


Vi 

G 


+  J 


4(l-f) 


2(l-f) 

l-2v 


+ 


1 

,l-2v 


2(l-f)J 


(117) 


This  relation  provides  a  very  precise  and  useful  characterization, of 
propellant  in  the  foamed  state  (i.e.,  after  dewetting).  \! i ch  v  =  l/ 4 , 
which  seems  to  fit  the  data  best, 

-  V +  J  (1_4f)  +  2fj2  (ll8) 

Xi 

For  this  particular  choice  of  Poisson's  ratio  and  for  various  values 
of  f,  the  relation  between  load  and  extension  in  tbe  poker  chip  has  been 
calculated  using  the  flat  model. 
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;7W»fl  VISCOELASTIC  CUMULATIVE  DAMAGE 

It  has  been  su*»*3ted  that  cumulate  damage  may  be  predicted  aa 
the  viscoelastic  creep  oi  a  suitably  mod?  ed  flaw  in  a  specimen  of  pro¬ 
pellant  subjected  to  vibratory  load  or  d:. placement  input.  This  problem 
was  fires  t  considered  for  a  vibratory  displacement  input  at  the  outer 
boundary  '*  a  spherical  shell  with  a  spherical  flaw. 

Neglecting  kinetic  effects  and  using  the  convolution  integral  to 
calculate  time-dependent  stress  resulting  from  a  time-dependent  strain, 
Eq.  116  becomes 

ia.tfLf  ieb<t'>  f*  W*  Lw) 


enr 


E„1(*-t')dCb(t') 


(120) 


for  linearly  viscoelastic  materials. 


Equation  119  reveals  that  if  the  quantity  in  brackets  increases 
with  tine,  a  (the  flaw  radius)  must  also  increase  to  maintain  the  product 
constant  and  equal  to  yA*  From  Eq.  120,  then,  associated  stress  for  a 
given  displacement  input  can  be  calculated.  The  question  of  the  behavior 
of  the  quantity  in  brackets  for  a  specified  displacement  input  then  arises. 

To  evaluate  the  integral  for  a  chosen  sinusoidal  input,  the  following 
behavior  of  the  relaxation  function  is  noted.  First,  propellant  data  are 
well  represented  by  the  form 

Er.l  '  Er  +  <E«  -  Er>  '  Er  <  Er.l  <  V*”  (Wl)  ' 

The  preceding  representation  is  actually  an  asympototic  expansion  of 
the  complete  function  known  as  Mittag-Leffler  (3l).  Data  are  fitted  with: 

n  ■  0.194 

a  -  0.00151  (t  in  minutes) 


With  this  representation; 


-  er  +  Mr>  rffej  kfe  <122> 


co  that 

~E  .  s  E^°\  for  all  k. 
re  I 

With  this  mild  approximation,  Eq.  119  may  be  simplified.  In  doing 
this,  the  integrals  must  be  separated  into  two  parts  to  remove  the  dis¬ 
continuity  associated  with  the  initial  strain  input.  In  carrying  out  the 
integration  by  parts,  no  assumption  is  made  that  Jt  is  large,  so  the  result 
reduces  to  the  correct  result  as  t  ■*  o  and  is  the  first  term  of  an  ordered 
expansion  that  may  be  made  as  exact  as  possible  by  including  higher  terms. 

Choosing  the  vibratory  input  m  the  form  gives: 

%  =  Cbm  ■1^U>t  ,  0  <  a  <  1  (123) 

whereupon  Eq.  119  becomes: 


3?  4 


2\^g  +  (^8a^  +  Cl  ^ UUt)  E^^t) 

"/  (Hu)2 


The  period  of  the  vibratory  input  is  given  by: 
uIT  =  2tt 

and  the  time  to  a  peak  is  given  by: 
t  =  (N  +  t)T 


(124) 


(125) 


(126) 


so  that 


t«t  =  2ttN  +  ^ 


(127) 
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Insertion  of  this  result  in  Eq.  124  yields: 


where 


N  ■  Number  of  cycles 


As  N  ->  0,  Eq.  128  reduces  to: 


Y«5 

2b 


-?  >  E 

6-  g 


(128) 


(129) 


which  is  the  correct  result  for  a  step  input  into  an  elastic  material. 

On  the  other  hand,  Erei  is  a  monotonically  decreasing  function;  so  Eq.  119 
predicts  that  the  flaw  radius  must  decrease  with  the  number  of  cycles. 

This  means  that  it  achieves  its  peak  on  the  first  cycle  and  stays  constant 
thereafter.  This  is  to  be  expected  on  the  basis  of  thermodynamic  be¬ 
havior  of  a  dissipative  material.  One  can  easily  take  a  simple  situation 
in  which 


E 


rel 


(130) 


and 

6^(1  -  cos  uut) 

2 

In  this  case  there  is  no  discontinuity  in  the  input,  and  the  integrals 
are  exactly  calculable.  That  the  peak  flaw  radius  (ya^/b^)  must  decrease 
linearly  with  N  is  obvious.  This  again  means  that  it  achieves  its  peak 
value  on  the  first  cycle  and  remains  constant  thereafter. 

In  summary,  a  vibratory  displacement  input  of  fixed  amplitude  does 
not  produce  cumulative  damage  under  the  assumption  of  Griffith  instability. 
To  achieve  same,  one  must  input  a  load  history  which  is  ever  higher  than 
the  value  to  which  the  load  will  relax  in  the  case  of  constant  strain 
input.  Then  and  only  then  will  cumulative  damage  ensue.  Conversely,  a 
vibratory  displacement  input  of  fixed  amplitude  is  known,  de  facto,  to 
produce  cumulative  damage,  which  strongly  indicates  that  other  mechanisms 
besides  Griffith  instability  are  involved. 
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To  compute  the  strain  output  produced  by  a  specified  load  input, 
th*  flow  radius  between  Eq.  119  and  120  must  be  eliminated.  Therefore: 

-3/5 


c.-lh  =  3  X 


"S*4  >  Jo  E«l(t  -‘"K*1") 


(132) 


where  Eq.  119  has  been  neglected  with  respect  to(bVa^)»  Equation  132 
in  general  will  have  to  be  integrated  numerically  because  of  the  nature 
of  the  problem  involving  a  moving  boundary.  In  fact,  previous  statements 
to  the  effect  that  the  strain  output  for  a  given  load  must  be  cast  in 
terms  of  moments  of  the  creep  function  are  meaningless  since,  for  the 
spherical  flaw  (neglecting  inertia): 


3-4 


(133) 


(134) 


(135) 


4(1  -  fr) 
b-> 

Now  the  strain  energy  plus  dissipation  rate  density  is  given  by: 

w  =  £  (oV-tsyej)  -  2  .(*  (V°r)dEe(t  ) 


(136) 
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where  the  pride  designates  the  value  at  past  time  t1  and  where 


J  (t-t')df(t') 
crp  ' 


(138) 


In  general,  Eq.  137  cannot  be  integrated  without  additional  input  since 
the  motion  of  the  boundary  as  expressed  by  a  (t,')  depends  not  only  on 
the  load  history  but  also  on  the  history  of  the  boundary. 


In  summary,  Eq.  119  and  120  are  the  only  two  equations  by  which 
relations  between  load,  displacement,  and  spherical  flaw  radius  can  be 
expressed . 
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APPENDIX  II 

PHOTOMICROSCOPIC  DOCUMENTATION 
OF  INTERNALLY  PRESSURIZED  CYLINDERS 
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APPENDIX  II 

PHOTOMICROSCOPIC  DOCUMENTATION  OF 
INTERNALLY  PRESSURIZED  CYLINDERS 

by 

J.  E.  Ifilzinger,  Rocketdyne  Research 

Investigations  of  the  mechanisms  resulting  in  the  initiation  and 
propagation  of  fracture  in  materials  are  continuous  attempts  to  reduce 
the  uncertainties  associated  with  failure  analyses.  The  analytical 
prediction  of  behavior  of  actual  structures  fabricated  frov,  come  classes 
of  these  materials  (i.e.,  brittle,  composite)  has  not  been  satisfactorily 
developed,  and  realistic  behavior  descriptors,  either  qualitative  and/or 
quantitative  are  needed  as  a  basis  for  practical  engineering  judgements. 

This  report  describes  a  study  conducted  by  the  author  on  the  behavior 
of  one  such  material,  a  composite  propellant,  under  structural  and  envi¬ 
ronmental  conditions  translatable  to  service  applications. 

OBJECTIVE 

The  primary  purpose  of  this  effort  was  the  photomicroscopic  docu¬ 
mentation  of  the  structural  behavior  of  the  Lockheed-furnished  STV 
composite  propellant  under  known  stress  and/or  strain  states,  (internal 
pressurization  of  small  hollow  cylinders  at  preselected  levels  of  stress 
or  strain  rate  and  temperature  provides  outer  surface  phenomena  indi¬ 
cative  of  the  type  and  extent  of  microstructural  damage  associated  with 
a  highly  filled  elastomeric  system.) 

To  achieve  this  objective,  the  following  were  accomplished; 

1.  Photographic  documentation  of  the  surface  behavior  of 
botn  machined  and  cast  surface  specimens  during  a  lineal 
increase  in  stress  of  100  psi/min  to  failure  at  ambient 
laboratory  conditions 

.  A  highly  detailed  record  of  constant  stress  rate  sur¬ 
face  mechanisms  for  machined  and  cast  surfaces  at 
relatively  low  stress  rates  and  ambient  conditions 


PROCEDURE 

The  preparation  of  hollow  cylindrical  test  specimens  from  the  CTPB 
propellant  took  one  of  two  forms:  machining  from  a  cast  block,  or 
sectioning  of  a  molded  cylindrical  tensile  specimen.  A  summary  of  the 
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cylinder  dimensions  (Table  XVII)  revealed  that  all  test  specimens  had 
nominal  l/8-inch  wall  thickness  although  minor  web  variations  between 
and  within  cylinders  resulted  from  the  internal  boring  operation.  Two 
deliberate  changes  in  cylinder  geometry  evolved  during  testing.  The 
initial  outside  machined  diameter  of  0.8  inch,  although  more  conducive 
to  good  photographic  practice,  was  reduced  to  0.5  inch  to  provide  a 
direct  comparison  with  the  cast  surface  specimens  available  only  in  the 
smaller  size.  In  addition,  the  length-to-diameter  ratio  was  increased 
from  1.5/1  to  2/l  when  the  calculated  possibility  of  substantial  end 
effects  reaching  the  center  hoop  stress  area  was  confirmed  photograph¬ 
ically. 

Propellant  cylinders  were  closed  at  each  end  by  Plexiglass  fixtures 
grooved  to  accept  the  web  thickness  and,  where  applicable,  the  internal 
0.0045-inch. latex  bladder.  An  ambient  cured  epoxy  resin  wes  used  to 
oond  the  propellant  into  the  end  fixtures.  Eight  cylinder  assemblies 
were  satisfactorily  fabricated  in  this  manner  after  an  initial  failure 
using  an  inadequate  end  cap  design. 

The  entire  test  series  was  conducted  using  dry  gaseous  nitrogen 
both  as  a  controlled  storage  medium  and  as  the  pressurant  for  all  stress 
levels.  Photographic  techniques  employed  ranged  from  time  lapse  methods 
with  35  mm  to  16  mm  motion  at  128  or  400  frames/second. 

RESULTS  AND  DISCUSSION 

Surface  phenomena  and  crack  propagation  of  five  hollow  propellant 
cylinders  was  photographed  while  the  cylinders  were  being  subjected  to 
internal  pressurization  at  a  rate  of  approximately  100  psi/minute. 

Three  of  the  specimens  bad  machined  outer  surfaces,  and  two  had  cast 
■  uter  surfaces.  The  behavior  of  the  machined  specimens  was  consistent 
and  differed  from  that  of  the  cast  specimens.  Clear  delineation  of 
this  difference  was  obtained  by  15  mm  photography  at  a  speed  of  400 
f 1  ames/second. 

Altbough  all  failures  were  centralized  longitudinal  cracks,  the 
machined  outer  surfaces  obviously  experienced  continual  structural 
degradation  during  a  substantial  period  before  failure.  Widespread 
blanching  and  localized  dewetting  are  visible  in  the  sequence  shown  in 
Fig.  39;  the  internal  failure  which  led  to  the  final  fracture  occurred 
at  the  interface  of  the  larger  oxidizer  particles.  An  example  of  a 
uear  surface  400+  micron  diameter  particle  initiating  fracture  of 
machined  cylinder  004  is  shown  in  detail  in  Fig.  40  .  The  extent  of 
the  surface  degradation  is  indicated  by  the  rapidity  of  the  crack 
propagation — less  than  0.0025  second,  or  one  camera  frame. 
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TABLE  XVII.  PROPELLANT  HOLLOW  CYLINDER  DESCRIPTIONS 
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12.5  Msec  Before  Failure  7.5  Msec  After  Failure 

Figure  39.  Cylinder  004:  Machined  Surface  with 
0.1-inch  Grid,  l6nm  Motion  Photo¬ 
graphy  at  400  Frames/sec 
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a.  Machined  Cylinder  004 


b.  Cast  Cylinder  C-002 
Figure  40  .  Fracture  Detail  with  0.1-Inch  Grid 
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Cast  surfaces,  however,  produce  distinctly  different  failure  phe- 
i.  No  surface  activity  was  observed  before  crack  initiation. 

r?*f oration  of  the  cylinder  wall  led  to  relatively  slow  (0,100  second) 
cra.;«  propagation,  and  extensive  analysis  of  the  various  films  indicated 
the  propellant  tearing  on  the  innerbore  proceeded  to  the  outside  surface 
where  the  crack  propagation  was  apparently  retarded  by  the  binder-rich 
outer  cast  surface.  This  indicated  behavior  is  shown  by  film  sequences 
in  Fig.  41. 

A  comparison  of  all  fracture  areas  for  the  100  psi/minute  pressur¬ 
ization  rate  (Fig.  42  )  also  indicates  another  possible  condition  which 
the  high-speed  camera  coverage  implies  but  cannot  confirm  because  of 
film  grain  limitations.  Cylinders  002  and  003  display  fractures  that 
apparently  pass  through  oxidizer  crystal  boundaries  in  contrast  to 
failure  zones  that  bypass  or  dewet  particles.  Consequently,  it  is  not 
unreasonable  to  expect  highly  localized  circumferential  strains  approach¬ 
ing  30 $  for  both  types  of  propellant  surfaces. 

Results  from  the  second  testing  pl-.se,  Detailed  Descriptions  of 
Local  Behavior,  clarify  this  particular  phenomena.  Two  cylinders,  one 
machined  and  one  with  a  cast  surface,  were  pressurized  at  a  relatively 
slow  rate  at  ambient  conditions.  A  small  portion  of  the  surface  was 
photographed  using  direct  magnification  of  approximately  20  power.  This 
technique  provided  a  viewing  area  of  0.05  sq  in.  which  included  a  repre¬ 
sentative  quantity  and  distribution  of  oxidizer  particles. 

The  machined  surface,  Fig.  43  and  44  under  comparable  conditions 
suggests  that  a  substantial  portion  of  the  exposed  oxidizer  crystals 
are  actually  damaged  during  machining.  Consequently,  cleavages  occur 
within  the  crystal  boundaries  at  relatively  low  stress  levels.  These 
vacancies  grow  with  increasing  stress  application  and  for  the  area 
observed,  this  form  of  stress  relief  prevails  over  "conventional" 
dewetting. 

Observation  of  the  cast  surface  (Fig.  45)  during  pressurization  to 
failure  (Fig.  46)is  not  particularly  revealing.  The  binder  rich  "skin" 
precludes  the  possibility  of  directly  viewing  matrix-filler  phenomena, 
although  a  number  of  distinguishing  surface  characteristics  are  avail¬ 
able  for  use  in  measuring  local  surface  strains. 

Data  from  these  tests  are  summarized  in  Table  XVIII. 
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2.5  Msec  Before  Failure  5  Msec  After  Failure  1C5  Msec  After  Failure 

Figure  42.  Cylinder  C-002:  Cast  Surface  with  0.1-Inch  Grid 
16mm  Motion  Photography  at  400  Frames/sec 


VIEWING  AREA 
0.2  IN.  x  0.25  IN. 


4 


Figure  43.  Machined  Surface  Crack  Location, 
Low  Stress  Rate,  Cylinder  005 


68.2  psi 

Machined  Surface  Detail  at  High  Magnification 
Low  Stress  Rate,  Cylinder  005 


* 


-S'  #* 


0  psi 


63.O  psi  18. 7  sec  Before  Failure 

Figure  45.  Cast  Surface  Detail  Cylinder  C-003  High 
Magnification,  Low  Stress  Rate 
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TABLE  XVIII.  SUMMARY  OF  TEST  RESULTS 
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APPENDIX  III 

SIMULTANEOUS  THERMAL  AND  MECHANICAL  STRAINING 

Two  tests  in  which  biaxial  specimens  were  stretched  and  cooled  very 
slowly  were  run  on  the  Cumulative  Damage  test  machine  to  simulate  the 
inner  bore  condition  of  a  strain  cylinder  during  "cooldown"  in  an  attempt 
to  induce  failure.  During  the  first  test,  the  specimen  was  cooled  at  a 
rate  of  3  P/hour  and  at  a  straining  rate  of  0.00143  in. /in. /hr.  Data 
collected  were:  load,  displacement  across  the  specimen  (inside  the  en¬ 
vironmental  box),  displacement  of  the  ram  (outside  the  box),  and  the 
temperature  of  the  specimen  and  air  inside  the  box. 

Figure  47a  is  a  plot  of  the  resuHs  obtained  in  the  lirst  test.  The 
test  progressed  as  expected  for  6l  hours  at  which  time  the  load  suddenly 
dropped  from  530  to  215  pounds.  A  careful  check  after  the  test  was  con¬ 
cluded  revealed  no  deficiencies  in  the  system,  and  the  specimen  itself 
was  not  damaged.  It  was,  therefore,  concluded  that  the  clamping  arrange¬ 
ment  loosened  as  the  temperature  dropped,  permitting  the  grip  on  the 
specimen  to  slip  until  the  clamping  bolts  bottomed  out  with  the  holes  in 
the  end  bar.  When  a  positive  mechanical  linkage  had  been  re-established, 
the  load  started  increasing  again.  The  test  was  allowed  to  run  until 
the  ramp  generator  completed  a  full-scale  run  (87  hours)  despite  the  doubt 
that  the  system  was  working  satisfactorily. 

The  true  strain  value  cannot  be  precisely  determined  because  of  the 
slippage.  However,  it  possible  to  estimate  the  amount  of  slippage  by 
simply  proportioning  the  loads  at  the  time  of  slippage  since  some  of  the 
tensile  load  remained  after  the  slippage  occurred.  Thus,  correcting  the 
final  displacement  figure,  the  final  strain  is  estimated  to  be  8. 8/6. 

The  peak  load  is  estimated  to  have  been  1450  pounds  as  indicated  by  the 
dotted  line,  Fig.  47a.  The  solid  line  connects  the  actual  data  points 
collected.  The  peel;  load  was  not  recorded  because  it  was  felt  that  very 
little  relaxation  would  occur  in  the  specimen  after  the  test  stopped. 

Since  the  test  terminated  in  the  middle  of  the  night,  it  was  assumed  that 
the  reading  the  next  morning  would  approximate  that  at  the  end  of  the 
ramp  run  since  the  machine  would  hold  the  final  strain  and  temperature. 

It  is  clear  that  considerable  relaxation  occurred  at  -100  F  as  indicated 
by  the  continued  drop  in  the  load.  Maximum  stress  at  87  hours  was  308 
psi.  The  stress  after  30  additional  hours  at  constant  strain  and  tempera¬ 
ture  had  dropped  to  133  psi.  The  test  was  discontinued  without  failure 
occurring. 
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In  the  second  test,  a  higher  strain  rate  was  utilized,  but  the 
cooling  rate  was  held  at  the  same  value.  The  strain  rate  in  this  case 
was  0.0029  in./in. /hr,  and  the  rate  in  temperature  change  was  3  F/hr. 

A  ponitive  pull  was  assured  by  removing  all  of  the  play  between  the 
clamping  bolts  and  the  holes  in  the  end  bar  before  starting  the  test. 
Test  results  are  plotted  in  Fig.  47b.  Failure  occurred  after  62.5  hours 
at  a  stress  of  170  psi  when  the  temperature  was  at  -37  F.  The  strain 
at  the  time  of  failure  was  18.2JG. 

These  data  are  quite  encouraging  since  precise  failure  data  can  now 
be  related  to  both  temperature  and  mechanical  strain  inputs  for  the 
critical  biaxial  stress  condition.  Thus,  the  critical  limits  for  these 
changes  as  they  affect  both  mechanical  properties  and  failure  can  be 
defined  in  terms  of  realistic  motor  failures. 


Figure  47-  Data  Traces  Taken  During  Simultaneous  Straining 
and  Cooling  Tests  of  Biaxial  Strip  Specimens 
13-5-  x  3/8-  x  2.8-inch  Gage  Length). 


Figure  47*  Concluded 
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APPENDIX  IV 

MACHINE  COMPLIANCE  AND  POKER  CHIP  TESTS 


To  better  understand  conditions  existing  at  the  boundaries  of  a  poker 
chip  specimen  during  a  pull  on  the  Instron  tester,  it  was  necessary  to 
determine  the  load-vs-deformation  curve,  or  compliance,  of  the  machine  and 
the  related  poker  chip  apparatus  as  one  unit.  These  tests  were  conducted 
with  the  same  hardware  (instron  tester  and  related  poker  chip  attachment) 
as  was  used  in  the  Cumulative  Damage  poker  chip  test  program.  The  only 
difference  was  that  instead  of  being  bonded  between  two  platens,  propel¬ 
lant  specimens  were  bolted  together  with  three  0.5-inch  bolts.  Two  dif¬ 
ferent  Instron  load  cells  (5  and  10K)  were  used  in  the  tests.  All  tests 
were  conducted  at  77  F  at  0.2,  0.2,  and  20.0  in. /minute  crosBhead  speed. 
Figure  48'  shows  related  hardware  in  position  on  the  Model  TTD  Instron, 

S/N  L572C  Figure  48  also  shows  the  position  from  which  the  crosshead 
def lecti^tf  was  measured. 

Tables  XIX  through  XXIV  summarize  data  obtained  at  the  three  condi¬ 
tions  of  crosshead  speed  for  the  two  load  cells.  Table  XXV  summarizes  data 
from  a  test  in  which  all  the  poker  chip  test  apparatus  was  removed  and  a 
direct  connection  was  made  between  the  Instron  10K  cell  and  the  crosshead 
(via  a  1.5-inch  length  of  1.25-inch  steel  bar  stock).  This  last  test  was 
conducted  to  ascertain  the  maximum  stiffness  obtainable  with  the  Instron 
machine  alone. 

Figures  49  through  51  are  plots  of  the  data  given  in  Tables  XIX 
through  XXIV.  In  Fig.  49  and  50,  the  relationship  between  load  and  de¬ 
formation  appears  to  be  independent  of  crosshead  speed  to  the  extent  of 
reproducing  the  hysteresis  loop. 

The  curve  for  the  machine  only,  Fig.  51,  is  based  on  the  0.02-in./ 
min  rate  data  only.  However,  tests  were  also  conducted  at  0.2  in. /min; 
but  since  all  previous  data  indicated  no  effect  from  rate  changes,  these 
data  were  not  reduced. 

The  load  deformation  curves  differ  considerably  from  those  expected. 

It  was  expected  that  the  stiffness  of  the  machine  alone  would  be  above 
1.0  x  10  lb/in.  The  value  obtained  is  approximately  110,000  lb /in,  over 
the  0-  to  3000-pound  range.  This  means  that  the  poker  chip  specimen  has 
a  much  higher  stiffness  than  the  machine  (3  x  10°  compared  to  1  x  10^  lb/ 
in.).  Therefore,  one  can  postulate  that  a  constant  rate  of  displacement 
at  the  boundary  is  not  being  obtained. 

Displacement  time  curves  for  the  machine  crosshead  can  be  obtained  by 
plotting  data  from  Tables  XIX  through  XXV.  In  general,  these  plots  show  a 
reasonably  repeatable  relationship  except  at  20.0  in. /min  crosshead  speeds. 

This  effort  was  sponsored  by  Rocketdyne  as  a  portion  of  a  larger  in- 
house  program.  It  is  reported  herein  because  of  the  pertinence  of  the 
data  in  support  of  the  triaxial  poker  chip  experimental  tests  conducted 
during  the  course  of  the  Cumulative  Damage  program. 
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TABLE  XIX.  SUMMARY  OF  BATA  FOR  POKER 
CHIP  COMPLIANCE  TEST* 


Instron 

Crosshead 

Time , 

Load , 

Deflection 

sec 

lb 

(LVDT) ,  in. 

30 

323 

0.0085^ 

60 

733 

0.0173 

CO 

O 

90 

1260 

0.0265 

bC 

120 

1865 

0.0360 

>-s 
f  CD 

150 

2530 

0.0440 

CO 

0) 

u 

180 

3270 

0.0524 

d 

M 

210 

4020 

O.O6O5J 

240 

4720 

0.0685^ 

be 

270 

4700 

0.0685 

d 

•H 

300 

4690 

0.0685 

'  0 

X 

330 

4000 

0.0625 

0.0544^ 

360 

3122 

390 

2370 

0.0460 

420 

1685 

0.0373 

a 

450 

1088 

0.0286 

'  -P 

480 

590 

0.019 

£ 

510 

212 

0.0098 

540 

0 

0  J 

*Test  Conditions:  5K  load  cell 
Instron  TTD  S/N  1572.  Crosshead 
0.02  in. /min  at  77  F  with  poker 
chip  apparatus  (see  Fig.  48) 


TABLE  XX.  SUMMARY  OF  DATA  FOR  POKER 
CHIP  COMPLIANCE  TEST* 


Instron 

Crosshead 

Time, 

Load, 

Deflection 

sec 

lb 

(LVDT),  in. 

3.2 

402 

0.010^ 

6.2 

826 

0.019  « 

9.2 

12.2 

15.2 

18.2 

21.2 

23.54 

27.2 

30.2 

33.2 

36.2 
9.2 

42.2 

45.2 

48.2 

51.2 

54.2 

57.2 

59.2 


1355 

1990 

2680 

3390 

4100 

4690 

4685 

4660 

4650 

4500 

3550 

2760 

2030 

1395 

835 

410 

102.2 

0 


0.0282  *"• 
9 

0.0372 1  .5 
0.0455  f  S 

u 

0.0538  a 

M 

0.0610 

0.0673 

0.0673^  u> 

0.0673  >| 

0.0673J  1 

0.0665^ 

0.0581 

0.0500 

0.0420  u 

>h 

0.0333  p 

0.0235  ^ 
0.0150 
0.0060 
0 


*Test  Conditions:  5K  load  celj 
Instron  Model  TTD  S/N  1572. 
Crossherd  speed  0.2  in. /min  at 
77  F  with  poker  chip  apparatus 
(see  Fig.  48  ) 
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TABLE  xn.  SUMMARY  OF  DATA  FOR  POKER 
CHIP  COMPLIANCE  TEST* 


Time, 

sec 

Instron 

Load, 

lb 

Crosshead 
Deflection 
(LVDT),  in. 

0.02 

71 

0.004 

0.04 

213 

0.0075 

0.06 

406 

0.0122 

0.08 

'84 

0.0184 

a 

0.10 

1070 

0.C255 

3 

0.12 

1588 

0.033 

{  c 

•H 

0.14 

2200 

0.040 

a 

oi 

0.16 

2860 

0.047 

0 

c 

W 

0.18 

3640 

0.056 

0.20 

4440 

0.063  > 

0.02 

4110 

0.071  "3 

0.04 

3800 

0.0672 

0.06 

3390 

0.0622 

0.08 

2890 

0.0575 

0.10 

2350 

0.0505 

c 

L 

0.12 

1770 

0.0428 

(  +* 

01 

Oh 

0.14 

1205 

0.0355 

0.16 

695 

0.0266 

0.18 

290 

0.0172 

0.20 

0 

0  J 

*Test  Conditions:  5K  load  cell 
Ins iron  Model  TTD  S/N  1572.  Cross¬ 
head  speed  20.0  in. /min  at  77  F 
with  poker  chip  apparatus  (see 
Fig.  AS  ) 
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TABLE  XXII.  SUMMARY  OF  DATA  FOR  POKER 
CHIP  COMPLIANCE  TEST* 


Instron 

Crosshead 

Time , 

Load, 

Deflection 

sec 

lb 

HSB 

in. 

30 

SSSflffflf^B 

60 

1 

90 

120 

1 

0.0240 

150 

0.030 

O 

180 

2640 

0.0355 

«> 

1.2 

210 

3400 

0.0405 

CO 

240 

4000 

0.0460 

u 

0 

a 

270 

4620 

0.0522 

M 

300 

3200 

0.0575 

330 

5820 

0.0635 

360 

6260 

O.O695J 

30 

5420 

0.0635> 

60 

4660 

0.0575 

90 

3980 

0.052 

120 

3320 

0.045 

150 

2700 

0.040 

a 

180 

2110 

0.035 

U 

>3 

r 

210 

1560 

0.029 

240 

1080 

0.0254 

270 

607 

0.0172 

300 

264 

0.0105 

330 

26.4 

0.005 

350 

0 

0  > 

*Test  Conditions:  10K  load  cell 
Instron  Model  TTD  S/N  1572.  Cross¬ 
head  speed  0.02  in. /rain  at  77  F 
with  poker  chip  apparatus  (see 
Fig.  48) 
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TABLE  XXIII.  SUMMARY  OF  DATA  FOR  POKER 
CHIP  COMPLIANCE  TEST* 


Instron 

Crosshead  ! 

Time, 

Load , 

Deflection  1 

sec 

lb 

(LTOT),t 

in. 

1.0 

65 

0.0048 

4.0 

606 

0.0145 

7.0 

1440 

0.0245 

0 

10.0 

2320 

0.0335 

60 

k.2 

13.0 

3190 

0.0415 

/  CO 

CO 

16.0 

4300 

0.0505 

u 

0 

a 

19.0 

5330 

0.060 

M 

22.0 

6250 

0.069 

0.0588*' 

3.0 

4830 

6.0 

3722 

0.0485 

9.0 

2700 

0.0400 

12.0 

1690 

0.032 

a 

u 

\  3 

15.0 

1002 

0.0228 

/+» 

£ 

18.0 

344 

0.0128 

21.0 

0 

0.002 

22.0 

0 

0 

> 

*Test  Conditions:  10K  load  cell 
Instron  Model  TTP  S/N  1572.  Cross¬ 
head  speed  0.2  in. /min  at  77  F 
with  poker  chip  apparatus  (see 
Fig.  48) 
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TABLE  XXIV.  SUMMARY  OF  DATA  FOR  POKER 
CHIP  COMPLIANCE  TEST* 


*Test  Conditions:  10K  load  cell 
Ir.stron  TTD  S/N  1572.  Crosshead 
speed  20.0  in. /min  at  77  F  with 
poker  chip  apparatus  (see  Fig. 
48) 


TABLE  XXV.  SUMMARY  OF  DATA  FOR  INSTRON 
ALONE  COMPLIANCE  TEST* 


Time, 

sec 

Instron 

Load, 

lb 

Crosshead 
Deflection 
(LVDT),  in 

30 

498 

0.0088 

60 

1295 

0.0172 

90 

2460 

0.0260 

120 

3810 

0.0350 

150 

5200 

0.0420 

180 

6520 

0.0500 

210 

8040 

0.0582 

216 

8320 

0.0600 

0 

8080 

0.0605 

30 

5520 

0.0522 

60 

3380 

0.044 

90 

1780 

0.035 

120 

615 

0.0265 

150 

50 

0.0175 

180 

0 

0.00850 

*Test  Conditions:  10K  load  cell 
Instron  TTD  S/N  1572.  Crosshead 
speed  0.02  in. /min  at  77  F.  Direct 
coupling  between  10K  load  cell  and 
crosshead  ( 1.5-inch  length  of  1.25- 
inch  steel  rod) 


Instron  Load,  lb 


0  10  20  30  40  50  60  70  80  90  100 


_3 

Crosshead  Deflection,  in,  x  10 

Figure  49.  Compliance  Curve  (Load  vs  Deflection) 
for  Poker  Chip  Apparatus  with  5K 
Instron  Load  Cell 


160 


Instron  Load,  lb 
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Figure  50.  Compliance  Curve  (Load  vs  Deflection) 
for  Poker  Chip  Apparatus  with  10K 
Instron  Load  Cell 
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Figure  51  .  Compliance  Curve  (Load  vs  Deflection)  for 
Instron  Model  TTD  S/N  1572  with  10K  Load 
Cell 
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APPENDIX  V 


MATERIAL  PROPERTIES 


This  appendix  is  included  to  familiarize  the  reader  with  the  sequence 
of  events  and  methods  used  to  obtain  the  material  property  results  pre¬ 
sented  in  this  report. 

Uniaxial  experimental  relaxation  data  was  furnished  to  Rocketdyne  by 
Lockheed  Propulsion  Company  (12);  and  limited  biaxial  relaxation  data,  also 
furnished  by  LPC,  was  shifted  at  Rocketdyne  to  determine  a  temperature¬ 
time  equivalence  relationship. 

Using  the  experimental  uniaxial  and  biaxial  relaxation  modulus  curves, 
the  creep  compliance  curves  were  computed  from  the  relation  (32) 


d  (tj  e  ,(t)  „  lis-sn 

crp  rel  niT 


(139) 


where  n  is  the  logarithimic  slope  of  the  relaxation  modulus  curve. 

Since,  the  theory  of  fracture  mechanics,  successive  integrals  of  the 
creep  compliance  and  relaxation  modulus  curves  are  necessary  to  compute 
Y,  the  characteristic  energy  release  rate  constant  for  different  kinds  of 
loading  histories,  a  short  discussion  of  the  method  used  to  obtain  them 
will  be  in  order. 

When  the  creep  compliance  curve  is  obtained  from  Eq.  l‘"9-  a  Prony 
series  is  written  for  the  curve  which  is  typically  of  the  f.)*m 


D  (t)  *  D  -  X  D  e-t/Ti 
crpv  e  .  1 


(140) 


This  gives  an  approximation  to  the  curve  which  can  be  easily  integrated 
term-by-term.  The  notation  that  will  be  used  for  the  first  time  integral 
of  the  creep  compliance  is 


*  /  °  D  (t) 

v  o'  Jo  crp'  7 


the  second  integral  is 


D^(t  ) 
v  o' 


■  fo  °  lo  Dcrp(a)  dCt  de 


(141) 


(142) 
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and  so  on.  Likewise,  a  Prony  series  for  the  relaxation  modulus  is  of 

the  form 

E  .  (t)  >  E  *  EE  e_t/Ti 
r*l  '  '  e  j  1 

and  the  integral  will  be  defined  in  the  same  way  as  before,  i.e., 


(143) 


e(1)<4„>  ■  f°  ^ 


3(t)  dr 


(1M) 


E  (a)  da  dB 


This  same  reasoning  can  be  extended  to  the  nth  integrals  defined  as 


Dw(t) 


-/0tD(n-l)( 


(145) 


E^t)  r  fo  E<“  ■  ^(a)  da 


The  results  of  some  integrations  of  the  Prony  series  are  as  follows 


E^(t  )  •  E  t  -EtE  L-to/Ti  -  lj 

0  e  °  i  1  1  L  J 

*(2)(t0)  -  V  *  ?  'S  ['’t0/Tl ->]^0<  TiEi 


(146) 


(147) 


^  -  *e  *4 *  <-  *  r w  [e_to/Ti  -  Vm ~r;  1C  =  > £ <-  (us) 


»(i)(t.)- v.  +  rv»  [e'Vi-1]  (149) 

»(2)(V  ■  »e  T-  -  t  Vi*  [•>  to/Tl  -  l]  -  t0  s  D1T1  (150) 

I>(n)(t  )  .  De  C  *  (-  l)n  *  1  I  D  t*  fe-to/Tl  -  ll  .  n  i  1  f  I  if-  (-  ])"-»  0^(151) 

0  nl  1  t  J  a  .  )  1  .  J 


(150) 


.XT'*-  V 


*  •***,' + 


The  power  law  for  the  relaxation  modulus  is 


Knl^  ”  (1  +  t/T„)n  +  KP 


(152) 


where 


Eg  =  glassy  modulus 
Ee  ■  rubbery  modulus 
tq  =  that  time  when  Ere](t)  =  Eg/2C 
n  =  the  slope  of  the  log  curve 

The  power  law  constants  for  uniaxial  and  biaxial  relaxation  formulas 
follow  in  Table  XXVI. 

TABLE  XXVI.  RELAXATION  FORMULA  CONSTANTS 


Eg,  psi 

E  ,  psi 
e 

n,  num. 

t  ,  sec 

0 

Uniaxial 

200,000 

100 

-0.194 

1.44 
x  10-15 

Biaxial 

470,000 

170 

-0.204 

5-13  i. 
x  10-15 

The  Prony  series  coefficients  are  given  in  Table  XXVII  for  Ere^(t)  and 
Dcrp(t)  for  both  uniaxial  and  biaxial  curves. 

The  plots  of  E^1),  e(2),  E^),  e(4),  d(1),  D^2) ,  d(5),  and  d(4)  for 
both  uniaxial  and  biaxial  cases  are  included  in  Fig.  52  through  55. 

They  have  been  normalized  by  multiplying  each  EW  or  DW  by  Si.  so  that 

they  can  be  plotted  on  the  same  graph.  In  addition,  the  shift  factor  (aj) 
for  both  uniaxial  and  biaxial  data  are  included  in  Fig.  56  and  57 > 
respectively. 

Based  on  the  biaxial  relaxation  modulus,  the  storage  (E1 )  and  loss 
(E")  moduli  were  derived  and  are  plotted  vs  frequency  in  Fig.  58. 

Other  physical  properties  furnished  by  Lockheed  pertinent  to  STV 
propellant  include: 

Glass  Transition  Temperature  -82  C  (-116  F) 

Coefficient  of  Thermal  Expansion  5*27  +  0.28  x  10~5  in. /in. /deg  F 
Density  0.06)3  +  0.00184  lb/cu  in. 


TABLE  XXVII.  FRONY  SERIES  COEFFICIENTS 


Uniaxial 

■r. 


Biaxial 


O.OOQs 


102 

19 

103 

0.001355 

103 

40 

!04 

0.0006 

101 

55. 

102 

0.00162 

102 

110 

IO3 

0.0008 

10° 

125 

101 

0.000946 

io1 

110 

io2*2 

0.0007 

10"1 

185 

10° 

0.000757 

10° 

273 

io1,2 

0.0006 

io"2 

320 

10-°- 8 

0.000488 

io"1  . 

400 

io0*2 

0.0007 

10-3 

470 

IO-1'6 

0.0005 

io"2 

900 

io-°-8 

0.0004 

o 

1 

it- 

790 

io"2,7 

0.0004 

10"3 

1,200 

io"1,7 

0.00035 

io-5 

...  1.100 

10-3- 7 

0.00023 

10~4 

2,400 

io-2-7 

0.00015 

io-6 

1.700 

io-4'7 

0.000159 

IO"5 

2,700 

10-3- 7 

0.00012 

O 

1 

•vj 

3,000 

10-5-7 

0.00012 

10"6 

5,800 

IO"4-7 

0.00006 

IO'8 

4,400 

0.000065 

io"7 

8,000 

10-5-7 

0.000046 

io"9 

6.700 

107’7 

0.000036 

10-8 

10,800 

IO’6-2 

0.000028 

io-10 

10,400 

10-8.S 

0.00003 

io"9 

19,000 

IO’7-7 

0.000019 

io"11 

17,300 

io'9-7 

0.000015 

10-10 

34,000 

io-8-7 

0.000009 

io-12 

27,000 

io"10-7 

0.000011 

io"11 

50,000 

IO-9,7 

0.000006 

h-» 

O 

1 

h-» 

40,000 

io-n.8 

0.000009 

io"12 

70,000 

10~10-7 

0.0000046 

h-> 

© 

1 

•t* 

40,000 

io-12,7 

0.000006 

io"13 

110,000 

10-n.7 

0.0000028 

io"15 

23,000 

l0-!3.8 

0.000001 

io-14 

110,000 

io-12-7 

0.0000015 

IT15- 

18,000 

IO"19 

0.0000012 

IO"13 

25,000 

io-13-9 

0.0000006 

io"17 

5,348 

0.0000006 

IO"14-8 

0.0000002 

io"17 

0.0000002 

IO’16 

0.0000001 
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Figure  53.  Biaxial  Stress  Relaxation  Modulus 
and  the  First  Four  Time  Integrals 
vs  Reduced  Time 


Figure  55 .  Biaxial  Creep  Compliance  and  the  First 
Four  Time  Integrals  vs  Reduced  Time 


Shift  Factor  Log 


Figure  56.  Uniaxial  Shift  Factors  for  Lockheed 
STV  Propellant 
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APPENDIX  VI 

CUMULATIVE  DAMAGE  TESTING 

The  unique  test  apparatus  furnished  by  Rocketdyne  was  designed  and 
built  by  MTS  Systems,  Minneapolis,  Minnesota,  specifically  to  Rocketdyne 
specifications  for  use  on  the  AFRPL  Cumulative  Damage  program.  This 
machine  is  considered  unique  because  it  provides  for  the  first  time  the 
means  of  determining  fatigue  life  for  strip  biaxial  specimens  of  solid 
propellant.  With  the  MTS  tester,  fracture  data  can  be  generated  over  a 
wide  range  of  conditions  including  various  combinations  of  displacements, 
frequency,  and  temperature. 

The  tester  consists  of  three  basic  units:  load  frame,  control  con¬ 
sole,  and  hydraulic  power  unit.  Figure  59  shows  the  control  console  and 
the  load  frame  with  a  biaxial  test  specimen  mounted  ready  for  testing. 


Figure  59 .  Cumulative  Damage  Tester 
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The  large  metal  base  supporting  the  load  frame  is  not  part  of  the 
original  design.  Originally,  the  load  frame  was  to  be  bolted  to  a  seismic 
mass  of  about  11,000  pounds  (a  cube  of  concrete  extending  below  the  labo¬ 
ratory  floor  lead).  However,  excessive  frame  movements  were  occurring  at 
50  cps  and  other  certain  critical  frequencies.  Excessive  frame  movements 
induced  inertia  loadings  which  generated  false  output  from  the  force  load 
cells.  This  made  it  impossible  to  interpret  the  force  data  at  higher  fre¬ 
quencies.  The  problem  was  eliminated  by  suspending  the  large  supporting 
frame  on  the  four  air  bellows.  These  bellows  exert  enough  force  upward 
to  lift  the  load  frame  as  well  as  the  original  seismic  mass.  Thus,  the 
load  frame  and  the  cube  of  concrete  are  floating  on  air.  This  technique 
very  effectively  eliminated  the  high-frequency  vibration  problem. 

Driving  force  for  loading  the  specimen  is  derived  from  a  double- 
acting,  hydraulically  operated  actuator  which  receives  pressurized  hydrau¬ 
lic  fluid  through  an  electronically  controlled  servovalve  that  meters 
fluid  to  either  end  of  the  piston,  Fig.  60.  Servovalve  position  is  con¬ 
trolled  by  an  electronic  circuit  that  puts  out  a  command  depending  on 
the  error  detected  between  a  desired  voltage  function  and  the  output 
from  one  of  three  transducers,  which  is  following  a  physical  change  in 
the  test  specimen.  Through  this  "closed  loop"  arrangement  it  is  possible 
to  make  the  actuator  or  driver  follow  any  desired  movement  within  the 
isechanical  and  hydraulic  limitations  of  the  driver  itself.  Thus,  it  is 
possible  to  vary  the  displacement  of  the  system  independently  of  the  load 
in  the  specimen  and,  conversely,  to  program  a  load  time  function  indepen¬ 
dently  of  the  displacement  or  strain  in  the  specimen.  These  operations 
can  be  conducted  to  control  specified  displacement  strain  and  load  func¬ 
tions  on  the  specimen  while  separate  static  transient  or  cyclic  tempera¬ 
ture  environments  are  induced.  The  temperature  control  point  can  be 
either  the  specimen  temperature  or  the  ambient  air  temperature  in  the  box. 

Characteristics  of  the  system  are 


Max  tensile  or  compressive,  4500  lb 
Max  stroke,  2.0  in. 

Max  ram  speed,  2600  in ./min 
Min  ram  speed,  4  x  10”5  in ./min 

Environmental  Box 

Max  temperature,  350  F 
Min  temperature,  -100  F 
Rate  of  change  in  temperature  (wide 
range  of  rates  possible  using  ramp 
generator  as  command) 


Electronic  Wave  Form  Inputs 
Sine 
Square 
Triangle 

Ramp  (Linear,  bilinear) 

Maximum  Ram  Stroke  at  Sine 
Frequency 

0.010,  150  Hz 
0.250,  50  Hz 
2.0,  1  Hz 
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Figure  60,  Sketch  Showing  Operating 
Principle  of  "Closed 
Loop"  MTS  Test  Equipment 
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APPENDIX  VII 


SINUSOIDAL  FORCE  LOADING  OF 
BIAXIAL  STRIP  SPECIMENS 

This  appendix  presents  experimental  failure  data  obtained  from 
biaxial  strip  specimen  tests.  The  boundaries  of  the  specimens  were 
subjected  to  a  varying  tensile  displacement  chosen  to  force  the  total 
load  to  vary  sinusoidally  with  time. 

The  13.25-  x  0.375-  x  2,80-inch  specimen  was  placed  between  two 
aluminum  platens  and  bonded  to  each  surface.  The  force  through  the 
machine  linkage  was  monitored  by  two  transducers  which  supplied  a  sig¬ 
nal  to  a  com  land  unit;  this  signal  was  then  compared  to  an  incoming 
signal  from  a  sinusoidal  function  generator.  The  voltage  resulting 
from  the  diffeience  in  these  two  signals  was  amplified  and  operated  an 
electromechanical  servo-valve  which  controlled  oil  flowing  from  a  res¬ 
ervoir  to  eitner  side  of  a  piston  which  was  moved  by  the  oil  pressure. 

The  harmonically  varying  force  boundary  was  applied  at  frequencies 
of  5,  15,  and  45  hertz.  At  each  frequency  level,  peak-to-peak  loads  of 
230,  200,  and  180  pounds  were  input.  The  load  form  was  displayed  on  an 
oscilloscope  and  the  number  of  cycles  was  recorded  electronically.  The 
specimen  was  observed  and  when  a  crack  appeared,  the  cycle  number  and 
time  was  recorded. 

The  data  and  a  typical  sinusoidal  wave  form  are  presented  in  Fig. 
6l  and  62.  The  theoretical  waveform  would  be 

F(t)  =  Fq  +  Fq  sin  ujt  (153) 

which  would  look  like  the  following: 

F 
F 

0 

t 

What  was  actually  imposed  was 

F(t)  -  F.^  +  FQ  sin  u)t 


(154) 


» 

giving  a  curve  such  as: 


So  the  machine  control  proved  to  he  incapable  of  imposing  a  step  load 
with  a  superposed  sinusoidal  load  of  exactly  the  same  amplitude.  However, 
observed  variations,  Fj  -  Fg,  were  always  less  than  5$  of  the  peak  force. 

Figure  6l  shows  the  number  of  cycles  to  failure  following  various 
stress  levels.  Stresses  were  calculated  by  dividing  peak  forces,  Fj  +  F0, 
by  the  nominal  specimen  cross-sectional  area. 

Figure  62  shows  a  typical  force-vs-time  history  from  one  of  these 
tests.  The  curve  varies  slightly  from  the  desired  sinusoidal  shape. 
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APPENDIX  VIII 

MEASUREMENT  OF  CRACK  GROWTH 

Althougb  recording  the  rate  of  crack  growth  across  the  specimen 
following  initial  failure  was  not  a  contractual  requirement  of  the  Cumu¬ 
lative  Damage  program,  it  soon  became  evident  that  unless  this  rate  of 
crack  growth  was  recorded  much  valuable  data  would  be  lost.  To  accomplish 
this,  it  was  necessary  to  photograph  the  specimen  at  some  precise  time 
interval  at  the  point  where  the  specimen  had  reached  its  maximum  extension 
during  the  cycle.  An  existing  Hulcher  70  mm  sequence  camera  was  modified 
to  expose  a  single  frame  upon  command.  Exposure  was  accomplished  by  first 
opening  an  auxiliary  shutter  then  discharging  an  electronic  flash  unit 
(having  a  15-microsecond  duration)  whose  trigger  circuit  was  energized 
by  the  voltage  output  from  the  appropriate  MTS  transducer  (stroke  or 
force)  at  the  right  point  on  the  wave  form.  To  accomplish  these  func¬ 
tions,  it  was  necessary  to  design  and  build  a  special  sequence  timer  and 
triggering  circuit.  This  activity  was  performed  by  Engineering  Test 
personnel  and  supported  totally  by  Rocketdyne. 

It  was  not  possible  to  obtain  a  photographic  record  for  each  fatigue 
test  because  testing  had  to  progress  simultaneously  with  the  building  and 
development  of  the  sequence  timer  and  triggering  unit.  However,  19  tests 
were  photographed,  and  the  data  were  reduced  (see  Fig.  63).  The  camera 
was  normally  used  to  monitor  the  specimen  for  crack  initiation — it  was 
set  to  take  photographs  at  relatively  infrequent  intervals.  The  time 
between  photographs  was  greatly  reduced  after  a  crack  initiated.  Crack 
initiation  varied  from  specimen  to  specimen.  When  the  crack  started  on 
the  back  side  of  the  specimen  (away  from  the  viewing  window)  it  progressed 
through  the  specimen  and  appeared  abruptly  on  the  fron  side  with  an  ap¬ 
preciable  initial  length  (~0.5  inch).  Such  crack  growth  was  measured 
during  the  rapid-growth  phase  as  depicted  in  Fig.  64d  through  64h.  When 
the  crack  initiated  on  the  front  (viewing)  side  of  the  specimen,  it  was 
noted  that  the  time  required  for  the  initial  surface  crack  to  work  through 
the  thickness  of  the  specimen  was  quite  significant  (Fig.  64a  through 
64d).  Front-surface  propagations  always  progress  in  the  manner  depicted 
in  Fig.  64a  through  65g.  All  specimens  tested  were  milled  specimens  and, 
therefore,  dod  not  have  a  binder-rich  layer  as  would  be  the  case  in  full- 
scale  motors.  Mechanical  damage  incurred  from  the  milling  process  pro¬ 
duces  many  sites  where  large  perchlorate  crystals  have  been  crushed,  and 
these  apparently  influence  the  failure  initiation. 
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(Times  between  figures  are  not  uniform.  Longer  periods  occurred  between 
a,  b,  c,  and  d  than  occurred  between  e,  f,  g,  and  h.  Sequence  camera 
did  not  photograph  small  initial  cracks  between  a  and  d.) 


Top  View  Looking  at  Cross  Section 
Going  Through  Crack 
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Rapid  Growth  Starts.  Note  Growth  in 
Each  Direction  Symmetrically  Away 
from  Original  Site 
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Figure  64.  Concluded 
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Replots  for  the  crack  growth  data  are  included  for  the  19  tests 
(Fig.  65  through  83 .  Seven  of  these  replots  reflect  the  low  initial 
rate  as  depicted  in  Figure  64b  .  The  remaining  twelve  curves  do  not.  In 
most  cases  the  crack  suddenly  appeared  from  the  back  side,  and  it  was 
not  possible  to  track  the  initial  small  surface  crack  in  these  cases,  as 
mentioned  earlier. 

Average  rates  shown  on  these  replots  are  based  on  the  growth  of  the 
crack  after  it  had  worked  completely  through  the  thickness  of  the  speci¬ 
men  from  the  initial  surface  crack  and  growth  had  started  toward  the 
edges. 


Figure  65.  Biaxial  Fatigue  Test  92 
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Figure  70.  Biaxial  Fatigue  Test  97 


Biaxial  Fatigue  Test  98 


Figure  72.  Biaxial  Fatigue  Test  100 
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Fatigue  Test  102 


Dual  Strain: 


Crack  Length,  in 


Crack  Length,  in. 
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47  cps 
Dual  Strain: 

6jt  Initial  (To  10,000  Cycles) 
To  Failure 

14.4  in. /min  Average  Bate 
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Figure  78.  Biaxial  Fatigue  Test  106 
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Figure  79.  Biaxial  Fatigue  Test  107 
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Figure  80.  Biaxial  Fatigue  Test  103 
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Figure  81.  Biaxial  Fatigue  Test  110 
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Figure  83.  Biaxial  Fatigue  Test  113 
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APPENDIX  IX 

UNIAXIAL  VS  BIAXIAL  PROPERTIES 


During  the  conduct  of  the  tasks  required  on  the  Cumulative  Damage 
Program  data  have  become  available  which  provide  additional  information 
for  a  direct  comparison  of  propellant  mechanical  properties  obtained  from 
both  uniaxial  and  multiaxial  specimens.  Engineers  within  the  industry 
have  long  questioned  the  validity  of  property  data  determined  under  one 
state  of  stress  for  predicting  behavior  under  other  states  of  stress. 

But  the  expense  and  technical  difficulty  of  performing  satisfactory  multi¬ 
axial  tests  has  limited  efforis  to  define  the  magnitude  of  the  discrep¬ 
ancies.  However,  Rocketdyne*s  past  experience  with  poker  chip  tests  coupled 
with  many  recent  biaxial  tests  conducted  over  a  broad  spectrum  of  condi¬ 
tions  provided  the  data  necessary  to  make  quantitative  conclusions  con¬ 
cerning  these  discrepancies. 

The  discrepancies  between  uniaxial  and  triaxial  tests  have  previously 
been  discussed  by  fiarbert  (33,  34),  This  appendix  is  limited  mainly  to  a 
discussion  of  discrepancies  between  uniaxial  and  biaxial  conditions  as 
found  in  a  more  recent  testing  sequence. 

VIBRATION  FATIGUE  TESTS 

Vibration  fatigue  tests  were  performed  using  the  Cumulative  Damage 
Testing  machine  with  biaxial  strips  13-23  x  3/8  x  2.8  inches  long.  In  all 
cases  the  specimen  surface  temperature  was  maintained  at  a  constant  140  F 
during  the  entire  test  duration.  Fatigue  loadings  were  harmonic  (in  time) 
extensions  of  the  form 

U(t)  .  Uq  +  Uq  sin  urt  (155)  ’ 

applied  to  the  2.8-inch  gage  length 

The  resulting  forces  and  the  ambient  temperature  were  both  monitored. 
Also,  the  time  a  crack  first  becomes  visible  and  the  growth  of  that  crack 
as  a  function  of  time  (or  number  of  cycles)  was  recorded.  The  force-time 
and  force-displacement  records  were  photographic  prints  of  oscilliscope 
traces  reflecting  force  transducer  and  linear  variable  differential 
transformer  (LVDT)  outputs.  Examples  of  these  traces  are  shown  in  Fig. 

84. 


a.  Lissajous  ifettern  after 
Only  a  Few  Cycles 


c.  Lissajous  Pattern  after 
50,000  Cycles 

Figure  84. 


b.  Displacement- time  and 
Force-time  Pattern 
after  Only  a  Few  Cycles 


d.  Displacement-time  and 
Force-time  Pattern 
after  50,000  Cycles 


Oscilloscope  Traces  Resulting  from  Biaxial  Fatigue  Tests 
on  Cumulative  Damage  propellant  (Maximum  extensions 
were  0.1875  inch  applied  at  15  cps.  Specimen  surfaces 
were  held  at  140  F.) 
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The  tests  vere  conducted  to  determine  (using  the  number  of  cycles 
to  failure)  apparent  values  of  the  characteristic  strain  energy  release 
rate  and  intrinsic  flaw  size.  These  parameters,  once  established,  can 
then  be  used  to  predict  times  to  failure  following  more  complicated  load¬ 
ing  histories.  The  idealized  model  upon  which  these  predictions  are  to 
be  based  consists  of  an  initial  spherical  (or  cylindrical)  flaw  in  an 
incompressible  spherical  (or  cylindrical)  body  of  linear  viscoelastic 
material.  Analytical  expressions  based  upon  fracture  mechanics  for  the 
time  to  the  growth  of  the  flaw  in  idealized  models  of  this  type  have  been 
derived  in  the  literature  for  several  typical  loading  histories  (20,  2l). 
These  histories  include  that  defined  by  Eq.  55- 

These  tests  were  all  performed  using  a  CTPB  propellant  furnished 
for  the  Cumulative  Damage  program  studies.  An  extensive  characteriza¬ 
tion  for  this  propellant  was  provided  by  Lockheed  IVopulsion  Company 
in  support  of  the  multicontractor  Cumulative  Damage  program  (35) .  Of 
particular  importance  to  the  Rocketdyne  fracture  mechanics  approach  was 
the  relaxation  modulus  reproduced  in  Fig.  83*  These  uniaxial  data  were 
obtained  using  an  Instron  machine  at  strain  levels  of  approximately  2.5$. 
The  specimens  had  the  standard  JANAF  dimensions  but  were  tab-ended.  Dis¬ 
placements  were  quickly  imposed  (2.0  in. /min  crosshead  Bpeed)  and  then 
held  constant  while  the  Instron  load  cell  monitored  the  force  deteriora¬ 
tion.  These  measurements  were  continued  for  20  minutes;  data  collected 
before  four  loading  times  had  elapsed  were  disregarded.  Test  temperatures 
ranged  from  -120  to  300  F  to  generate  the  curve  of  Fig.  85.  The  time- 
temperature  shift  factors,  a?,  used  for  the  curve  are  shown  in  Fig.  86. 

To  apply  the  relaxation  (or  transient)  modulus  of  Fig.  85  to  the 
harmonic  displacement  input  of  the  fatigue  tests  it  became  necessary  to 
convert  these  data  to  the  dynamic  format  implied  when  one  writes 

o(t)  =  E*  e(t)  (156) 

If  the  displacement  (and  in  turn,  the  strain)  is  a  specified  harmonic 
function  in  time 

e(t)  -  e  eiu)t  (157) 

then  the  resulting  stress  may  be  written 

o(t)  .  „o.  i(urt  +  fi)  (158) 
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where  6  ia  the  angle  by  which  the  resulting  stress  leads  the  applied 
strain.  From  these  one  can  deduce  the  value  for  a  dynamic  modulus,  E* , 
by  writing 


0  +  6) 

E*  -  -2  „  V -  - 

C  coelU,t 

•o 


(159) 


If* 


0o  16 
■—  e 
e 

o 


~ ~  (cos  6  +  i  sin  6) 

v 

o 


(160) 


The  stress  was  decomposed  into  two  vectorial  components,  one  in-phaBe 
with  the  strain  and  the  other  90°  out  of  phase.  Correspondingly  the 
modulus  E*  has  an  in-phase  component,  E' ,  and  an  out-of-phase  component, 

IT ,  and  is  written  as  the  complex  number 

E*  -  E*  +  i  E"  (161) 

where  from  Eq,  160  the  real  and  imaginary  parts  are  seen  to  be,  respectively1 


E* 


~  cos  6;  E"  *  —  sin  6 
o  Co 


(162) 


Noting  the  decaying  character  of  the  relaxation  modulus,  Schapery 
(l()  has  suggested  the  use  of  a  finite  Prony  series  to  analytically  rep¬ 
resent  the  experimentally  measured  time  function.  The  series  of  terms 
is  written 


E„l(t)  -  %  .  I  Vt/Tk  (163) 

where  the  several  E^' s  and  T^' s  can  be  selected  by  collocation.  Further, 
the  term-by-terr  Laplace  transformation  of  the  functions  yields  the  image 

_  E,pT 

®r.l(p)  -  %  +  Z  OM) 


(s'2  *  p'2}  ^  -  "A 
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Stress  Relaxation  Modulus  vs  Reduced  Time 
(UOC’j  -6 IE  STV  Propellant — Furnished  by  LPC 
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Figure  86.  Time-Temperature  Shift  Fact 
to  the  Uniaxial  Relaxation 
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The  bar  implies  that  the  function  has  been  transformed.  The  inter¬ 
relationship  of  the  dynamic  properties  with  these  so-called  transient 
properties  can  be  determined  simply  by  replacing  the  transform  parameter 
p  with  +  iu).  The  function  becomes  a  complex  number,  the  real  part  of 
which  corresponds  to  E’  and  the  imaginary  part  E".  By  following  this 
line  of  reasoning  as  suggested  in  (’36),  when  p  =  +  iu) 


N  +iuu  t  E 

E(+  i«j)  =  E  +  S  Z. - 

v—  '  r  1  +  iuu  T 

n=o  ~  n 


(165) 


E(+  iuu)  =  E  +  S 


N  (cu  t  )  E 


n  n 


ip  T  E 
_ n_n 

,  x2 


n=o  (u)  t  )  +1  n=o  (w  t  )  +1 
'  n'  '  n 


4“  (166) 


=  E'  +  i  E" 


(167) 


In  polar  rotation 


E(+  iuu)  = 


E_  +  E 


N  (u)  T  )  n 


n=o  (tw  t  )  +1 
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N  10  T  E 

£  n  P 

n=o  (to  t  )  +1 
'  n' 


XT  tu  T  E 

N  n  n 

.  .  -1  ^  (uut)^+1 

+  i  tan  n=o  v  n _ 

e  N 

E  +  7!  (uu  T  )2E 

r  v  n'  n 

n=o - 5 — 

(uu  Tn)  +1 


(168) 


Comparing  Eq,  156  and  168  one  has 
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The  relaxation  modulus  of  the  Lockheed  propellant  (shown  in  Fig.  85) 
was  numerically  expanded  into  a  Prony  series  (Eq.  163)  of  23  terms  in 
addition  to  the  rubbery  modulus,  Ee.  The  best  value  for  the  coefficients 
and  are  shown  in  the  table  on  the  left  side  of  Fig.  85.  These 
coefficients  have  been  used  to  calculate  both  the  corresponding  storage 
modulus,  E',  and  loss  modulus,  E",  as  indicated  by  Eq.  166  and  168. 

These  functions  of  the  vibrating  frequency,  U),  are  shown  plotted  in 
Fig.  87- 

It  was  intended  that  these  data;  Erei  (t),  E"  (id),  and  E*  (uj);  would 
be  used  with  the  viscoelastic  fracture  mechanics  equations  for  the  growth 
of  a  spherical  or  cylindrical  flaw.  The  remaining  parameters  could  be 
chosen  to  give  the  best  fit  to  the  number  of  cycles  to  failure  observed 
during  the  fatigue  tests.  These  parameters  would  then  permit  the  use  of 
the  theory  for  predicting  failure  following  more  complicated  loading 
histories;  i.e.,  a  numerical  evaluation  of  the  cumulative  damage.  But 
first  a  cursory  comparison  of  the  properties  determined  from  uniaxial 
relaxation  tests  with  those  deduced  from  the  responses  of  the  dynamic 
testa  of  biaxial  strips  was  made.  This  step  was  performed  to  confirm 
the  validity  and  check  the  accuracy  of  the  calculations  utilizing  the 
uniaxial  relaxation  properties. 

COMPARISON 

The  test  data  used  for  these  calculations  was  that  collected  from 
t*  biaxial  strip  specimen  subjected  to  a  harmonic  tensile  displacement 
with  a  frequency  of  94.5  radian/sec.  The  equation  for  the  platen  move¬ 
ment  was 

u(t)  =  0.09375  +  0.09375  sin  94. 5t 

where  t  is  time  in  seconds  and  u(t)  the  displacement  in  inches.  This 
proved  to  be  a  trouble-free,  valid  test. 
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Figure  88  shows  replots  of  the  hysteresis  loops  shown  as  photographs 
in  Fig.  84.  But  also  shown  are  the  predicted  stresses  and  strains  corre¬ 
sponding  to  the  uniaxial  relaxation  modulus  and  the  dynamic  properties 
of  Fig.  87.  The  differences  are  seen  to  be  the  largest  at  a  low  number 
of  cycles  even  though  this  is  where  one  would  logically  believe  the 
results  should  be  best  because  the  nonlinearity  due  to  dewetting  from 
cycling  has  not  fully  come  into  play. 

Because  the  material  properties  are  a  vital  part  of  the  equations 
derived  from  the  fracture  theory,  it  was  felt  that  some  resolution  of 
this  discrepancy  should  be  accomplished  before  the  effort  continued. 

This  work  resulted  from  the  attempt  to  understand  these  phenomena  and 
intelligently  reconcile  theory  and  test  results. 

BIAXIAL  CONSTANT  STRAIN  RATE  TEST 


The  discrepancies  between  the  measured  dynamic  response  of  the 
biaxial  specimen  with  the  predicted  response  using  properties  deter¬ 
mined  from  static  tests  of  uniaxial  specimens  caused  us  to  attempt 
to  determine  just  which  new  condition  our  analytics  could  not  handle, 
the  two  dimensional  stress  state  or  the  dynamic  loading.  To  do  this 
we  ran  a  single  sawtooth  displacement  test  where  a  linear  cross-head 
rate  first  tensed  and  then  relieved  a  biaxial  specimen  (see  Fig.  89a). 
The  resulting  measured  forces  are  shown  in  Fig.  89b.  Those  forces  pre¬ 
dicted  using  the  modified  power  law  as  derived  from  the  uniaxial  relaxa¬ 
tion  data  are  also  shown.  This  four-parameter  representation  is  of  the 
form 


Erel^) 


E  -  E 

(i  f  t/ Tq  )n  +  Ee 


(170) 


and  the  procedure  for  using  it  to  predict  the  response  of  constant 
strain  rate  tests  is  discussed  in  (12).  The  results  of  this  procedure 
are  shown  by  the  curve  labeled  "Predicted"  in  Fig.  89b. 


RELAXATION  TO  CONSTANT  RATE  TESTS 


There  are  two  possible  sources  for  the  discrepancies  described  in 
the  previous  paragraphs.  One  is  that  the  relaxation  modulus  (the  ratio 
of  the  time-dependent  stress  to  strain  in  constant  strain  test)  for  the 
propellant  cannot  be  used  to  predict  constant  strain  rate  behavior  even 
under  a  like  stress  state.  The  other  is  that  the  uniaxial  properties  are 
not  appropriate  for  biaxial  or  other  stress  states.  To  check  the  im¬ 
portance  of  the  first  of  these  differences  for  the  Lockheed  propellant, 
predictions  of  the  time-dependent  force  due  to  a  linear  strain  rate  were 
made.  These  predictions  were  then  compared  to  the  forces  as  reported  by 
Lockheed  in  their  constant  rate  characterization  tests  of  uniaxial  specimens. 
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U)  = 

13*25 

u(t) 


140  F 

94.5  rad/sec 
x  3/8  x  2.8  in.  Specimen 
=  0.09375  (1+  sin  94. 5t) 


Measured 


The  predictions  were  made  using  equations  of  linear  viscoelasticity 
and  the  superposition  principle  beginning  with  the  convolution  integral 
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(171) 


which,  when  evaluated  for  a  constant  strain  rate  history  becomes 

0  -  (at)  \  f Erol  (?)  d?  (172) 

■Jo 

which  when  Rt  =  e(t)  is  substituted  becomes 
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(173) 


The  predicted  forces  are  compared  to  those  experimentally  determined  in 
Fig  90  for  strains  of  less  than  8/6.  It  will  be  noted  that  the  varia¬ 
tions  (especially  at  the  lower  rates)  are  small.  So  it  is  concluded 
that  a  significant  portion  of  the  errors  discussed  in  this  report  are 
due  to  the  state-of-stress  rather  than  the  kind  of  loading  history. 


BIAXIAL  VS  UNIAXIAL  RELAXATION  MODULUS 

As  a  consequence  of  the  comparisons  outlined  in  the  preceding  para¬ 
graphs  a  biaxial  relaxation  characterization  was  requested  of  LPC.  In 
j  response  they  furnished  the  results  of  a  series  of  biaxial  relaxation 

tests  with  the  specimens  subjected  to  an  average  strain  of  2.5J6.  Three 
j"  points  over  two  decades  of  time  at  each  of  seven  temperatures  were  re- 

[  *  ported.  The  raw  plots  are  shown  in  Fig.  91. 

I  ' 

!  To  compare  the  relaxation  modulus  indicated  by  these  tests  wi+h  that 

predicted  using  the  previously  reported  uniaxial  tests,  the  series  of 
[  curves  of  Fig.  91  were  shifted  using  factors  corresponding  to  the  curve 

;  of  Fig.  86.  The  results  are  indicated  in  Fig.  92  as  is  the  comparison 

with  the  uniaxial  prediction.  While  at  first  glance  the  prediction  looks 
good,  a  closer  inspection  shows  that  for  some  temperatures  discrepancies 
as  large  as  3 0)6  are  present.  Also  it  is  evident  that  except  for  the 
hottest  temperature  the  stiffness  prediction  based  on  uniaxial  tests 
are  lower  (more  compliant)  than  are  observed  in  the  actual  biaxial  tests. 
This  result  implies,  of  course,  that  uniaxial  data  give  unconservative 
I  predictions  of  biaxial  stresses. 
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0.0125  0.025  0.0375  0.05  0.075 

e(t)  =  cot,  in. /in. 

Figure  90.  Comparisons  of  the  Stresses,  Predicted  and  Measured, 
Resulting  from  Uniaxial  Constant  Strain  Rate  Tests — 
(Predictions  Based  on  Uniaxial  Relaxation  Data) 

CONCLUSIONS 

Analytical  predictions  of  the  stress  levels  in, realistic  stress 
fields  based  on  material  properties  deduced  from  tests  ^f  uniaxial 
specimens  give  erroneous  results.  For  simple  biaxial  strip  vibration 
and  tensile  tests  the  predicted  force  response  was  less  than  two-thirds 
of  that  measured.  The  comparison  of  biaxial  and  uniaxial  relaxation 
test  results  confirm  this  trend.  Such  inaccuracies  in  predicting  the 
responses  of  carefully  controlled  and  monitored  laboratory  tests  make 
the  prediction  of  real  motor  stress  responses  an  apparent  impossibility. 
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Figure  92.  Comparison  of  the  Biaxial  Relaxation  Data 
to  the  Relaxation  Modulus  as  Determined 
from  Uniaxial  Tests  (Biaxial  Data  has  been 
Shifted  According  to  the  a  's  of  Fig.  85  ) 
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This  variability  of  the  composite  propellant  properties  vith  stress- 
state  has  previously  been  noted  by  others.  Schapery  has  discussed  the 
propellant  nonlinearities  resulting  from  the  method  of  tests  (37)  and, 
as  noted  earlier,  Herbert  has  discussed  the  differences  he  has  discovered 
between  uniaxial  and  triaxial  moduli.  But  the  magnitude  and  importance 
of  this  effect  may  have  been  minimized  to  facilitate  arriving  at  early 
answers . 
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